General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



M83- 1783 1 


(H AS A-CH- 17071 1 ) DYNAMIC BALANCE 
IMPROVEMENT ERCGRAM Final Report 
(Rocke tdyne) 3 15 p HC A14/MF AO 1 CSCL 13K 

Duel as 

G3/37 02732 




Rockwell International 



FFR i p P3 
RECEIVED 
NASA Sn FACILITY 3/ 
ACCESS DEPT. .\y/ 



Rockwell International 

Rocketdyne Division 

6633 Canoga Avenue 
Canoga Park, California 91304 


RI/RD82-285 


DYNAMIC BALANCE IMPROVEMENT PROGRAM 
FINAL REPORT 

PHASE 1 


Contract NAS8-34423 


PREPARED BY 


M. F. Butner 
Technology Projects and 
Mechanical Elements 



APPROVED BY 


A. Csomor 
Project Engineer 
Technology Projects and 
Mechanical Elements 


&7y{./6uU^ 

T3 tl7 rr * 1 v 


F. M. Kirby 
Program Manager 
Advanced Booster Propulsion 


FOREWORD 


This final report is submitted to the National Aeronautics 
and Space Administration, Marshall Space Flight Center, in 
fulfillment of Contract NAS8-34423, Dynamic Balance Improvement. 


ABSTRACT 

A 12-month program was conducted to address the reduction of 
residual unbalance in the SSME high-pressure turbopump rotors. 
Elastic rotor response to unbalance and balancing requirements, 
multiplane and in-housing balancing, and balance-related rotor 
design considerations were assessed. Conclusions are drawn 
from the study results, and recommendations presented for 
near-term improvement of the SSME balancing, as well as future 
study and development efforts. 
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SUMMARY 


A 1-year program was conducted beginning in September 1981 to define the feasi- 
bility of improving the residual unbalance condition of the SSME high-pressure 
turbopumps. The program was conducted by the Rocketdyne Division of Rockwell 
International and Mechanical Technology Incorporated (MTI) under the direction- 
of NASA Marshall Space Flight Center, with Mr. George Von Pragenau as the Program 
Monitor. 

Digital computer programs were used in conjunction with elastic rotordynamic 
models of the high-pressure turbopumps to define the rotor balancing requirements. 
The programs were used to describe the response of the turbopump assemblies to 
residual unbalance. These response descriptions were then used to assess the 
ability of a multiplane rotor balancing computer code developed by MTI to produce 
effective correction recommendations. It was shown that if corrections are 
applied at all major masses, low-speed balancing is very effective in producing 
low response throughout the operating speed range. It also was shown that varying 
degrees of response are obtained with corrections applied at fewer planes. In 
general, more balance runs at higher speeds are required as the number of correc- 
tion planes is reduced. 

The capabilities and application of unbalance detection instrumentation were 
assessed. As a result, it was recommended that proximity measurements located 
away from the bearings be used in combination with force measurements at the 
bearing location. This combination is concluded to permit detection of response 
to both first and second modes and can be used as diagnostic tools for monitoring 
operating turbopumps for balance-related anomalies. 

It was concluded that in-housing balance corrections can be applied only to those 
planes accessible from the turbine duct or pump inlet. Access to the interior 
impellers for correction and cleaning is not considered practical. Implementation 
of balance correction by laser beam is considered premature due to the low material 
removal rate and the difficulty entailed in cleaning the pump passages. 
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Assessment of various rotor joints for balance maintenance capability indicated 
that pilots and stack fasteners that maintain or increase preloads are necessary; 
however, since these are already incorporated into SSME designs, the options for 
improving the balance maintenance consist of reductions in the number of joints or 
improvements in piloting. Means were shown for reducing the number of joints in 
the HPFTP and for improving the preburner impeller and shim piloting in the HPOTP. 

Bearing tilt was found to exert little influence on residual unbalance levels. 
Bearing wear can change response to unbalance by causing a reduction in bearing 
stiffness. 

Turbine blade erosion is likely to be distributed evenly over the blade complement 
and, therefore, does not have a strong effect on residual unbalance. 

Recommendations for improvement of the dynamic balance of the SSME high-pressure 
turbopump rotor include: 

1. Hvtljv i'erm: 

a. Initiate multiplane balancing 

b. Incorporate trim balance correction at the HPOTP preburner 
impeller 

c. Eliminate dead bands during balance operations; investigate 
hydrostatic dead-band eliminators for both balancing and 
engine operation 

d. Incorporate rotordynamic instrumentation into all high-pressure 
turbopumps 

2 . Long-Term: 

a. Study modification of SSME high-pressure turbopumps to. 
incorporate in-housing .balancing capability 

b. Incorporate in-housing balance feature."- into all future high-speed 
turbopump designs 
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INTRODUCTION 


In the SSME and in future advanced engines, turbopump perrermance and life improve- 
ments can be realized if rotor unbalance can be reduced. In present practice, the 
rotating components generally are detail balanced, then built into assemblies upon 
which a rigid rotor, two plane balance is performed. Because turbopump housings 
do not have longitudinal parting planes, the rotor assembly must be broken down to 
components prior to final turbopump assembly. As a consequence, relocation of the 
rotating component? even with interference-fit pilots produces some rotor unbalance. 
In-housing rigid rotor trim balance procedures have been used to minimize the resid- 
ual unbalance prior to actual operation. Rigid retor balancing will not achieve 
the optimum balance condition for flexible rotors that operate at speeds above one 
or more critical speeds. Correcting rotors through modal or influence coefficient 
methods at high speed couid be employed to accommodate their deflected shapes. 
However, if disassembly and reassembly are required, the potential for introducing 
unbalance through relocation still exists. The best balancing procedure involves 
high-speed, in-housing balancing that (1) treats the flexible characteristics of 
the rotor and (2) eliminates the relocation unbalance. Such a process requires 
access to the rotor for detecting unbalance and to effect corrections, as well as 
a means of driving the rotor to high speed and lubricating the bearings. In addi- 
tion, propellant cleanliness must be preserved. All these requirements are diffi- 
cult to achieve in the densely packed configuration typical of a high-pressure 
turbopump, where housing penetration may have a serious stress-increasing effect. 

The purpose of this study program was to define the requirements and methods for 
achieving satisfactory balancing of high-speed turbopumps. The studies conducted 
were focused on the SSME high-pressure fuel turbopump (HPFTP), Fig. 1, and the 
high-pressure oxidizer turbopump (HPOTP), Fig. 2. 

The elements required to achieve satisfactory balance are: 

1. Define acceptable residual unbalance magnitudes 

2. Detect unbalance 

3. Specify the correction necessary to achieve (1) 

A. Perform the correction 

5. Maintain the correction 
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To establish feasible means of Implementing these requirements for high-speed 
elastic rotors, three main technical tasks were performed: 


Task 1: b l ast lc Rotor Dynamic Requirements 

A. Adapt or develop rotordynamlc models of HP FTP and HPOTP 

B. Define balancing requirements for HPFTP and HPOTP 

C. Define measurement techniques and Instrumentation to detect unbalance 

and associated vibrations for ambient and cryogenic environments 

D. Define bearing load measurement techniques 

Task 11: In-housing Balancing Syst em 

A. Develop a multiplane flexible rotor balance program 

B. Assess the effectiveness of applying balance corrections at a limited 
number of planes to correct for unbalances distributed along the rotor 

C. Assess correction methods adaptable for In-housing balancing 

Task III: Balance Consi s tent Rotor Designs 

A. Assess conventional rotor joints for their ability to maintain a balanced 
condition 

B. Assess and develop design criteria for rotor Joint configurations that 
tighten in operation 

C. Determine the effect of bearing tilt and wear on residual unbalance 

D. Assess the unbalance caused by turbine blade erosion 

Current balancing practice for the SSME high-pressure turbopump rotors Includes 
sequential two-plane balance so that each plane Is corrected. Subsequent to the 
last two-plane balance, the rotor must be disassembled prior to assembly into the 
turbopump casing. After assembly, they are subjected to a low-speed (5000 rpm) 
tr 4 m balance. 
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The HP FTP trim balance is conducted using slave bearings. This permits the use 
of lubricants unacceptable in operational bearings. Corrections are determined 
fcr the outboard (first stage) turbine disk and the slinger at ‘.he pump inlet end. 
Final assembly includes new bearings with high points of eccentricity oriented to 
the same angle relative to the rotor as were the balance bearings. Because of 
their outboard position, the bearings can be replaced without disturbing the main 
rotor masses. 

Because of the inboard location of its bearings, the HPOTP trim balance is executed 
on its operational bearings, which are cooled during spinning with GN^. Limits are 
placed on permissible total balance spin duration to preserve bearing conditions. 
Corrections are applied to the turbine end only by exchanging turbine bolts. 

More detailed descriptions of these balance procedures are contained in Appendix D. 
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CONCLUSIONS AND RECOMMENDATIONS 


It was concluded that an Improvement in SSME high-pressure turbopump balance 
quality can be realized by instituting a multiplane balance procedure using the 
Influence coefficient method. Ideally, balance data should be taken at high 
speeds and corrections should b’ applied to inboard pump impellers, as well as 
the accessible outboard rotor components. From a practical standpoint, however, 
applying balance corrections to the interior rotor parts is not considered feasible 
in the near term. Therefore, Initiation of a multiplane out-of-housing rotor 
balance procedure is recommended for both HPFTP and HP0T1 . Trim balance of the 
HPOTI preburner impeller is recommended in addition to the present procedure. The 
derails of the required procedure are to be determined experimentally. A separate, 
but related, effort would explore the detail requirements of implementing in-housing 
balance with corrections applied to interior balance planes. 

CONCLUSIONS 

Conclusions drawn from the studies conducted under Task 1, Elastic Rotor Dynamic 
Requirements are: 

1. Trim balancing is required to obtain satisfactory response over the 
operating speed range of the HPFTP and HPOTP. 

2. A two-plane, low-speed, trim balance procedure will produce satisfactory 
response of the HPFTP. 

3. A second correction plane must be used in trin balancing the HPOTP to 
produce satisfactory response as the speed approaches 34,000 rpm. 

4. The nonlinearity introduced by bearing dead bands may result in higher 
or lower bearing loads than those predicted by linear models. The dead 
bands must be eliminated during balance operations for maximum accuracy. 

5. Hydrostatic dead-hand eliminators are superior to plunger types, as they 
maintain axial freedom of bearing races. 
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b. The high-pressure turbopumps can be instrumented to gather dynamic data 

useful for balancing during operation. Proximity measurements at stations 
away from the bearings are most useful for first mode information. Force 
measurements by stra'n guges or sensitive proximity sensors near the 
bearings are more useful for balance-related Information in the operating 
speed range near the second modes. 

Task II, In-Housing Balancing System Studies, led to the following conclusions: 

1. Practical multiplane balancing computer codes are available which are 
capable of determining necessary corrections tor the high-pressure 
turbopump. 

2. Low-speed multiplane balancing procedures, which Include corrections at 
each rotor mass, can produce satisfactory response for both HPFTP and 
HPOTP over their required operating speed ranges. 

3. As the number of balance correction planes is reduced, the number of 
speeds at which balance data must be taken increases. 

4. The most practical reduced-plane balance of the HPFTP requires correc- 
tions applied to the third-stage pump impeller and the first-stage 
turbine; balancing data must be taken at 14,000 and 29,000 rpm. 

'». The high-speed (*"34,000 rpm) response of the HPOTP can be improved if 

corrections are applied to the turbine plane and the preburner impeller. 

6. Access to interior rotor masses for correction is not currently practical 
on the SSME high-pressure turbopumps. The advantages would not outweigh 
the potential reliability degradation which would result if ports were to 
be incorporated into highly stressed casings. In addition, due to the 
multiple wall construction, a high level of complexity would be introduced 
in attaining sealing of different pressure zones in the access path. 

7. The optimum method of in-housing correction employs exchangeable weights, 
for example: 

a. Bolts or washers of selectable weight 

b. Locking set screws attachable by means of threaded holes incorporated 
into rotor components. 


O 
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Task III, the study of Balance-Consistent Rotor Designs led to the following 
conclusions: 

1. To maintain a balanced condition, rotor Joints must not permit relative 
shift to rotor components. Detail dimensions of rotor elements must be 
selected to maintain tightness at all operating conditions without 
yielding at any others. Elastic, thermal, and centrifugal aspects must 
be included in the analyses. 

2. Alternative HPFTP rotor configurations with fewer joint alignment sur- 
faces could potentially reduce relocation variations. Two configurations 
were presented. 

3. Improvement in the location of the HPOTP preburner impeller can be 
obtained by providing an outboard pilot surface; two concepts were 
presented. 

4. Bearing tilt does not significantly alter the dynamic response of the 
SSME high-pressure turbopump rotors. 

5. Bearing wear generally will not affect the balance of a rotor, but can 
affect the rotordynamic response by reducing bearing stiffness or intro- 
ducing nonsynchronous forcing functions. 

6. Turbine blade erosion of the magnitude and character experience on the 
SSME high-pressure turbopump does not introduce significant unbalance. 

RECOMMENDATIONS 

Recommendations for improving the balance procedures for the SSME high pressure 
turbopunips are: 

1. Near-term SSME balance procedure improvement: establish procedures for 

both HPFTP and HPOTP which include multiplane rotor-only balance using 
the influence coefficient method. The program should include: 

• Determination of minimum speed and number of runs required to 
product acceptable response for all speeds; initial balance 
speed = 6000 rpm 
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• Determination of the carry-through balance quality attained by 
conducting a comprehensive study of balance condition and response 
after repeated teardowns and reassemblies 

• Dead-band elimination 

• Use actual SSMF. rotor hardware modified so that reversible balance 
corrections suitable for high-speed testing can be made 

2. Incorporate trim balancing provisions into the HPOTP preburner impeller 

3. Develop hydrostatic dead-band eliminators for use in trim balancing and 
possibly in operation; investigate the effect of added damping at the 
bearing locations 

4. Incorporate rotordynamic sensing instrumentation into the high-pressure 
turbopumps, including: 

• Proximity sensors at: 

HPFTP liftoff seal sleeve 
HPOTP purged turbine seal 

• Force sensors - select either strain gages or high resolution 
proximitors at: 

HPFTP pump inlet end bearings 
HPFTP turbine end bearings 
HPOTP preburner end bearings 
HPOTP turbine end bearings 

5. Incorporate an outboard pilot into the HPOTP preburner impeller mounting 

6. Long-Term SSME Balance Improvement: Modify both high-pressure turbopumps 

to incorporate in-housing balance capability; modify the HPFTP casing 

to permit access to the third-stage pump impeller for correction; modify 
the HPOTP to permit access to the main pump impeller; incorporate 
exchangeable balance weights into both these impellers; concurrently, 
incorporate balance-monitoring instrumentation into the turbopump 
assemblies. Modify both high-pressure turbopump rotors to reduce the 
number of joints and minimize the use of nonindexable parts. 
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Incorporate in-housing balancing features into future turbopump designs 
including: 

• Balance Plane access 

• Exchangeable balance weights 

• Dead-band eliminators 
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TASK I: ELASTIC R0T0RDYNAM1C REQUIREMENTS 


The study of elastic rotordynamic requirements for the SSME high-pressure turbo- 
pumps entailed the following: 

1. Modeling the SSME HPFTP and HPOTP assemblies for dynamic response 
studies 

2. Defining the balancing requiresnents for the high-pressure turbopumps 

3. Identifying measurement techniques for dynamic response measurement. 

DYNAMIC MODELS 

The turbopump rctor casing elastic systems were modeled for finite-element analy- 
sis as vibrating beams with mass and elastic properties simulating bending stiff- 
nesses. These beams were connected to elastic casing models through compliant ele- 
ments simulating bearings. Unbalance was simulated by applying synchronous forcing 
functions to the rotor model. Fluid forces due to seal and turbine action are 
included as applicable. 

The dvnamic model stations or nodes are listed in Table 1-1 and shown in Fig. 1-1 
(HPFTP) and 1-2 (HPOTP). For the SSME HPFTP, the rotor model contains 13 nodes 
representing the positions of the major masses, bearings, and seals. The PPOTP 
model is similarly constructed using 14 nodes. In jach model, at each node, two 
radial orthogonal degrees of freedom are used to describe rotor response at that 
station. Adequate model response was ensured by including frequencies to three 
times maximum operating speed. The dynamic models used for this study include 
confirmation obtained by rap testing actual rotors and casings to verify natural 
frequencies and mode shapes. Gyroscopic effects are included for each major mass 
group. 

To simulate the effects of operation, the models used for the "wet" response 
studies were altered to reflect the effects ol operating temperature on the elastic 
materials properties, the additional mass of toe working fluid in rotor passages, 
ar.d the damping introduced by the presence of fluid. Bearing dead-hand effects 
were assessed using a simulation model to account for the nonlinearities introduced 
by these clearances. 
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Figure 1-2. HPOTP Rotor Model Node Locations 
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TABLE 1-1. ROTOR MODEL NODE LOCATIONS 


FUEL PUMP AXIAL LOCATIONS 


NUMBER LOCATION 


BEARING 1 
BEARING 2 
IMPELLER I 
INTERSTAGE I 
IMPELLER 2 
INTERSTAGE 2 
IMPELLER 3 
RING SEAL 

SECOND-STAGE TURBINE 
TURBINE INTERSTAGE 
FIRST-STAGE TURBINE 
BEARING 3 
BEARING 4 


LOX PUMP AXIAL LOCATIONS 


NUMBER LOCATION 


PREBURNER IMPELLER 
PREBURNER LABYRINTH 
BEARING 1 
BEARING 2 

LEFT MAIN IMPELLER INDUCER 
LEFT MAIN IMPELLER TIP 
RIGHT MAIN IMPELLER TIP 
RIGHT MAIN IMPELLER INDUCER 
BEARING 3 
BEARING 4 
RING SEAL 

SECOND-STAGE TURBINE DISK 
TURBINE INTERSTAGE 
FIRST-STAGE TURBINE DISK 


A description of th<» linear and nonlinear rotor casing system models and the equa 
tions used in their analysis is given in Appendix A. Influence coefficient 
matrixes were generated using the dynamic models for ufcfe in evaluating multiplane 
and in-housing balancing, and in evaluating the effects of correction locations. 

BALANCING REQUIREMENTS 

The definition of those elements that are req'> J red to achieve satisfactory dynamic 
response for elastic rotors was approached in tne iollowing steps: 

1. Permissible response limits were established for the SSME high-pressure 
turbopump. Allowable bearing loads and rotor deflections were defined 
using bearing fatigue life and seal clearances as criteria. 
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2. Maximum anticipated residual unbalance magnitudes were established using 
manufacturing tolerances and component masses. 

3. Linear elastic rotordynamic models of the HPFTP and HPOTP rotors and 
housings were used to predict responses to the residual unbalance magni- 
tudes arising from manufacturing tolerances and procedures. 

4. The effects of bearing dead bands were illustrated using a nonlinear sim- 
ulation model. Meats of linearizing the rotordynamic system by eliminat- 
ing dead bands were described. 

5. rhe capabilities and limitations of rigid rotor, multiplane, and flexible 
rotor balancing procedures were reviewed and compared for adapatation to 
rotor-only and in-housing balancing of the SSME high-pressure turbopumps. 

Permissible Rotor Response 

Limitations for bearing loads and rotor deflections were defined and are shown in 
Fig. 1-3 for the HPFTP and Fig. 1-4 for the HPOTP. The criteria for these limita- 
tions were as follows. 

Bearing Loads . Bearing radial loads due to unbalance were limited so as not to 
exceed the axial preload applied to the bearing. Design values for the preloads 
in the SSME high-pressure turbopumps are: 

• HPFTP (all bearings) - 700 pounds 

• HPOTP (preburner end) - 800 pounds 

• HPOTP (turbine end) - 1100 pounds 

The preceding limits apply at all speeds within the operating range. Radial loads 
at steady-state speeds were limited to those corresponding to a fatigue life of 
7.5 hours, assuming normal axial preload is maintained. Bearing life relations 
are shown in Appendix B. 

Rotor Deflections . Since balancing operations may event-ally be performed at high 
speed at room temperature, radial clearances between parts with potential for rub- 
bing were listed for both ambient and chilled operating conditions. Ambient 


RI/*D82-285 

1-4 
















clearances would be applicable to dry dynamic responses as defined in the rotor- 
dynamics studies, while the chilled clearances would be associated with the wet 
responses. 

Res idua l Un balance Magn itudes 

Maximum expected unbalance values were established by adding the residual unbalance 
tolerances to those resulting from assembly runout limits for the HPFTP. It was 
assumed that these runout limits represent worst-case relocation variations. The 
resulting unbalances are shown in Fig. 1-5. The total runout limits for the HPOTP 
based on drawing tolerances were extremely low; to establish realistic total 
unbalances, a TIR of 0.001 was assumed for all components. When added to permis- 
sible residual unbalance (prior to in-housing balance), the maximum total unbalance 
predicted agrees with unbalance values actually found in the in-housing balance 
operation for the turbine end. These are shown in Fig. 1-6. 

Unha 1 anc e Re spouse 

The linear response of the high-pressure turbopumps to unbalance was calculated 
over the operating speed range for both "wet" operation, t.e., with working fluids 
present, and for "dry" running, as would occur in a balancing procedure. 

Modes. The predicted critical speeds of the high-pressure turbopumps are indicated 
in Fig. 1-7 (HPFTP) and 1-8 (KPTOP). The rotor mode shapes arising at these 
resonant speeds are shown in Fig. 1-9 and l- 10. Observation of the mode shapes 
indicates that best observation of the first mode would be mid-span for the HPFTP 
and at the turbine end on the HPOTP. Second mode detection would be enhanced at 
the bearing locations. 

Controllability . Optimum control of unbalance response requires balance correction 
or weight adjustment at the most active planes. For the HPFTP, corrections at the 
second- or third-stage pump impellers is best for control of first critical speed 
response, while the turbine and inlet end would affect second mode response more 
strongly. In the HPOTP, present operation is below, but close to, the second 
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Figure 1-5. HP FTP Unbalance Values 




Figure 1-6. HPOTP Unbalance Values 
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Figure 1-7. HPFTP Critical Speeds 
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critical speed. Corrections at the turbine disks provide effective control of the 
first mode; however, corrections at the main impeller are desirable for attenuating 
response at the second mode. Accessibility of any but the outboard components in 
the existing turbopump assemblies requires penetration of high-pressure housings 
for access to perform corrections. Secure attachment of balance weights, mainten- 
ance of propellant clean conditions, and security of closure ports to prevent high- 
pressure propellant leaks are absolute requirements for balancing procedures 
involving corrections at inboard locations such as second-stage turbine disks or 
main pump impellers. 

Modal Truncation . The SSME response models include modes resulting from three 
times the maximum spin speed. This practice ensures that significant response 
content from higher modes is present and, therefore, that adequate accuracy of 
the analytical results is obtained. 

Bearing Loads Due to Unbalance . The effect of unbalance at each of the balancing 
planes is contained individually in the influence coefficients. For example, 

Fig. 1-11 shows the load on the No. 3 and 4 bearing pair as a function of speed 
for a unit (1 gram-inch) unbalance at plane III, the second-stage turbine wheel. 

When the effects of individual unbalances at each plane are summed, the total 
bearing reactions in response to any given unbalance distribution can be obtained. 

To define the maximum anticipated bearing loads, worst-case unbalance distributions 
were used. 

Worst-Case Unbalance Distribution . The combination of unbalance vector phase 
angles resulting in the highest peak bearing loads was designated the worst-case 
unbalance distribution, and was identified in the following manner: 

1. The maximum anticipated residual unbalance magnitudes estimated previously 
(shown in Fig. 1-5 and 1-6) were arranged in various combinations of 
phase angles of 0 and 180 degrees, listed in Table 1-2 (HPFTP) and 

Table 1-3 (HPOTP). 

2. The resulting unbalance distributions were input to both the dry and wet 
linear .ponse programs. Speeds to 45,000 rpm were included in the 
respo..^ study to illustrate high-speed effects; actual operation at 
this speed would be precluded by strength limitations. 
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TABLE 1-2. HPFTP UNBALANCE ORIENTATIONS, UNBALANCE RESPONSE STUDY 






BALANCING PLANE 

CASE 

1 

2 

3 

4 

5 

1 

♦ 

4 

4 

4 

4 

2 

- 

♦ 

4 

4 

4 

3 

♦ 

- 

4 

4 

4 

4 

♦ 

♦ 

- 

4 

4 

5 

♦ 

4 

4 

- 

4 

6 

♦ 

4 

4 

4 

- 

7 

- 

- 

4 

4 

4 

8 

♦ 

- 

- 

4 

4 

9 

4 

4 

- 

- 

4 

10 

♦ 

4 

4 

- 

- 

11 

- 

4 

- 

4 

4 

12 

♦ 

- 

4 

- 

4 

13 

4 

4 

- 

4 

- 

14 

- 

4 

4 

- 

4 

15 

♦ 

- 

4 

4 

- 

16 

- 

4 

4 

4 

- 


— 4 

7 

7 

5 

5 — 


MAXIMUM ANTICIPATED UNBALANCE 
MAGNITUDE. G-IN. 


h CORRESPONDS TO A PHASE ANGLE OF 0 DEGREES 
- CORRESPONDS TO A PHASE ANGLE OF 180 DEGREES 







TABLE 1-3. HPOTP UNBALANCE ORIENTATIONS, UNBALANCE RESPONSE STUDY 



BALANCING PLANE 

CASE 

1 

2 

3 

4 

1 

4 

4 

4 

4 

2 

- 

4 

4 

4 

3 

♦ 

- 

4 

4 

4 

♦ 

4 

- 

4 

5 

4 

4 

4 

- 

6 

- 

- 

4 

4 

7 

4 

- 

- 

4 

8 

4 

- 

4 

- 


— 2 

5 

5 

5 — 


MAXIMUM ANTICIPATED UNBALANCE 
MAGNITUDE. G-IN. 


♦ CORRESPONDS TO A PHASE ANGLE 
OF 0 DEGREES 

- CORRESPONDS TO A PHASE ANGLE 
OF 180 DEGREES 
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3. The combination producing the most severe response was identified as the 
worst-case unbalance phase angle distribution. 

This combination was subsequently used to evaluate balancing procedure improvements. 

Peak loads were predicted for the unbalance distributions listed in four tables: 


Table 

Turbopump 

Condition 

Response 

Figure 

1-4 

HPFTP 

wet 

1-12 

1-5 

HPFTP 

dry 

1-13 

1-6 

HPOTP 

wet 

1-14 

1-7 

HPOTP 

dry 

1-15 


Considering the full operating speed range. Fig. 1-12 shows that the bearing loads 
resulting from the worst-case unbalance distribution exceeded the load criteria 
previously established (Fig. 1-3); the transient radial load at 22,000 rpm exceeds 
the 700-pound iimitation previously selected. In addition, the steady-state 
radial loads at RPL and FPL exceed the load/life criteria. Figure 1-13 indicates 
that extremely high transient shaft deflections can result if the HPFTP with 
worst-case unbalance is spun dry at the first bending resonant speed of approxi- 
mately 14,000 rpm. This result could be intuitively inferred by considering that 
in Case 1, all the unbalances are in the same direction and would reinforce the 
mode shape at 14,000 rpm (Fig. 1-9). From a practical standpoint, dry operation 
would be entailed only in a balancing operation and the excessive loads would be 
avoided by making progressive balance corrections before increasing speed to the 
resonant condition. 


Balance Improvement 


Although nonlinearity introduced by the bearing clearances or "dead bands" can 
alter the turbopump's actual response from that indicated by the linear model, 
insights concerning the effectiveness of balance improvement can be gained from 
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TABLE 1-6. WORST -CASE UNBALANCE STUDY 


WET LOX PUMP RESPONSE 


MAXIMUM 

UNBALANCE RESPONSE 

BEARING NO. 

BEARING LOAD, 
POUNDS 

SPEED, 

RPM 

1 

2938 

35,000 

1 

1732 

35,000 

1 

1775 

45,000 

1 

2696 

35,000 

1 

3116 

35,000 

1 

2878 

35,000 

1 

1925 

35,000 

1 

1500 

35,000 


* WORST CASE 


TABLE 1-7. WORST-CASE UNBALANCE STUDY 


DRY LOX PUMP RESPONSE 


MAXIMUM 

UNBALANCE RESPONSE 

BEARING NO. 

BEARING LOAD, 
POUNDS 

SPEED, 

RPM 

1 

4180 

37,000 

1 

5900 

45,000 

1 

8249 

45,000 

1 

3791 

35,000 

1 

5083 

35,000 

1 

4700 

35,000 

1 

8040 

45,000 

1 

6108 

45,000 


*W0RST CASE 
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ROTOR SPEED. RPM 


Figure 1-12. HPFTP Wet Response, Worst-Case Unbalance 
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HPf TP -uncorrected rotor response 



Figure 1-13. HPFTP Dry Response, Worst -Case Unbalance 



estimates based on linear response analysis. Using this approach, values for the 
permissible unbalance, potential benefit of in-housing balance, and the effects 
of data accuracy and instrumentation-associated noise were made. 

Permissible Unbalance . It is evident that the HPFTP dry rotor has a more highly 
peaked response to worst-case unbalance distribution than does the turbopump con- 
taining fluid. For example, using the wet response, if a linear relationship of 
load to unbalance is assumed, the worst-case unbalance magnitudes must be reduced 
to 73% of the maximum anticipated unbalance, magnitudes to meet operating bearing 
load criteria (Table 1-8). By contrast, if the dry turbopump is spun to 14,000 
rpm, rubbing will occur at the turbine seal, station 10, unless the unbalance is 
reduced to 10% of the worst-case magnitude. This ratio falls to 8% if overloading 
of the No. 3 bearing is to be avoided during a dry spin. This value is probably 
conservative because the bearing outer diameter radial clearance at the dead band 
will absorb some of the bearing deflection, reducing the loading. 

Applying similar criteria to the HPOTP, it was found that bearing load and shaft 
deflection limits were not exceeded at speeds below 32,000 rpm for either the 
wet or dry responses (Fig. 1-14). Table 1-9 lists the radial deflections for seals 
at stations 1, 2, 11, and 13. The bearing stiffnesses were used to estimate an 
allowable radial deflection for comparison to the predicted deflection from the 
dry response profile (Fig. 1-15). 

Benefit of In-Housing Tr>m Balance . A beneficial reduction in bearing loads can 
be achieved if a trim balance can be performed, correcting the individual masses 
of the rotor while not disturbing the assembled state of the turbopump. For pur- 
poses of illustration, it was assumed that, in an optimized balance process, an 
assembly with an initial as-assembled worst-case balance condition was corrected 
to 1 gram-inch unbalance at each plane. 

For the HPFTP, the wet response was reduced from that shown in Fig. 1-12 to the 
lesser values of Fig. 1-16. The peak transient load at bearing No. 4 at 22,000 rpm 
was reduced to 178 pounds from 920, while the maximum steady-state load at 32,000 
rpm was lowered from 620 to 122 pounds, well within the residual balance limit 
criteria shown in Fig. 1-3. 
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TABLE 1-8. ESTIMATE OF PERMISSIBLE HP FTP UNBALANCE 


A. WET HPFTP - BEARING LOAD CRITERIA 

LOAD LIMIT, POUNDS PEAK LOAD, POUNDS 

TRANSIENT STEADY STATE TRANSIENT 
STATION (ALL RPM) (32,000 RPM) (22,000 RPM) 



PERMISSIBLE + WORST CASE 


TRANSIENT I STEADY STATE 


B. DRY HPFTP - SEAL CLEARANCE AND BEARING LOAD CRITERIA 



ROTOR DEFLECTION LIMIT, 
INCH 

(ALL RPM) 

MAXIMUM DEFLECTION, 
INCH 

(14,000 RPM) 

0.0115 

0.022 

0.0115 

0,033 

0.006 

0,033 

0.003 

0.029 

0.00175 

0.022 


PERMISSIBLE i WORST CASE 
0.52 
0.35 
0.18 
0.103 
0.Q8 


TABLE 1-9. ESTIMATE OF PERMISSIBLE HPOTP UNBALANCE 
DRY HPOTP - SEAL CLEARANCE AND BEARING LOAD CRITERIA 



ROTOR DEFLECTION LIMIT, 

MAXIMUM DEFLECTION, 



INCH 

INCH 


STATION 

(ALL RPM) 

(14,000 RPM) 

PERMISSIBLE * WORST CASE 

1 

0.0070 

0.0012 

>1 

2 

0.0075 

0.0012 

>1 

3 ) 

0.0016 

0.0013 

>1 

4 1 




9 ) 

0.0022 

0.0016 

>1 

10 1 




11 

0.0146* 

0.0065 

>1 

13 

0.012 

0.0082 

>1 


INDICATES TOTAL OF SHAFT AND HOUSING CLEARANCES 


RI/RD82-285 


















RI/RD82-285 

1-25 


3.XI0 


H.X10 




i 


rst-Case Unbalance, 









ORIGINAL FAG- i- 
OF POOR QUALITY 



ance 



ORIGINAL PAG'c [3 
OF POOR QUALITY 



ROTOR SPEED, RPM 


Figure 1-16. HPFTP Wet Response Optimized Balance, 
Bearing No. 4 
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The wet HPOTP optimized balance response (Fig. 1-17) shows a reduction in load at 
the preburner end bearings at 34,000 rpm from the limit (Fig. 1-4) of 1600 pounds 
to 440 pounds, also well within the limits established for bearing load (Fig. 1-4). 

Accuracy. Proximity probe manufacturers cite resolution of better than 50 micro- 
inches. Applied to the component weights shown in Fig. 1-5 and 1-6, this resolu- 
tion should permit detection of less than 1 gram-inch for any plane in either SSME 
high-pressure turbopump. Rationale for this estimate is shown in Appendix C. 

Noise . Extraneous signal components may cause a background level of noise which 
will obscure the desired response signal. The noise components arise from thermal 
source«' in the system, power supply ripple, electrical or mechanical runout, and 
stray electrical signals conducted or induced into the signal conductors. Noise 
interference can be reduced by utilizing a band-pass filter to exclude those fre- 
quencies outside the range of interest. For a turbopump with a synchnonous spin 
frequency of 500 Hz, an 800 Hz low-pass filter is typically used between the 
signal source and the recording device. A tracking filter is often applied to the 
taped data to emphasize the synchronous content. Synchronous noise such as mech- 
anical and electrical runout must be treated by analytic subtraction, as discussed 
in Task II. 

Dead-Band Effects 


"Dead band" is the terra often used to describe the radial clearance between the 
ball bearing outer race and the mating housing bore (Fig. 1-18). The effects of 
these clearances were assessed by comparing responses to unbalance as determined 
using a linear rotordynamic model with those responses indicated by a nonlinear 
rotordynamic model. The results (summarized in Table 1-10) indicate that: (1) 

the superposition principle, which is central to linear analysis, is not valid for 
systems with radial clearance; (2) peak bearing load may be reduced by the intro- 
duction of small radial clearances; and (3) the peak bearing load may be increased 
as the radial clearance is enlarged. The term may be is used to modify these con- 
clusions because the solution found may not be unique; multiple solutions result 
from the nonlinear characteristic introduced by the clearance. 
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Figure 1-17. HPOTP Wet Response, Optimized Balance, 
Bearings No. 1 and 2 
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TABLE 1-10. DEAD BAND EFFECTS, HPFTP ROTOR RESPONSE 


A. LINEAR RESPONSE (NO DEAD BAND) 


CASE 

UNBALANCE, 
G-IN . 

A 

1 

B 

2 

C 

3 



NO. 4 BEARING 
PEAK LOAD, 
POUNDS 



B. NONLINEAR RESPONSE (DEAD BANDS INCLUDED) 


CASE 

UNBALANCE, 

G-IN. 

CLEARANCE, 

INCH 

NO. 4 BEARING 
RMS LOAD, 
POUNDS 

NO. 4 BEARING 
PEAK LOAD, 
POUNDS 

D 

1 

0.0 

30.0 

42.0 

E 

1 

0.00025 

23.2 

33.0 

F 

1 

0.00075 

62.7 

89.0 

G 

1 

0.00150 

106.0 

150.0 

H 

2 

0.00025 

33.1 

47.0 

I 

2 

0.00075 

64.6 

91.3 

H 

2 

0.00150 

125.0 

177.0 

II 

3 

C 00025 



n 

3 

0.00075 

69.6 

98.0 

M 

3 

0.00150 

126.3 

179.0 
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Figure 1-18. Bearing Clearance or Dead Band 


Although it Is not clear that elimination of dead bands will improve turbopump 
dynamic performance, their presence is considered deleterious to accurate detec- 
tion of unbalance because balance methods are based on linear assumptions. Means 
of eliminating dead-band effects during balancing operations include use of a 
spring-loaded plunger (Fig. 1-19) pressurized plunger (Fig. 1-20), or hydrostatic 
film (Fig. 1-21). Of these methods, only the hydrostatic film preserves the axial 
freedom of motion of the outer races. Therefore, it is the only aieans adaptable 
to turbopump operation. An additional potential benefit of the hydrostatic film 
is addition of damping at the bearings. 


Balancing Procedures 

In evaluating means to achieve improved ualance for the SSME high-pressure turbo- 
pumps, the various aspects of alternate methods were considered before arriving 
at the recommended procedur. 


RI/RD82-285 

1-31 





ORIGINAL PAGE IS 
OF POOR QUALITY 





For the near-term improvement of SSME high-pressure turbopump balancing, it was 
concluded that the optimum method is a multiplane, moderate speed (6000 rpm) , 
influence coefficient-based balance performed on the rotor prior to turbo- 
pump assembly. Hydrostatic or other means of dead-hand elimination should be 
employed. Assembled check and trim balance procedures are to be retained with 
modifications. Improved balanced conditions are expected because, when correc- 
tions are applied to the location where tin* unbalance previously existed, no 
internal moments remain. The present practice of applying two-plane balance 
corrections to assemblies that may be incomplete or which have fixtures sub- 
stituted for rotor components produces valid rigid rotor balancing. However, 
internal moments may remain which may augment rotor deflections, thus reinforcing 
rotor responses. 

For the HPFTP, the relocation variation causing a large proportion of the total 
anticipated residual unbalance arises (Fig. 1-5) from the substitution of a 
sleeve for the second- and third-stage Impellers In tin* first balance assembly 
(see Appendix D) . A reduction In tin* magnitude of this variation is expected 
when all planes are balanced in the full rotor assembly. The in-housing check 
balance is to be continued as a means of detecting anomalies arising during the 
assembly process, with the change that slave bearings would not be substituted 
for the operational bearing set. This would eliminate. the relocation variation 
caused by differences in bearing inner race eccentricity. 

The HPOTP trim balance has been effective in reducing the rotor response at 
speeds to FPL. This result is due to relocation variations of the first-stage 
turbine wheel and turbine blades which are the major cause of residual unbalance 
shift between balance assembly and final assembly. Correcting all major masses 
at the same time, prior to final assembly, should reduce internal moments, thus 
improving the results obtained bv the current superposition of separate two-plane 
balances. Preburner impeller corrections incorporated into the trim balance 
procedure would improve the balance condition attained. 

Elimination of dead bands during trim balance operations will enhance the validity 
of the balance procedure which is based on linear response. 
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The near-term improvement attainable is limited by practical considerations of 
access and correction in the assembled turbopump. Assuming solution of these 
constraints, alternate balance procedures with potential for higher quality balance 
results were considered as outlined below. 

Ho t -F iring Based Ba lance . Potentially, the most accurate balance procedure 
incluues corrections based on dynamic operational data. These corrections may be 
made through the turbuopump casings, eliminating disassembly after balance, thus 
avoiding relocation variations. Unbalance data collected during hot firing could 


contain 

all operational 

factors including: 

1 . 

Speed 


2. 

Power 


3. 

Temperature 


4. 

Synchronous fluid forces 

5. 

Destabll iz ing 

influences 


• Seal forces 

• Impeller clearance forces (Alford forces) 

(Nonsynchronous factors would normally be filtered out of 
the signal .) 

Ft m a practical standpoint, however, unbalance detection is more difficult during 
engine operation: 

1. The number of Instrumentation sensor locations is restricted. 

2. Invasive instrumentation increases leakage and ignition hazards. 

3. The dynamic response of the rotor/housing system is damped by the 
working fluids. 

4. Multiple tests for progressive correction and repeatability checks 
are prohibitively expensive. 
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In-housing corrections to inboard components will be difficult to implement. Assum- 
ing the availability of unbalance detection methods and balancing programs (soft- 
ware) capable of identifying the required corrections, the structures of the SSME 
high-pressure turbopumps preclude easy access to rotor elements for implementation 
of the correction. In both high-pressure pumps, correction access to main impellers 
entails penetrating highly stressed pressure containments with attendant sealing 
and security requirements. In addition, multiple walls must be penetrated to reach 
the HPFTP impellers. Further, assuming that a satisfactory method of access is 
developed, a means of implementing the correction while maintaining LOX-clean con- 
ditions is necessary. 

H igh Speed In -Housing Ba lance . A second method of improved balancing includes a 
high-speed, in-housing trim balance conducted at ambient temperature after assembly. 
The assembled turhopumps will be electrically driven in an evacuated environment. 
Advantages of this procedure are: 

1. Elastic response to speed is included. 

2. Response is undamped, increasing unbalance signal magnitudes. 

3. Sensor probe environments are moderate. 

4. Sensor probes can be removed prior to hot-firing operation. 

5. Balancing process control is enhanced 

• Speed can be controlled within small limits. 

• Multiple balance runs are feasible. 

Balancing performed in this way lacks complete simulation of operating conditions, 
among which are: 

1 . Temperature 

2. Power effects 

• Torque 

• Pressures 

3. Synchronous fluid effects 
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An additional requirement is that the bearings be cooled and lubricated to avoid 
damage during balance operation. 

In-housing correction poses the same difficulties noted for the hot-fire test 
method: access to the rotor through the high-pressure turbopump casings for 

correction would require extensive redesign of the turbopumps and would increase 
the hazard of leakage and contamination. 

Multiplane Balancing . Mulrimass rotors such as those of the HPFTP and HPOTP can be 
balanced so that internal moments are eliminated if corrections are made to each 
major mass for its individual unbalance. Standard balance machines cannot perform 
this operation; however, the required corrections can be calculated from influence 
coefficients which are input to a multiplane balancing computer program. As sub- 
sequently shown in Task II, low response can be obtained over the entire speed 
range with a multiplane balance carried out at moderate (below resonant) speeds. 

The accuracy of the corrections can be improved by utilizing data taken at higher 
speeds, which includes elastic properties of the rotor. Because it requires access 
to each major mass, this method is difficult to adapt to in-housing balancing, but 
is applicable to complete multiplane rotor assemblies mounted in fixtures outside 
the pump casings. The disadvantage of the system is that, whilt internal moments 
are minimized, potential relocation variations arise when the rotor is disassembled 
for reassembly into the turbopump casing. 

MEASUREMENT TECHNIQUES 

Balancing techniques based on influence coefficients can utilize phase-correlated 
signals from essentially any source with linear response to shaft motions or forces. 
Accelerometers, def lectomers , velocity or proximity pickups, force measuring devices, 
ox* possibly even pressure transducers are potentially useful. The signal process- 
ing block diagram shown in Fig. 1-22 is adaptable to outputs from any of these 
potential sources, and can be applied to determination of the following conditions 
by inference. 
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Condition 


Evaluation Criteria 


Unbalance 

Resonance 

Damping 


Modes 

Whirl type 
Dead bands 


Structural changes 
Bearing defects 


Stack loosening 


Rubbing 


Synchronous content of composite signal 

Peaked response at resonant frequency 

Shifting or suppression of resonant peaks. 
Measurement requires knowledge of forces and 
relative velocities at interface where damping is to 
be quantified. Fourier analysis of signals contain- 
ing these elements can yield stiffness and damping 
coefficients. 

Comparison of relative response of several sensors 
as a function of their axial position. 

Comparison of whirl frequency and peak displacement 
Phase to that of a synchronous reference. 

Presence of dead band may be inferred from nonlinear 
response. Analysis of the signal from a sensitive 
radial deflection sensor could be used to assess the 
presence and magnitude of dead bands. For example, 
the signal from a fiber optic sensor (see Fig. 1-231 
observing the outer race deflections would contain: 

1. Local deflections of the race caused by indi- 
vidual ball contacts due to axial, radial, and 
centrifugal loading. These deflections will be 
in the microinch range. 

2. Motions of the outer race across its dead band 
will superimpose relatively much larger modula- 
tion on the individual ball load deflections at 
the frequency and phase of the predominant 
radial load. Thus, a fixed radial load will 
produce a d-c shift of signal, while a rotating 
radial load will add a large a-c component at 
the load's frequency. The magnitude of these 
large shifts could be related to the dimensions 
of the dead band. 

Changes in response data independent of spin speed 

Changes in response independent of spin speed. 
Analysis of signals from accelerometers, acoustic 
emission devices, or def lectometers can be used to 
identify the cause. 

Changes in response independent of shaft speed. 
Harmonics of shaft speed appear. See Fig. 1-24 for 
accelerometer response to a loosened HPFTP tie bolt. 

Increase in amplitudes of harmonics of shaft speed 
over normal amplitudes (Fig. 1-33). 
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Motion and Force Sensor Concepts 


Data most adaptable to balancing and rotordynamlc diagnostic information consist 
of rotor motion and bearing loading. Consideration of the rotordynamics response 
of the SSME high-pressure turbopumps led to the selection of the sensor classifi- 
cations and locations listed in Table 1-11 and shown in Fig. 1-25 and 1-26. The 
mode shapes of the rotors at resonances suggest that rotor deflection at stations 
away from the bearings is most useful for detection of response at the rotor first 
critical speed, while the bearing reactions are more responsive near the second 
critical speed. Therefore, for the operating speed ranges of the high-pressure 
turbopumps, bearing loads appear to be the most desirable measurements to describe 
rotor response. Several measurement techniques may be adaptable to bearing load 
determination including proximity (shaft radial deflection is related to bearing 
load by bearing stiffness) strain gages, piezoelectric cells, or photoelectric 
sensors. 

Pulsed X-rays have been considered as a noninvasive measurement technique to deter- 
mine rotor motions. Testing using a J-2 turbopump indicated that X-ray energy 
levels of 2 MV is required to produce images of internal components. However, 
resolution was concluded to be unsatisfactory for measurement of unbalance, as it 
was estimated to be no better than ±0.125 mm (±0.005 in.) in vibration amplitude 
(Ref. 2). 

All sensor signals must be phase-related to produce valid rotordynamlc data suit- 
able for balancing. This is normally accomplished by providing an index on the 
rotor; proximity sensors can be phase related by providing a depression of limited 
depth (~0.005 inch) and angular extent (~45 degrees). The depth of depression also 
provides a calibration for the sensor's scale factor. Photoelectric probes can be 
phase referenced by providing a contrasting surface reflectivity. 
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TABLE 1-11. POTENTIAL SENSOR LOCATIONS 
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Figure 1-25. SSME HPFTP Instrumentation for Improved Rotor Balancing 




for Improved Rotor Balancing 



Probe Types for Me asur ement of R adial Shaft Displacement . Proximity sensing tech- 
nique^ suitable for a high-speed rotating shaft can utilize the following sensor 
types : 

• Capacitive 

e Eddy current 

• Magnetic reluctance 

• Optical 

In principle, all of these techniques can work with the sensor and target immersed 
in liquid 0^ or H,^. Any of these sensor techniques can be made to expose a prac- 
tically acceptable (i.e., minimal) amount of organic materials to the propellants. 
For balancing purposes, the static (zero) stability required of the sensors is of 
secondary concern. In Table 1- 1 2, a qualitative summary and rating is given of 
sensors using the listed techniques. The numbers given for frequency response and 
resolution are approximate and a function of the probe size. 

Appendix E provides a description of the four probe types mentioned above with 
additional background information on inductive sensors. 

Practical considerations rule out positioning of probes in labyrinth seal regions 
due to the surface configuration. In addition, measurements in the hot high- 
pressure region of the turbine could only be considered for special investigations. 
For example, if the bending mode of the LOX pump rotor at 11,028 rpm could only be 
detected with sufficient sensitivity through a measurement of the turbine disk 
excursions, and if that measurement must be made on the operating pump, then high- 
temperature/high-pressure resistant sensors must be used. In the following listing, 
however, the turbine region is not represented in the pressure snd temperature 
specif icat ions . 

The target specifications for the displacement sensor are: 

• Frequency response 

• One per revolution (700 Hz) for balancing 

* Multiple per revolutions (5 to 20 Hz) for diagnostics 
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TABLE 1-12. FEATURES OF PROXIMITY SENSORS 
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• Displacement range: 0.020 

• Linearity: ±2X 

• Resolution: 0.0001 inch 

• Static stability: 10% of reading 

• Standoff (mean gap): 0.020 inch minimum 

• Target area available: 1/8 to 1/4 inch; 1/2 inch in a few locations 

• Cable size: 1/8 inch maximum 

• Probe size: 3/8-inch-diameter (max)-l/2 inch in a few locations; 

1 inch long (max) 

• Pressure resistance: 300 psi maximum 

• Temperature range: 15G R to room temperature in LO 2 

40 R to room temperature in 

• Vibration (typical): 10 Hz at 20 g to 2000 Hz at 400 g 

The rating order in Table 1-12 is based on the following criteria in order of 
importance . 

• Effect of medium (liquid versus gas) in the gap on output 

• Noise caused by magnetic fields 

• Size (diameter and length), and target size needed 

• Simplicity, ruggedness 

• Frequency response 

• Dynamic resolution 

• Linearity 

• Sealability 

• Static stability 

The conventional eddy current probe is rated most applicable to cryogenic turbo- 
pumps because its output is not affected by the medium in the gap or by magnetic 
fields. In general, the eddy current probe is a pancake coil that requires a 
dielectric support to avoid eddy currents in back of the coil. This requirement 
determines the minimum length of this type of probe. The dielectric space in 
back of the coil must be equal to 1 to 2 diameters of the coil. Around the probe 
tip, a dielectric space (free from metal) of about 2 coil diameters is required. 
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The fiber optic probe is also very attractive for this application because of its 

insensitivity to magnetic field noise, ability to acquire data from small targets, 

simplicity of design, and ruggedness. In additon, the fiber optic probe requires 

a minimum of space, except in the axial direction (of the fibers). The presence 

of potting resin in the probe face (over about 15% of the face area) is probably 

-4 

acceptable (in the LOX pump) as the absolute area is very small (4 x 10 inch). 

The fiber optic probe provided with an ("extender") lens combines a long standoff 
distance with high sensitivity. However, this probe is affected by changes in 
the gap medium such as cavitation and therefore is rated second below the eddy 
current probe. If it is certain that cavitation will not occur, fiber optics 
would be the probe type chosen. 

The capacitive probes are rated third and fourth because of greater size, greater 
complexity, and high insulation requirements. The highly linear guarded probe 
system has poorer dynamic response and, therefore is rated below the unguarded 
(less linear) system. 

The magnetic reluctance probe is rated least desirable because of the deleterious 
effects of magnetic fields on the output. These probes have a magnetic core, 
which channels the flux very effectively, but the nonlinear magnetization prop- 
erties of the core material result in spurious signal output when stray magnetic 
fields are present. If the target is magnetic, the effect of stray fields is 
increased. The most attractive feature of this probe is that it can be made 
hermetic. It tolerates a substantial thickness (0.010 to 0.040 inch) of non- 
magnetic material in the gap whereas the eddy current probe will allow only 0.001 
to 0.002 inch of covering material in the gap. 

Locations for Displacement Measurements . The most desirabl measurement locations 
are those with maximum shaft amplitude for the mode under study. A compromise must 
be made between the optimum locations for various mode shapes and practical 
constraints. 


RI/RD82-285 

1-50 


Candidate proximity probe locations are shown in Fig. 1-26 at location F3 and In 
Fig. 1-25 at locations 03 and 04. These locations are also marked in the mode 
shape diagrams (Fig. 1-9) and In response charts (Fig. 1-13 and 1-15). For mea- 
surements at the first critical speed, locations F3 and 04 show amplitudes of 
more than 50Z of the maximum . For the second mode, locations F3 shows an ampli- 
tude of 25% of the maximum; while 03 coincides with over 50Z of maximum; thus, 
at all of these locations, a usable balancing signal can be expected. At 04, the 
operating temperature range of 480 to 640 F requires a medium temperature capabil- 
ity sensor. The commercially available proximity probes (of the types listed in 
Table 1-12) cannot meet this temperature requirement, but special probes can be 
obtained (Ref. 2). The radial and axial space available at location 04 between 
the two seal rings can accommodate a button-type eddy current proximity sensor. 

The gaseous helium pressure at 04 is less than 250 psia, which requires no special 
provisions in the probe design. In addition, signal leads can be routed through 
drain passages. 

At location F3, (Fig. 1-25), the maximum pressure Is 5400 psi, which may pose 
stress and leakage problems for any of the proximity probes and the connecting 
leads or fiber bundle. Hard-line cable sheets can be developed to resist high 
pressures by routing them through a compression fitting, as shown in Fig. 1-27. 
Small-diameter cable (0.040 to 0.060 inch) must be used to obtain required 
f lexibil ity . 

For fiber optic sensors, the fiber bundles inside cables must be left unattached 
to the sheathing if the cable is to be sufficiently flexible for Installation. 

Only at the ends of the cable can the fibers be potted to prevent leakage. To 
withstand the high pressure, special self-sealing tapered optical "lead" termina- 
tions must be developed, i.e., the one sketched in Fig. 1-28. Class fibers also 
can be fused thermally, but this process lias not been developed. 
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Figure 1-28. Pressure-Resistant Fiber Optic Probe (Concept) 

Axial motion of the rotor can be monitored by proximity probes installed in the 
pump casing, as shown in Fig. 1-29. Routing and sealing of the signal conductor 
can be accomplished in a manner similar to that for the radial motion sensor. 

Bearing Load Measurement . The most direct means of determining bearing loads is 
application of strain gages to the bearing mountings. The disadvantages of strain 
gages lie in their fragile nature and the complication of the assembly procedure 
resulting from the need to route conductors to the casing exterior through avail- 
able passages. Compatibility with an LC^ environment can be obtained by coating 
the gages with a compatible protective material. Other means of inferring bearing 
load from rotor motions and bearing race deflections can utilize sensitive proxi- 
mity measurements. However, the motions permitted by dead bands must be taken 
into account in interpreting the data. 
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Figure 1-29. Axial Motion Sensor Concept 


RI/RD82-285 


An alternate means of semi-direct load measurement requires incorporation of 
piezoelectric cells into the bearing housing. This method requires sensor loca- 
tion at the site of load application and, therefore, is less effective for moni- 
toring rotating radial load. Since the magnitudes of dead hands can he large 
compared to elastic deflections caused hy load, the data from piezoelectric 
sensors is subject to large nonlinear indications. 

Fiber Optic Sensors . Direct radial load isolation is not possible because 
the main axial load is carried by balance pistons. Transient axial loads on the 
HPFTP rotor can be isolated by measuring deflection of the thrust bearing support. 
The latter concept is illustrated in Fig. 1-30 where the axial displacement of 
the transient thrust bearing is measured with a fiber optic probe and a lens 
The use of the lens allows a large standoff of the probe, which provides space 
for a window with pressure seal. The displacement resolution and frequency 
response of the fiber optic probe is considered sufficient for this application. 

Figure 1-31 shows a radial load measurement concept that uses a dual fiber optic 
probe. The probe senses the radial motion in two mutually perpendicular direc- 
tions. The principle of operation of the probe is illustrated in Fig. l-3lb 
which shows three fibers in section and a narrow target (0.005- inch thick). The 
center fiber illuminates the target, and the outer fibers receive reflected light. 
The difference in the intensity (Ij-I.,) received by the two outer fibers is a 
function of latei\il displacement (X) of the target, as sketched in the respon^ 
in Fig. l-31b. The gap in this measurement is assumed to be constant and is 
adjusted to be at the optical peak of the response of the sum (1^ + I^) of the 
two outer fibers. Thus, small changes in the gap will have only second-order 
effects or. the output of the probe. The ratio of the difference and sum: 

(Ij - I )/(Ij - Ij + 1^) will be a signal that is compensated for gap effects. 

If necessary, this type of gap compensations can be incorporated in the probe 
readout. To increase the signal-to-noise ratio, and to eliminate microscopic 
effects of target reflectivity, the three fibers shown in Fig. 1-31 are expanded 
into lines of fibers containing typically 50 to 100 fibers each. These three 
rows of fibers, and the target that they face, are shown in Fig. l-31c). 


RI/RD82-285 


1-55 


ORIGINAL PAGE IS 
OF POOR QUALITY 



Liquid H 2 
p • 200 psia 
T - 40*R 


Sensor housing 


Lens 


Image of sensor on 
reflective target 


Thrust Load 


Window clamp 


Window with 
pressure seal 


Figure 1-30. Transient Axial Load Measurement Concept - HPFTP 
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Figure 1-31. R.ulial Bearing Load Measurement Concept, HPFTP 
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The probe Is shown installed in the housing in Fig. l-31a, and the cross-shaped 
target is shown attached to the bearing housing. The target must be raised 
above the surrounding material a sufficient amount (0.010 to 0.020 is suggested) 
to avoid reflections from the background. The purpose of using a proximity 
probe with transverse sensitivity is to allow Installation in an area with very 
limited radial access. The purpose of a dual probe is to reduce the number of 
separate probes to be installed. The probe is used here to sense the radial 
motion. 

A similar concept can be used to sense axial motion between stops of the turbine 
end bearing cartridge in the LOX pump, as illustrated in Fig. 1-32. This con- 
cept is not suitable for measuring loads greater than those causing bottoming 
of the cartridge against its axial stops. 

The fiber optic probe concepts described are capable of response that is suitable 
for hearing load monitoring. The accuracy and frequency response of fiber optic 
probes exceed the capabilities of eddy current or capacitive probes. 

Piezoelectric Load Cells . Some design concept effort was expended during 
November 1977 - March 1978 in which a radial load measurement concept for the 
turbine-end bearings of the LOX pump employing piezoelectric sensors was inves- 
tigated. The concept made use of piezoelectric load cells installed in the 
bearing cartridge, either in direct contact with the bearing outer race, or in 
contact with the next radial interface in the cartridge. The cells were to be 
typically 3/16- to 1/4-inch thick, and installed rigidly in the cartridge. The 
cells would sense only that part of the total bearing load that would pass 
through them, with the percentage of the bearing load sensed depending on the 
amount that the cells protrude beyond the supporting surface. Since the load 
cell surfaces must contact the bearing outer races, sliding contact would occur. 
In addition, the protrusion of the load cells would modify the radial clearance, 
possibly interfering with bearing operation and altering the turbopump's dynamic 
response. For these reasons, implementation of piezoelectric load cells for 
radial bearing load measurements is not recommended. 
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Figure 1-32. Axial Load Measurement Concept, HPOTP Turbine Bearing Cartridges 
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Bearing Diagnostics 


Methods used to detect Incipient bearing failure fall into two major categori s: 
vibration analysis and wear debris analysis. Because rocket engine turbopumps 
have large coolant flows that are not recirculated, wear debris analysis is not 
considered a practical option for bearing monitoring. On the SSME, a system of 
accelerometers mounted on the turbopump casings is used for this purpose. In 
addition to using peak acceleration magnitudes as redline shutdown criteria, 
signals from the accelerometers are analyzed after testing to detect operational 
anomalies including bearing faults. The example of vibration signature of 
engine operation shown in the PSD chart in Fig. 1-33 illustrates the synchronous 
speeds of both turbopumps and indications of rubbing. In addition to this gen- 
eral approach now used to bearing monitoring, some specialized techniques are 
available for detection of bearing defects. The following section presents 
characteristics of several bearing fault detection systems. 

Incipient Failure Detection (IFD) . A number of techniques have been developed 
to detect and assess the condition of rolling element bearings based on the 
dynamic signals issued from the bearing. Many of these are currently categorized 
under the general heading of incipient failure detection (IFD). The common 
aspects of most IFD methods are: 

• The detection is passive. All signals are generated within the bear- 
ing due to its own dynamic condition. 

• Condition judgment is done by comparative evaluation. Damage is 
detected by observing large changes in signal levels. 

• M'st sensors use piezoelectric elements either in the form of accelero- 
meters or crystal detectors. 

Users of IFD often disagree about the mechanisms involved, the comparative 
effectiveness of the different techniques and even the classification of generic 
terminology; e.g., acoustic emission, shock pulse, etc. 
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Figure 1-33. Vibration Signatures from Accelerometer on 
High Pressure Oxygen Turbopump 
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In this discussion, IFD techniques will be categorized in four broad groups 
based on the apparent nature of the signal detected from the bearing. These 
groups and their characteristics are presented as follows: 

• High-Frequency Resonance Technique (HFRT) - detects the high-frequency 
resonant vibrations of the bearing that occur as a result of the 
impacting (the detected signal is already amplified). Generally, only 
one of several resonant frequencies is selected and filtered out for 
analysis . 

• Shock Pulse Method (SPM) - uses a specially designed piezoelectric 
accelerometer to measure the mechanical impacts or shock pulses 
generated in the bearing as a result of defects. These shock pulses 
excite the accelerometer at its resonant frequency and thus amplify 
the signal. 

• Acoustic Emission Technique (AE) - detects and amplifies the traveling 
stress waves that occur when a race or ball surface irregularity is 
impacted. The signal is amplified by the resonant response of the 
transducer . 

• Low-Frequency Vibration Technique (LFVT) - uses nonresonant signals 
from the bearing generally ranging in frequency from sub rotational to 
harmonics of the highest bearing defect frequencies. 

High-Frequency Techniques . High-frequency diagnostic techniques appear to 
be applicable to the SSME turbopumps provided transducers can be installed in 
the right locations to pick up the bearing information. The principle problems 
presented are: (1) the cryogenic temperatures and the immersion in tne liquids, 

(2) the need for both small size yet robust construction, and (3) the require- 
ment for direct or near-direct contact with the bearing outer race. The third 
item may be the key issue because of the dead band or slight radial clearance 
that exists between the outei rcc^s and the carrier rings. 
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Sensor availability . Eight accelerometer manufacturers were contacted to 
determine what sensors were available that would be capable of operating in I.H^ 
or LO^ while meeting both the general size and frequency requirements. The 
salient features of the candidate models are given in Table 1-13. Most of 
these units would require some design modification for the application. 

Endevco clearly has the experience edge. Their model 7704 is being supplied for 
SSME monitoring and has a proven performance record at LH^ temperatures. In 
addition, model 2272 is rated to -450F and models 2221 F and 2271 A also have 
been shown to operate at LH^ temperatures. The model 2221 F has a toroidal 
design and is attached by a bolt through the center hole. The design appears 
very compact and particularly rugged. 

All of the companies contacted, with the exception of Kulite and Wilcoxon, had 
some experience in building accelerometers capable of operating in LN^ and felt 
there would be no serious drawback to modifying existing designs for 
temperatures. Several manufacturers suggested using a hardline integral cable 
and a mounting thread size of no less than 10/32. 

Acoustic emission sensors are very similar to high-frequency accelerometers, use 
essentially the same materials and, therefore, are also likely to be suitable 
for cryogenic use and should be available. 

Size considerations . Size plays an important role, as the sensor must be 
small enough to fit in the desired location in the turbopump. The sensor’s 
mass determines its natural frequency; generally, the smaller the mass, the 
higher the frequency. Conversely, units with small mass also have low 
sensitivities. 

For the HFRT, the accelerometer should have a resonance of greater than 50 kHz. 
This condition is generally met by units weighing 2 to 7 grams and having cylin- 
drical dimensions of 0.1 to 0.4 inch in diameter and 0.4 to 0.5 inch high, not 
including the connector or connecting cable. Smaller units could be used 
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provided their sensitivities are sutflclent. For the shock pulse method 
employing resonant frequencies In the 30 kHz to 35 kHz range, slightly larger 
units are applicable. These are in the 25-gram mass range and have dimensions of 
O.b inch in diameter by 0.9 inch high. 

Transducer mounting . The placement of the accelerometer may also be diffi- 
cult considering the very tight spaces around the bearings in the pumps and the 
need to piovide a direct path of transmission from bearing to sensor. 

Bearings in both pumps operate with dead bands — small but finite radial 
clearances between the outer race and its housing. These are necessary to 
permit the bearings to slide axially. The turbine end bearings in the oxygen 
pump are axially restrained by a soft spring that permits the balance piston to 
control the shaft axial position without ap ing high axial loads to the ball 
bearings. Signal intensity will be attenuated in passing through these dead 
bands. 

Another aspect of the dead band's presence is that bearing vibrations are 
transmitted to the carrier ring only during intervals when the outer race is 
forced haru against the carrier. This, in turn, occurs only over a narrow 
arc as the force vector passes the sensor location. 

Mounting a sensor directly on each outer race would produce optimum signal 
transmission. Since this is impractical, an alternate approach might employ 
a steel rod as a wave guide to connect the sensor outside of the carrier 
assembly to one of the outer races, as shown <n Fig. 1-34 and 1-35 applied to 
the HPOTP bearings. The accelerometer-rod assembly must be kept very light 
weight and spring loaded to follow the motion of the outer race. The effects of 
preload on the accelerometer operation would require assessment prior to imple- 
mentation of this concept. 

A more conventional approach is to mount the piezoelectric sensor on the 
carrier ring as close as possible to the bearing contact points. This would 
provide good transmission from the carrier to the accelerometer. 
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Figure 1-34. Accelerometer Mounting Concepts, HPOTP Pump F.nd 
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Provided sufficient signal is transmitted across the clearance, reasonable 
detection is possible. Both radial and axial accelerometer mountings are shown. 


Low-Frequency Vibration Techniques (LFVT) . This system is considere. 
least applicable to the SSME turbonurops for bearing fault diagnosis. While this 
technique works well for detecting and diagnosing many machinery anomalies such 
as unbalance, misalignment and others, serious problems arise when it is direc- 
ted toward rolling element bearing monitoring. These include: 


1. Vibrations indicating other phenomena usually obscure low-frequency 
bearing defects; for example, vibrations from the rotor transmitted 
through the bearing are not discriminated from those arising within 
the bearing. 

2. Bearing vibrations in the low-frequency range are very complex and 
require a great deal of experience for proper interpretation. Each 
analysis must be carefully done and requires a fair amount of time. 
Because of the signal comolexity, computerizing the diagnostics is 
extremely difficult. 

3. Because the technique relies heavily on changes in often unique 
spectral patterns, baseline profiles must be generated on each machine, 
often after each rebuild. 

4. Since small signal changes are important, transducers and instrumen- 
tation must have flat response and nonvarying gains and sensitivities. 

A high degree of consistency in data processing and analysis are 
fundamental requirements for successful application of LFVT. 

5. Vibration patterns often change significantly with relatively small 
changes in a machine operating parameter such as speed or pressure, 
possibly exceeding preset limits of response. 

6. Bearing defects often show up as sidebands or modulators of other 
frequencies. Signal analysis to separate defect indicators requires 
high resolution frequency analyzers and is time consuming. 


o 
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Shock Pulse Method . Shock Pulse Method (SPM) was utilized early in the 
SSME program to monitor the turbopump bearings. To achieve maximum effective- 
ness, the signal-sensing accelerometer must be coupled as directly as possible 
to the bearing outer race. Each cylindrical fit with clearance that exists 
between the bearing and the sensor greatly attenuates the signal. In complex 
structures like the SSME turbopumps, direct solid metal paths from bearings to 
accelerometer mounting locations generally do not exist. Applied to a simple 
structure, the HPOTP bearing tester, the SPM was found to provide indications 
of bearing deterioration. In the commercial SPM unit marketed by SKF Industries, 
the voltage representing shock pulse intensity is a logarithmic function, so 
that at high signal levels, a small change in voltage corresponds to a large 
increase in shock intensity. Interpretation of the shock level could be 
enhanced by modifying output to emphasize incremental signal changes at the high 
background level typical of rocket engine turbopumps. 

Optical Def lectometer . A fiber optic def lectometer can be incorporated into an 
active system that provides a means of monitoring the condition of individual 
bearing balls and races, and of relating bearing performance to operating loads 
and clearances. The condition of a ball bearing is monitored with the deflecto- 
meter by means of a fiber optic bundle mounted near the surface of the bearing 
outer race. A light beam directed at the race surface is reflected from it to 
a detector via the fiber optic bundle. Variations in light intensity provide 
a signal corresponding to the localized, cyclic deformations on the race caused 
by the passage of the bal* s (Fig. 1-36). 

The signals produced bv the rar’ surface motions contain rotor motion informa- 
tion as well as bearing condition indications. Figure 1-23 shows the potential 
application of a def lectometer to the HPFTP No. 2 bearing. The high pressure 
and relative inaccessibility of the turbine end of the HPFTP would preclude use 
of def lectometers of the same form at that location. Application to the HPOTP 
will require consideration of the oxygen compatibility of the probe materials. 
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TASK II: MULTIPLANE AND IN-HOUSING BALANCE SYSTEMS 


In Task I, the application of multiplane balancing to the SSME high-pressure turbo- 
pumps was analytically found to effect substantial reductions in bearing loads and 
shaft deflections. Task II addressed the means of performing multiplane balancing, 
including requirements for software and hardware, number of correction planes nec- 
essary to achieve acceptable balance, and physical requirements for in-housing bal- 
ance corrections. 

Conceptual and analytic studies conducted by MTI Included two major elements: 

1. Computer programs were adapted or developed to interrogate multiple sen- 
sors, cancel runout, filter noise, and separate different whirl motions. 
The required monitoring and data storage systems concepts were defined. 

2. Balance simulations were conducted to assess balance procedures of vary- 
ing complexity. The number of sensors and the number of speeds from 
which balancing data were derived, as well as the number of correction 
planes used, was varied to identify the least involved procedure yield- 
ing satisfactory results. 

Assessments were made of the physical requirements and consequences of implementing 
in-housing balancing including access, correction means, maintenance of cleanliness, 
and bearing 1 mage prevention. 

COMPUTER ANALYSIS OF SIGNALS 

This section presents the results of studies directed specifically toward in-hous- 
ing balance systems potentially applicable to the SSME high-pressure turbopumps. 
Software codes were implemented to interrogate multiple sensors, filter noise, can- 
cel runout, and separate different whirl motions. In addition, alternative concepts 
wore generated for integrating the software with computer hardware, sensors, and 
additional software to i: lement an in-housing balance system. 
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Software 


Software to perform the following functions has been implemented and tested: 

• Interrogate multiple sensors 

• Filter noise 

• Separate whirl motions 
® Cancel runout 

• Balance (linear least squares influence coefficient method) 

Each software component is described in subsequent paragraphs. 

Interrogate Multiple Sensors . The software to interrogate sensors is hip depend- 
ent upon the hardware chosen for the data acquisition process. For examp the 
software necessary to control the analog to digital converter must consider items 
such as: 

• The hardware register addresses 

• The metnod of handling device interrupts 

• The meaning and use of the various hardware control registers 

• The type of A/D triggering to be used 

A code is presented in Appendix F that will acquire data from multiple sensors; 
however, this code was developed specifically for a DEC PDP 11/03 microcomputer, 
the DEC RT-li operating system, and a DATA TRANSLATION analog to digital converter. 
While this software is highly dependent on the comptuer configuration, the soft- 
ware format illustrates the general seeps required to acquie data from multiple 
sensors. 

For example, the inputs to a generalized A/D sampling routine must include: 

« A/D multiplexer address code 

• Sampling rate 

• Number of desired samples 

• Mode control flags 


RI/RD82-285 


The output of the sampling routine is the set of sampled data points and a status 
flag that indicates whether any errors occurred during the sampling process. To 
enhance its generality, the sampling software must also convert each sampled data 
point from the A/D format, such as 12-bit offset binary format, to a more common 
format such as volts in a 32-bit floating point format compatable with high-level 
programming languages. 

Filter Noise . The filtering software implements low-pass, high-pass, and band-pass 
Butterworth-type digital filters. The software for each type of filter consists 
of a routine to design the filter and a routine that uses the calculated filter 
coefficients to filter a set of data. The routines and their functions are: 



LOPAS 

- Low-pass 

filter a set of i 

da«.a 


LPDES 

- Design a 

low- pass filter 



HIPAS 

- High-pass 

filter a set uf 

data 


HPDES 

- Design a 

high-pass filter 



BANPAS 

- Band-pass 

filter a set of 

data 


BPDES - 

- Design a band-pass filter 



The source code for each of these routines is contained in Appendix G. 

To test the software, a 512 point set of data was generated from the equation: 

Y(t) * 2*SIN(w*t) + SIN(3*w*t) + SIN(5*w*t). Figure 2-1 illustrates the resulting 
time series data set. The data were input to routine LOPAS together with the 
parameters that specified a four-pole filter with a corner frequency of 2*w. 

Figure 2-2 shows that only the 2*SIN(w*t) component of the data was passed by the 
low-pass filter. Next, the original data set was input to routine HIPAS together 
with the parameters that specified a four-pole filter with a corner frequency of 
4*w. Figur' - 2-3 shows that only the SIN(5*w*t) component of the data was passed 
by the high-pass filter. Finally, the data set was input to routine BANPAS 
together with the parameters that specified a four-pole filter with a pass band 
from 2*w to 4*w. Figure 2-4 shows that only the SIN(3*w*t) component of the data 
was passed by the band-pass filter. 
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Figure 2-1. Unfiltered Data 
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Figure 2-4. Band-Pass Filtered Data 



1 


Examination of Fig. 2-2 through 2-4 indicates that a phase shift has been intro- 
duced into the filtered data set. This phase shift Is expected. Additionally, 

Just as with hardware filters, software filters require a finite amount of time 
to initially "charge up." This explains the amplitude distortion near the begin- 
ning of the filtered data set. 

Separate Whirl Mo t ions . The software to separate whirl motions consists of a Fast 
Fourier Trarsorm (FFT) routine. The routine, entitled FFTR, accepts as input two 

[ ' It 

arrays (real and imaginary) of digital data, optionally applies a Hanning window 
function to the data, performs the FFT computation, and returns the FFT coefficients 
(real and imaginary) together with a scaling factor. The calculated FFT coeffici- 
ents overwrite the original data arrays. The FFT routine implements the Cooley- 
Tukey FFT algorithm using time decomposition and input bit reversal. The Hanning 
window used is defined by the equation: H(t) ■ 0.5 * (1 . -COS (w*t ) ) * F(t), where 

F(t ) is the original data set. A source code listing for routine FFTR is contained 

in Appendix H. 

To test the routine, the data set of Fig. 2-1 was input to routine FFTR. Figure 
2-5 shows the resulting spectrum. As expected, the data set was found to possess 
frequency components of w, 3w, and 5w at amplitudes of 2, 1, and 1 respectively. 

Cancel Runout . The software to cancel runout consists of the routine SUBRUN which 
accepts as input two spectra calculated by FFTR and subtacts one of the spectra 
from the other spectrum. The subtraction is performed one spectral line at a time. 

A source listing for routii e is contained in Appendix I. 

] 

To test routine SUBRUN, two time series data sets were generated as follows: 

Runout Speed: R(t) - 2*SIN(w*t) + SIN(3*w*t) + SIN(5*w*t) 

Balance Speed: B(t) - 5*SIN(w*t) + 0 5*SIN(2w*t) + 2 .5*SIN(3w*t) + 

1. 5*SIN(5w*t) 

Routine FFTR was used to compute the FFT spectrum for each data set. Figures 2-6 
and 2-7 show the spectrum for each data set. The two spectra were input to routine 
SUBRUN and subtracted. Figure 2-8 shows the resulting runout corrected spectrum. 
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Figure 2-7. balance Speed FFT Spectrum 





Balance . Program BALANCE Is a mult i speed, multiplane balancing program that 
implements the linear least-squares influence coefficient method of balancing. 

The inputs to program BALANCE include: 

• Number of balance speeds 

• Number of balance planes 

• Number of balance probes 

• Vibration data at runout speed 

• Vibration data at balance speed 

• Magnitude and location of installed trial weights 

• Trial weight vibration data sets 

The outputs from program BALANCE include: 

• Influence coefficient matrix 

• Balance weight set 

• Predicted residual vibration after balancing 

The software consists of a main program and a routine to invert a complex matrix. 

The main program prompts for the inputs, performs the balancing calculations, and 
outputs the results. The routine CMATEQ is invoked by the main program and per- 
forms the inversion of the influence coefficient matrix. Appendix J contains list- 
ings of both routines. 

Figure 2-9 is a listing of a sample run of program BA1.ANCE. The data input to the 
program represents actual vibration daca acquired from an MTI balancing test rig. 

The results output by program BALANCE agree with the results of other in-house 
balancing tools. 

System Integration Concepts 

To create a fully functional in-housing balancing system, the software components 
described above must be integrated with computer hardware, other software components. 
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u 

ENTER NUMBER OF BALANCE SPEEDS I 1 

ENTER NUMBER OF BALANCE PLANES: 2 

ENTER NUMBER OF PROBES: _J_ 

•OUT OF ROUNDNESS DATA* 

PROBE 1 (RADIUS, PHASE): .2606.248.2 

PROBE 2 (RADIUS, PHASE): .4687,246.3 

•BASELINE DATA* 

**SPEED NUMBER 1** 

PROBE 1 (RADIUS, PHASE, WEIGHTING FACTOR): .6718.84.7.1.0 

PROBE 2 (RADIUS, PHASE, WEIGHTING FACTOR): 1.5801.98.5.1.0 

•TRIAL WEIGHT DATA* 

••PLANE NUMBER 1** 

ENTER TRIAL WEIGHT (MAGNITUDE, PHASE): 1.55.0 

ENTER VIBRATION DATA (SPEED NUMBER 1) 

PROBE 1 (RADIUS, PHASE): .2458,50.4 

PROBE 2 (RADIUS, PHASE): 1.6650.167.9 

••PLANE NUMBER 2** 

ENTER TRIAL WEIGHT (MAGNITUDE, PHASE): 1.55,0 

ENTER VIBRATION DATA (SPEED NUMBER 1) 

PROBE 1 (RADIUS, PHASE): .5812,249.6 

P .JBE 2 (RADIUS, PHASE): .9836,255.6 

•INFLUENCE COEFFICIENT MATRIX* 

••SPEED NUMBER 1** 

( 6. 105E-02, -3.094E-01) (-1.707E-01.-7.830E-01) 

(-8 . 996E-01 ,-7 . 831E-01) (-7 . 134E-03 ,-l . 623E+00) 

•CALCULATED VALUES OF CORRECTION WEIGHTS* 

PLANE NUMBER - 1 WEIGHT - 0.106 ANGLE - 68.5 

PLANE NUMBER - 2 WEIGHT - 1.154 ANGLE - 0.3 

•PREDICTED RESPONSE AFTER BALANCING* 

••SPEED NUMBER 1** 

PROBE AMPLITUDE PHASE 

1 0.261 248.2 

2 0.469 246.3 

Figure 2-9. Program Balance Weighted Least Squares Balancing Program 
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and the sensors that measure the vibration levels on the unit to be balanced. MTI 
has generated two alternative concepts for an integrated system: 

• A microcomputer-based single user system 

• A minicomputer-based multiuser system 

These two generic system concepts are discussed in subsequent paragraphs froa both 
a hardware and software configuration perspective. 

Hardware Configuration . The microprocessor-based single-user balancing system is 
composed of the following hardware components: 

• 16-bit microprocessor 

• Memory 

• Mass storage device 

• Auxiliary mathematics processor 

• A/D converter 

• Programmable real-time clock 

• CRT terminal 

• Other output devices 

• Signal conditioning electronics 

The single user system is configured around a 16-bit microprocessor to provide the 
•cmputat ional power necessary to perform the FFT and balancing computations and to 
provide reasonable system speed. Sixty-four kilobytes of main memory will be suf- 
ficient to meet the needs of the system. The mass storage system will be used to 
store the computer operating system software, the balancing system applications 
software, and the acquired sensor data. At least 1 megabyte of disk storage will 
be necessary. A dual-drive disk system with at least one removable media is recom- 
mended to provide the capability to back up the acquired data or to transfer data 
to other computer systems. The auxiliary math processor provides the hardware 
floating point arithmetic support that will enhance the system's computational cap- 
abilities and speed. The A/D converter is used to digitize the analog sensor sig- 
nals. At least an 8-channel, 12-bit, 100 kHz throughput rate A/D converter is 
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recommended. The p" ogrammable real-time clock is used to acquire speed data and 
also is used to provide the a/D with a sampling trigger signal. The CRT terminal 
is used as the operator's control console. Other output devices might include a 
system printer or a graphics device. The signal-conditioning electronics hardware 
is used to condition the sensor output signals to make them compatible with the 
input requirements of the A/D converter. Suitable anti-alias filtering is a 
required part of the signal conditioning. Figure 2-10 illustrates the configura- 
tion of the single user balancing system. 

The multiuser balancing system is configured around a general-purpose, 16-bit mini- 
computer. This system is composed of many of the same types of hardware components 
as the single user system, including: 

• Central processor 

• Memory 

• Mass storage device 

• Floating point processor 

• Line printer 

• Graphics output devices 

A larger, more powerful processor is used. The system possesses more memory to 
support multiple users and a timesharing operating system. More mass storage is 
provided (at least 10 megabytes) and fixed media, Winchester-type drives are used. 

A floating-point processor is included to provide sufficient computational power. 
The line printer and graphics devices supply the necessary output resources. 

To make the system capable of supporting multiple balancing users, the following 
additional components are added: 

• Several CRT terminals 

• Several signal acquisition processors 

• Several signal conditioning electronics packages 

• An optional array processor 
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Since each balancing user has a separate CRT control console, a separate signal 
acqulsiton processor, and a separate signal conditioning electronics package, many 
users may use the balancing system simultaneously. The signal acquisition proces- 
sor may be a simple A/D converter or may be an Intelligent device with its own 
microprocessor. An array processor can be added to provide more computational 
power for such tasks as FFT calculation and matrix inversion. Figure 2-11 illus- 
trates the multiuser balancing system configuration. 

Software Configuration . The software for the balancing system is composed of three 
major components: 

• The computer operating system software 

• The balancing system applications software 

• Other software support packages 

The operating system software for the microprocessor-based balancing system would 
be a simple, single user software system. One control console is supported. 
Applications programs would be stored on disk and could be called up for execution 
one at a time. The operating system for the multiuser balancing system would sup- 
port many control consoles and would use a timesharing algorithm that would allow 
many applications programs to be executed at the same time. A priority-based 
scheduling algorithm would be used to ensure that time-critical tasks were given 
preference for execution time. Applications program development could be performed 
concurrently with balancing system operation. 

The balancing applications software package would be very similar for the two sys- 
tems, but the method of implementation of the software might vary. The functions 
supported by the software would include: 

• Data acqulsiton 

• Data stc .age 

• Filtering 

• r FT calculation 

• Runout cancellation 
e Balancing 
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• Report generation 

• Graphic output 

Figure 2-12 illustrates the various software modules and the flow of data among 
those modules. 

Generally speaking, the only other software support packages required would be the 
software to support the graphics output devices configured into the system. If 
desired, an intercomputer system communications software package could be added to 
provide the capability to link the balancing system to other computer systems. 

BALANCING SIMULATIONS 



Utilizing the response characteristics from the results of Task I, a study of 
multiplane balancing effectiveness was conducted at MTI. This study evaluated the 
applicability of linear influence coefficient - based balancing tools for the SSME 
fuel and oxidizer turbopump balancing processes. Since a tabrication-associated 
balancing process is being considered, the "dry" rotor response in the absence of 
working fluids was used. Specific factots considered included: 

• Balance planes 

• Probe location 

• Balance speed (s) 

Referring to Fig. 1-13 and 1-14, both the HPFTP and the HPOTP rotors are seen to 
have easily identifiable balance planes at the impeller and turbine disks. How- 
ever, each of the locations noted may not necessarily be needed and only the out- 
board components are accessible. Of prime concern in performing in-housing balance 
plane and probe selection is that the location be both responsive and accessible. 
Similarly, balance speed selection must take into account the rotor response maxi- 
mum operating speed expected and its relation to any resonant speeds. 
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The balance simulation studies showed that: 


1. The balancing code correctly identified the unbalance distribution con- 
tained in the taped response to worst-case unbalance distributions. 

2. Both HPFTP and HPOTP rotors display satisfactory response over the entire 
operating speed range when corrected in accordance with a multiplane 
balancing program based on low-speed response data. The significance of 
this result is that an improvement in balance quality can be obtained by 
multiplane balancing contrasted with sequential two-plane balance pro- 
cedures that may result in residual moment unbalances. The final turbo- 
pump assembly balance, therefore, is based ot. a mere ideal set of correc- 
tions and is expected to have reduced dynamic response. 

3. Effective in-housing balancing can be obtained from a limited number of 
correction planes, sensors, and balancing runs. In general, better 
balance condition quality is obtained as each of these factors is increased. 
Low-intensity response can be attained for the high-pressure turbopumps 
with the minimum impact balancing procedures listed in Table 2-3. The 
term, "minimum impact," describes the least complicated balance procedure 
capable of producing low-intensity response and is not intended to imply 
applicability to the present SSME high-pressure turbopump configurations. 

The practical constraints involved in gaining access to the in-board 
masses are not considered in this discussion. 

Unbalance Input 

As input to the balance simulation study an array of unbalances was selected that 
approximates the worst-case distribution identified in Task I (Tables 1-5 and 1-7). 

The phase angles of the unbalances were shifted slightly from the zero and 180- 
degree values of the worst-case distributions to provide a better assessment of 
the balance system's accuracy. The theoretical responses to the distributions 
shown in Table 2-1 were recorded on magnetic tape by Rocketdyne and transmitted to 
MTI. The taped response was the only definition of the unbalance input furnished 
to MTI. 

Studies Conducted 


Rocketdyne-supplied linear influence coefficients and the response to the unknown 
unbalance distribution was reviewed, and candidate balancing speeds and planes 
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TABLE 2-1. BALANCE SIMULATION STUDY-UNBALANCE DISTRIBUTION 





UNBALANCE AT PLANE 








SUBJECT 

TURBOPUMP 


1 

2 

3 

4 

5 

HPFTP 

MAGNITUDE, g-IN. 

4 

7 

7 

5 

5 


PHASE ANGLE, 
DEGREES 

37 

353 

4 

332 

18 

HPOTP 

MAGNITUDE, g-IN. 

2 

5 

5 

5 



PHASE ANGLE, 
DEGREES 

15 

156 

351 

22 



selected. A series of balancing trials was run using a linear least-squared 
influence coefficient balancing method. Based on the rotor mode shapes, as in 
Fig. 2-13, the most responsive high-pressure fuel turbopump (HPFTP) balancing 
planes are the third-stage impeller (III) and the first-stage impeller (I) for the 
14,000 and 29,000 rpm modes, respectively. The first turbine stage, however, is 
itself very responsive to the second mode and is accessible, while the first stage 
impeller is not, making the turbine a more desirable location for balancing. For 
the high-pressure oxidizer turbopurap (HPOTP), the most responsive planes for the 
two modes at 11,028 and 34,379 rpm are the first turbine stage and the main impeller, 
respectively. While the turbine stage is readily accessible, the same is not true 
for the main impeller. Thus, an alternate possibility, the preburner impeller, 
was considered since it also is responsive at the second mode (refer to Fig. 1-10 
which shows HPOTP mode shapes). This second balance plane was included to assess 
its applicability for contemplated increases in operational speeds beyond the 
current 30,000 rpm. Table 2-2 lists the conditions studied in the balance simula- 
tions, while the responses obtained are presented in Fig. 2-13 through 2-61. 

H PFTP Ba lancing Results . Even though 14,000 and 29,000 rpm were identified as the 
optimal balancing speeds for the HPFTP, several alternate yet still responsive 
speeds were considered. These alternate speeds were included to identify the low- 
est possible speeds for balancing, given the fact that in-housing balancing is to 
be performed under poorly or possibly even unlubricated corn... ions. As can be seen 
in Fig. 2-15 through 2-20 the HPFTP rotor can be balanced at 6000 rpm provided all 
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TABLE 2-2. BALANCE SIMULATION STUDY 


ORIGINAL PAGc 
OF POOR QUAL 


13 

Itv 



CO 

SI 

1 

c 

2-15 TO 2-20 
2-21 TO 2-25 
2-26 TO 2-30 
2-31 TO 2-34 
2-36 TO 2-39 

2-40 TO 2-43 
2-44 TO 2-48 
2-44 TO 2-48 
2-49 TO 2-53 
2-54 TO 2-58 
2-59 TO 2-63 

BALANCING PLANE CORRECTIONS 


o 

s 

2 s s a a 

N *- p- 


G-IN. 

4.6 

45.7 
5.0 
2.32 

19.8 



o 

3 

s 2 

jb 3 £ 2 3 8 

z 

1 

o 

7.3 

4.8 

8.87 

4.9 

4.99 

8.69 

9.06 

8.84 


DEG 

8 CO CO O' 

174 

171 

z 

1 

o 

5.9 

7.2 

22.2 

34.9 

CM r— 

© O 

vO 60 


DEG 

167 

173 

337 

336 

z 

1 

o 

~ 00 

CO vO 

4.99 

5.00 


o 

UJ 

Q 

222 

217 

196 

195 

179 

320 

315 

Z 

1 

o 

66" E 
86E 

2.01 

2.00 

1.77 

5.51 

3.96 

PROBE/STATION NUMBER 

(NJ (M M M N 

o o 

CO 00 

>- 

O CO VO CO 00 

9 

1 IX 
9 
9 
9 
9 


>» 

o> 


X 

a* 

BALANCING 
SPEEDS. 
RPM X 1000 

6 

6 

14 29 

14 22 29 
6 14 

7 

7 

14 35 

6 

14 34 

14 35 


NUMBER 

OF 

PLANES 

vo vo cm cm 

- vt 9 CM CVJ M 


NUMBER 

OF 

PROBES 

o cn vo co o') 

CM ** CM CM CM CM 


NUMBER 

OF 

SPEEDS 

*— *— CM CO CM 

r- r- CM r- CM CM 

o. 

1 

HYDROGEN 

OXYGEN 


RI/RD82-285 



HPfTP -UNCORRECTED ROTOR RESPONSE 



Figure 2-13. HPFTP Uncorrected Rotor Response, Worst Case Unbalance 



HPOTP UNCOBBECTED BOiOB BESPONSE 



Figure 2-14. HPOTP Uncorrected Rotor Response, Worst Case Unbalance 





Figure 2-15. HPFTP - 1-Speed, 6-Probe, 5-Plane Balance (Probe 7) 






AFTER BALANCE 



Figure 2-17. HPFTP - 1-Speed, 6-Probe, 5-Plane Balance (Probe 13) 









Figure 2-18. HPFTP - 1-Soeed, 6-Probe, 5— Plane Balance (14,000 rpm Mode) 
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available balance planes are used. Due to the distributed nature of the 
correction weights, this low speed balance will provide for smooth operation 
throughout the entire HPFTP speed range under the assumed operating conditions. 

A second balancing configuration where the number of speeds, probes, and balance 
plane combinations were reduced to a minimum of one plane and one speed was 
evaluated for acceptability. As can be seen in Fig. 2-21 through 2-25, using 
plane V at 6,000 rpm with three probes, does not provide an acceptable balance 
at operating speeds. The postbalance results are worse than prebalance response 
and, in many instances, exceed acceptable seal clearances over the length of 
the rotor. 

Based on the results obtained and shown in Fig. 2-26 through 2-30, it is evident 
that the use of all five balance planes also will provide a high-quality balance 
when data are collected at high speeds. In this third case, all five balance 
planes were included as was the influence coefficient data for balancing speeds 
of 14,000 and 29,000 rpm. In addition to determining the optimum balance attain- 
able using all five available balance planes, the results of this case show that 
the tape-encoded unbalance distribution was accurately detected. The balance 
corrections were: 3.99 g-in. at 37 degrees, 6.87 g-in. at 357 degrees, 7.2 g-in. 

at 4 degrees, 4.8 g-in. at 330 degrees, and 5.0 g-in. at 18 degrees, respectively, 
for balance planes I through V; these values agree closely with those of Table 2-1. 

The remaining two cases are presented to show the minimum balance plane and 
speed combinations necessary to balance for smooth operation over the operating 
speed envelope. 

As can be seen in Fig. 2-31 through 2-34, when data are collected at both the 
first and second critical speeds (14,000 and 29,000 rpm), and a speed between 
the two criticals that appeared responsive (22,000 rpm), acceptable rotor 
response can be attained over the rotor's operating speed envelope through the 
use of only two balance planes (III and V). 
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Figure 2-22. HP FTP Rotor, - 1-Speed, 3-Probe, 1-Plane Ealance (Probe 13X) 
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Figure 2-24. HPFTP Rotor, - 1-Speed, 3-Probe, 1-Plane Balance (29,000 rpm) 
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Figure 2-26. HPFTP - 2-Speed, 6-Probe, 5-Plane Balance (Probe 7X) 
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Figure 2-27. HP FTP - 2-Speed, 6-Probe, 5-Plane Balance (Probe 13X) 
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Figure 2-30. HPFTP - 2-Speed, 6-Probe, 5-Plane Balance (38,000 rpm) 
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Figure 2-32. HP FTP Probe 13 (3-Speed, 3-Proh« Planes III and V) 
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Figure 2-33. HPFTP Prebalance Response 
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Figure 2-34 
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Having determined that balance planes III and V can be used to control rotor 
unbalance response when data are acquired at the first and second critical 
speeds, an attempt was made to see if appropriate balance weights for these two 
planes could be calculated from data acquired at lower speeds. The predicted 
correction weights, as presented in Table 2-2, were determined from influence 
coefficient data generated for 6,000 and 14,000 rpm. Figures 2-35 through 2-39 
show that this approach was unsuccessful. While the response at speeds below 
14,000 rpm are improved, at speeds above 17,000 rpm, the response is generally 
worse than for the uncorrected rotor with vibration amplitude reaching a peak 
of 15 mils at approximately 29,000 rpm (Fig. 2-38). 

Based on this preliminary linear evaluation, it is evident that a minimum of two 
balancing speeds (near 14,000 and 29,000 rpm) and two independent balancing 
planes (III and V) are required to obtain minimum HPFTP rotor response within 
the operating speed envelope. As stated previously, both balance planes I and V 
are responsive and either can be used to control vibrations at the second 
critical speed. However, balance plane V (the first turbine stage) was selected 
due to its accessibility for balancing. Attempts at using balance planes IV and 
V to control the first and second modes proved unsuccessful due to the close 
proximity of these planes to one another. This close proximity results in a 
lack of independent or redundancy of balance planes. Therefore, the third-stage 
impeller (plane III), which is active during the first mode, is required if 
high-speed balancing is to be effective. 

HPOTP balancing results . Figures 2-40 through 2-63 present the results of the 
HPOTP linear balancing study. For current maximum operating speeds, a single- 
plane, low-speed balance using the first-stage turbine shoulc provide acceptable 
response, as shown in Fig. 2-40 through 2-44. However, if operating speeds 
increase beyond 33,000 rpm, vibrations begin to rise rapidly as the second 
critical speed at 34,379 rpm is approached. Since at the second critical speed, 
maximum displacement occurs near the main impeller, a second balance plane is 
needed to control this mode. 
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Figure 2-36. HPFTP - 2-Speed, 3-Probe, 2-Plane Balance (Probe 13) 
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Figure 2-37. HPFTP - 2-Speed, 3-Probe, 2-Plane Balance (14,000 RPM) 
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Figure 2-38. HP FTP - 2-Speed, 3-Probe, 2-Plane Balance (29,000 RPM) 
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Figure 2-39. HP FTP - 2-Speed, 3-Probe, 2-Plane Balance (38,000 RPM) 
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Figure 2-40. HPOTP - 1-Speed, 2-Probe, 1-Plane Balance (Probe 6Y) 
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Figure 2-42. HPOTP - 1-Speed, 2-Probe, 1-Plane Balance (13,000 RPM) 










To determine the best balance quality possible using all balance planes, two 
trials were conducted. For the first case, the balance speed selected was 

7.000 rpm while. In the second case, 14,000 and 35,000 rpm (the first and second 
critical speeds) were chosen. As seen in Table 2-2, the results of balancing in 
each case are almost Identical. Predicted response throughout the speed range 
is shown in Fig. 2-45 through 2-48. In addition to giving an indication of the 
quality of balance attainable when using all available balance planes, these two 
test cases also were used to determine the total distributed unbalance used to 
predict the rotor response for this balancing study. As seen in Table 2-2, the 
single-speed and two-speed cases yield essentially the same results for required 
correction weights. 

The required balance plane/sneed combination for operation to speeds up to 

35.000 rpm were illustrated in the final three test cases. For operation at 
speeds up to 35,000 rpm, a low-speed (6,000 rnm) , two-plane (preburner and first 
turbine) balance does not provide acceptable response. Vibrations remain the 
same both before and after balancing at 35,000 rpm (Fie. 2-49 through 2-53). 
These influence coefficients do not take into ac> _>unt the effects of the second 
mode and, thus, are ineffective in reducing the rotor response to acceptable 
levels above 32,000 rpm, when only two planes are used for low-speed balancing. 
Therefore, to ensure that both modes are taken into account, influence coeffi- 
cient data calculated for 14,000 and 34,000 rpm were used to predict the two- 
plane correction weights and residual response after application of calculated 
correction weights. As can be seen in Fig. 2-54 through 2-58, rotor response is 
generally improved in the operating speed envelope as well as through speeds up 
to approximately 43,000 rpm. However, it should be noted that at 30,000 rpm, 
preburner impeller response increased by 1.1 mil and that response at 35,000 rpm 
is still quite high. This is due, in part, to the fact that the 34,000 rpm 
balance speed is not exactly at the critical speed where optimum balance data 
are obtained. Therefore, the response will not be sufficiently reduced when 
passing through the critical speed. For this reason, a test case in which the 
balancing speeds were changed to 14,000 and 35,000 rpm was included. 
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Figure 2-45. HPOTP - 1-Speed, 4-Probe, 4-Plane Balance (Probe 14X) 





























Figure 2-48. HPOTP - 1-Speed, 4-Probe, 4-Plane Balance (30,000 RPM) 
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Figure 2-52. HPOTP - 1-Speed, 2-Probe, 2-Plane Balance (30,000 RPM Mode) 
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Figure 2-54. HPOTP - 2-Speed, 2-Probe, 2-Plane Balance (Probe 6X) 
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Figure 2-55. HPOTP - 2-Speed, 2-Probe, 2-Plane Balance (Probe 14X) 
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Figure 2-56. HPOTP - 2-Speed, 2-Probe, 2-Plane Balance (20,000 RPM Mode) 
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Figure 2-57. HPOTP - 2-Speed, 2-Probe, 2-Plane Balance (30,000 RPM Mode) 
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Figure 2-58. HPOTP - 2-Speed, Probe, 2-Plane Balance (34,000 RPM Mode) 
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Figures 2-59 through 2-64 show that with this change in the second balancing 
speed, rotor response is considerably improved at 35,000 rpm without sacrificing 
lower speed operation. This does point out, however, the need for accurate data 
acquisition and proper speed selection when balancing for operation above or 
near a critical speed. 


Minimum Impact In-Housing Balance . The in-housing balance procedures producing 
satisfactory responses, and which involve the minimum number of balance runs and 
correction planes, are summarized in Table 2-3 for the HPFTP and HPOTP. Access 
to the HPFTP third-stage impeller is required for this procedure, while 
corrections appl ?d to externally accessible planes will suffice for the HPOTP. 


TABLE 2- j. MINIMUM IMPACT IN-HOUSING BALANCING REQUIREMENTS 

FOR LOW-AMPLITUDE RESPONSE 



BALANCE SPEEDS, 

CORRECTION 

DATA SENSOR 

RESPONSE SHOWN 

UNIT 

RPM 

PLANES 

LOCATIONS* 

(FIG. NO.) 

HPFTP 

14000 

III 

F ! 

2-31 


AND 

AND 

THROUGH 


29000 

V 

F 3 

2-34 




F 4 


HPOTP 

14000 

I 

°o 

2-59 


AND 

AND 

L 

THROUGH 


35000 

V 

°4 

2-64 


♦Refer to Fig. 1-25 (HPFTP), 1-26 (HPOTP). 


RI/RD82-285 

2-74 




SPEFO (*1000 8PM) 


Figure 2-59. HPOTP - 2-Speed, 2-Probe, 2-Plane Balance 
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Figure 2-60. HPOTP - 2-Speed, 2-Probe, 2-Plane Balance 
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Figure 2-61. HPOTP - Deflected Rotor Shape (20,000 RPM) 
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Bearing Loads . In the Balance Simulation Studies, rotor response results were 
obtained in terms of radial displacements. For purposes of comparison only, 
estimates of load t eductions obtained by the minimum impact corrections for the 
worst-case unbalances were made. These load values are the products of the 
deflections at the bearing stations and a nominal bearing stiffness of 0.5 x 
10^ lb/in. The precorrection loads listed in Table 2-4 are, in most instances, 
substantially reduced by the balance procedure as shown in Table 2-5, which 
presents the post-balance load estimates. The loads listed should be used for 
comparison only because of the approximations used in their calculation. In 
addition, the following factors should be applied when considering these 
load values: 

1. The dry response applies only to a balancing procedure. 

2. No dead-band effects are included, If these clearances were present, 
the response may be attenuated as they will absorb some of the rotor 
radial motions without loading the bearing. 

IN-HOUSING CORRECTION METHODS 


Performing balance weight corrections to an assembled turbopump requires access 
to the correctioii plane, a means of modifying the mass distribution of the 
balance plane, and a means of maintaining cleanliness of the flow passage. 

Access 


Accessibility of turbopump rotor components diminishes as the assembly progresses 
toward final installation on •'he engine. At the turbopump component level, 
access is feasible to the outboard turbine disk and the inlet of the outboard 
pump impeller. It is at this stage that the present in-housing trim balance is 
performed. Access to inboard rotor components or access to any rotor components 
at a later stage of assembly requires an entry port positioned properly with 
respect to an area on the rotating member that has been adapted to weight change. 
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TABLE 2-4. PREBALANCE RESPONSE 


UNIT 

SPEED. 

RPM 


BEARIN 

3 POSITIO 

H 

1 

2 

3 

4 

HPFTP 


STATION 

1 

2 

12 

13 


14,000 

DISPLACEMENT, MILS 

5.5 

8.5 

22 

20 



LOAD, POUNDS 

2750 

4250 

11000 



22,000 

DISPLACEMENT, MILS 

0.72 

0.78 

0.22 

0.18 



LOAD, POUNDS 

360 

390 

110 

90 


29,000 

DISPLACEMENT, MILS 

1.8 

1.8 

2.2 

2.2 



LOAD, POUNDS 

900 

900 

1100 

1100 


38,000 

DISPLACEMENT, MILS 




0.8 



LOAD, POUNDS 




400 

HPOTP 


STATION 

3 

4 

9 

10 


20,000 

DISPLACEMENT, MILS 

0.4 

0.25 

0.55 

0.6 



LOAD, POUNDS 

200 

125 

275 

300 


30,000 

DISPLACEMENT, MILS 

0.2 

0.25 

0.5 

0.5 



LOAD, POUNDS 

100 

125 

250 

250 


34,000 

DISPLACEMENT, MILS 

3.8 

4.2 

3.0 

2.5 



LOAD, POUNDS 

1900 

2100 

1500 

1250 
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UNIT 

SPEED, 

RPM 

HPFTP 

14,000 


22,000 


29,000 


38,000 

HPOTP 

20,000 


30,000 


34,000 


T ABLE 2-5. MAXIMUM BEARING LOADS - POST-BALANCE RESPONSE 


BEARING POSITION 


STATION 


DISPLACEMENT, MILS 
LOAD, POUNDS 

DISPLACEMENT, MILS 
LOAD, POUNDS 

DISPLACEMENT MILS 
LOAD, POUNDS 

DISPLACEMENT, MILS 
LOAD, POUNDS 


STATION 

DISPLACEMENT, MILS 
LOAD, POUNDS 

DISPLACEMENT, MILS 
LOAD, POUNDS 

DISPLACEMENT, MILS 
LOAD, POUNDS 
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A single-wall casing can be penetrated for correction by providing a cover or 
threaded plug, as shown in Fig. 2-64, which depicts an adaption of a Mark 3 fuel 
pump to In-housing balance. Penetration of the casings of the SSME high-pressure 
pumps, however. Is undesirable because raising stress levels in pump casings 
lowers the rupture safety factor, and the access port Increases the potential for 
propellant leakage. Access to inboard impellers also requires crossing more than 
one wall or parting plane, adding to the risk of leakage associated with the 
access port. 

Co rrect ion M ethods 

The correction methods can be divided into several categories: weight exchange, 

material removal by means of mechanical processes, solutions, electrical dis- 
charge, and laser beams. 

Exchange able Weights. Exchangeable weights are used where cutting of the rotor 
is detrimental, as in the case of the gold-plated SSME turbine wheels. In the 
HPOTP, bolts of differing weights are installed to achieve correction, while 
washers are inserted under turbine wheel securing nuts in the HPFTP. If repeated 
and varying balances are to be performed, as in spin arbors, socket drive head- 
less set screws are often inserted, as needed, in a circle of threaded holes in 
the arbor. Since structural material is not removed, weight exchange has the 
advantage of low contamination potential; the risks involved include weight 
migration and dlslodgement , loss of weights that may fall into the casing during 
installation or removal, slivers of material stripped from screw threads, and 
contamination introduced through the access port on the correction weight or on 
the installation tool. A concept for set screw weight correction of the HPOTP 
preburner impeller is shown in Fig. 2-65. An alternate method, shown in Fig. 2-66, 
utilizes captured pins in the hub area. 

A high-priority requirement of any weight addition correction process is the 
security of the weight's fastening. This is generally provided by captured 
locking pins or, in the case of set screws, bvittons of material enclosed in 
recesses in threaded weights that prevent screw rotation by friction. Loss of 
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Figure 2-66. Preburner Impeller Exchangeable Weight Captured Pin 






weight- at installation or removal can be prevented by insertion of a guard tube 
through the access port to contain any dropped weights. 

Although not strictly exchangeable, corrections can be made by welding weights 
to rotors. This technique is most useful when a large correction must be made, 
and, since fine adjustments are difficult to make in this manner, a subsequent 
balance procedure is a necessity. 

Movable Weights . Weights whose position can be changed to affect balance 
include headless screws in threaded radial holes in the rotor; the balance is 
changed by varying the weights' ef fective radius. An alternate method requires 
two nesting annular masses with an eccentric interface. Changing the relative 
angular orientation of the masses alters their effective eccentricity. 

Mechani cal Material Removal. Material removal is probably the safest correction 
method, since the weight adjustment cannot change in amount or location. Balance 
adjustment is often effected by removal of material from an appropriate location 
on the rotor by rotary file, milling cutter, or twist drill. 

This method is well adapted to out-of-housing balancing since the component can 
be effectively cleaned to remove the chips generated. When applied to in-housing 
correction, however, the debris generated by mechanical removal of material 
poses an almost insoluble contamination problem, rendering this methoa unappli- 
cable to the pump impellers of the SSMF turbopumps. Even though a vacuum 
collector could be incorporated into the cutter mechanism, it is impossible to 
guarantee that no chips escape and enter the propellant passages. Application 
of mechanical removal means to turbine areas may be acceptable, since propellant- 
clean conditions may be less important in these areas. 


R1/RD82-285 

2-87 


Solution-carried Material Removal . Electrochemical machining (ECM) adapted to 
material removal for balancing is shown in concept in Fig. 2-67. Advantages of 
this system are: 

1. There are no moving parts at the machining site. 

2. The cavity formed is smooth-surfaced and spherical in shape. 

3. The amount of material removed is accurately controllable since it is 
a function of the product of the electric current and the time of 
exposure. 

4. The removed material leaves the work site in solution, and thus con- 
tains no particulate matter. 

5. The work area is never opened to the pump interior when working fluid 
is present. The removal process can be fol loved b” a flush of the work 
area with a neutral fluid to remove all traces of electrolyte. 

For ECM, the working fluid must be conductive. Representative electrolytes a.e 
KOH (potassium hydroxide) and NaOH (sodium hydroxide); the specific fluid 
selected will depend upon the material to be removed. An electric potential of 
8 to 15 Vdc between the workpiece and the fli *d delivery channel is required. 

An alternate material removal process, chemical milling, requires a chemically 
active working fluid such as HE (hydrofluoric acid) or HNO^ (nitric acid), but 
requires no electric potential. Use of these strong fluids will require 
neutralization and flushing. 

Electrical Discharge Machining (EDM) . This process uses electric discharge 
through a nonconductive (dielectric) fluid to melt small globules from the work 
piece. These particles are carried away from the work site by the circulating 
fluid, which may be oil or deionized water. The fluid enters the work area 
through a tube, as shown in the concept in Fig. 2-68. 
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SPHERICAL CAVITY FORMED 






Figure 2-68. EDM Pat tide Removal 



There is an electric potential difference of approximately 100 V between the 
tubular electrode and the work piece. The gap between the electrode and the 
work piece is maintained at approximately 0.001 inch by a servomechanism that 
advances the electrode as material is removed. 

As in the case of ECM, cleanliness can be maintained by flushing the flow path 
with propellant-compatible fluid, then drying the passage with dry, inert gas 
prior to disengaging the sealed fitting from the work piece surface. 

Material Removal By Laser . Balance correction based on material removal by 
laser beams has been utilized experimentally (Ref. 3-1). The energy of the beam 
is sufficient j locally melt and/or vaporize the work piece material. A poten- 
tial advantage of this correction method is that it can be automated to perform 
corrections while the rotor rotates if the laser discharges are timed to impinge 
repeatedly on the desired area. Disadvantages of laser correction lie in the 
small amount of material removed per shot, and the potential contamination of 
the machine interior by melted globules of rotor material. The following dis- 
cussion and assessment of the application of balance correction by laser to the 
SSME high-pressure turbopumps was made by MTI . 

The rate of material removal attainable with a laser system is dependent on the 

beam energy and the work piece area exposed to the beam pulse. A laser system 

-3 

with an output of 25 J and a pulse duration of 1 x 10 seconds has a corre- 
sponding peak power of 25,000 W. A typical beam divergence ot this type of 
system is less than 15 mrad. If a lens of 4-inch (10 cm) focal length is used 

to focus the energy, the spot area exposed to the focused laser beam becomes 

-32 -32 

0.49 x 10 in. . A 25,000 W beam focused onto an area of 0.49 x 10 in. 

6 2 

results in a density of 51 x 20 W/in. . A power density of this magnitude is 
sufficient to vaporize any known material. 

Not all of the material is removed by evaporation, however. As shown in 
Fig. 2-69, laser machining is basically a high-speed ablation process. The 
evaporation of a very small portion of liquid metal takes place so rapidly 
under the high intensities of a focused laser beam that a substantial impulse 
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Is transmitted to the liquid. Material leaves the surface not only through 
evaporation but, also. In the liquid state at a relatively high velocity. 


The amount of energy needed to raise a volume of material auch as stainless 
steel to its vaporization point can be calculated approximately as the energy 
required to raise the metal to its vaporization point plus the latent heats of 
fusion and vaporization, as shown. The energy required for vaporization of 
1.0 gram of metal requires: 


1. Heating from room temperature to melting point: 

E, - C(T - T ) ■ 0.11(1535 - 20) - 167 calories 
1 mo 

2. Changing from solid to liquid at T : 

m 

F "2 * l.j = 65 calories 

3. Heating from melting point to boiling point: 

E, * C(T. - T ) = 0.11(3000 - 1535) = 161 

3 b m 

4. Changing it from liquid to vapor at T^: 

E . = L = 163d calories 

4 v 

and Ej + E,, + E. ( + = 2023 calories = 8500 joules per 1 gram 


where : 


m 


specific heat in calories/gram 

ambient temperature in°C 

mel.ing temperature in°C 

boiling temperature in°C 

heat of fusion in calories/gram 

heat of vaporization in calories/gram 
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If the laser output energy is on the order of 35 to 40 joules: 


35 

8500 


0.0041 gram removed to 


40 

8500 


0.0047 gram removed 


This gives close agreement with the experimental results of between 4.1 to 
5.1 mg (Ref . 3) . 

However, under balancing conditions, the shaft will be rotating. This will 
reduce the time the focused beam impinges the shaft or wheel surface. Thus, 
removal rates for laseis of the size discussed will be less than that calculated 
for the static condition (Fig« 2-70). 


Father, the removal of material also is dependent on reaching the threshold 
energy level for vaporization. Since the laser power supply wattage is limited, 
reducing the pulse duration to accommodate a rotating shaft results in decreased 
metal removal rates. Also, as the depth of the beam penetration increases, the 
amount of material removed decreases. Generally, this limit of depth at 1 msec 
pulse duration is on the order of 3 to 4 mm, since the beam becomes unfocused at 
greater depths. 

Thus, to take advantage of the laser's capabilities, the balance zon' surface 
speed should be minimized. This required a procedure where data are acquired 
at or near rotor system critical speeds, followed by a reduction in shaft speed 
for the balancing process. In light of the HP0TP and HPFTP bearings running 
unlubricated or nearly so, this speed reduction would appear beneficial. 

The major drawback, however, is that even running at slow speeds, the bearings 
may be subjected to heat and radial loads for unacceptable periods of time. For 
example, if 5 grams of material were to be removed from one balance plane, using 
a laser pulse duration of 1 msec with a recharge cycle of 2 seconds and removal 
of 4 mg per shot would require 1250 shots over a 42-minutes period. This time, 
of course, can be reduced in a number of ways. First, multiple lasers could be 
controlled and fired simultaneously to increase the removal rate per shot. 
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Figure 2-70. Rate of Material Removal per Shot 



Secondly, the beam(s) could be focused at the maximum radial distance possible, 
thus maximizing the effect of the metal removed. Material removal rates can be 
boosted if the power of the laser is increased. For instance, in laboratory con- 
ditions, up to 2 oz/min can be removed by a 200-joule laser with a shot rate of 
8 shots/sec. Actual removal rates for balancing would be subject to a development 
program where the removal rate could be optimized while varying the energy release, 
pulse rate, rotational speed, and distance from lens to work piece. 

Another consideration relative to the SSME turbopumps is that the vaporized and 
liquefied metal that is removed during each shot is deposited subsequently on 
adjacent pump interior surfaces, thereby contaminating the assembly. Develop- 
ment of some means of collecting this debris is prerequisite to use of laser 
correction for in-housing balancing of the SSME turbopumps. 

Maintaining Cleanliness 

Turbopump interiors must be maintained at the high levels of cleanliness 
specified for propellant handling systems. In order for balancing operations 
to comply with this requirement, the correction sequence must entail the 
following steps: 

1. Protect the environment of the access port to prevent entry of con- 
taminants. In field conditions, a plastic barrier or tent may be 
erected around the turbopump and supplied with filtered air delivered 
by fan. 

2. Personnel must wear nonlinting caps, gloves, and clothing. 

3. The external surfaces of the pump casing adjacent to the access port 
must be flushed prior to opening with compatible solvent such as 
Freon Tf. 

4. Tools used to contact all interior purts 'f the pump casing must be 
flushed with similar solvent. 

5. Correction must be performed in accordance with a noncontaminating 
method. 

a. Determine angular orientation of correction site. 

b. Perform correr.rli. s by exchanging weights, or removing material. 

c. Flush and purge access port area. 
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6. Replace Access covers, securing closure and perform leak checks of 
closure seal. 

Post-Correction Cleaning 

In the event that contaimination of the turbopump interior has occurred, return- 
ing the assembly to certifiable propellant-clean condition requires disassembly, 
detail cleaning of components, and reassembly in a clean room meeting stringent 
cleanliness standards. Ultrasonic agitation has been tested as a means to dis- 
lodge particulate matter within assemblies; it was intended that the particles 
gravitate to some lower position, be collected, and removed. In practice, 
however, some particles always remain trapped. Bubble nucleation also has been 
investigated as a means of dislodging all particles from injector assemblies. 

In this method, the assembly to be cleaned is heated above the boiling point of 
the flushing medium, which may be a solvent such as Freon Tf. LN^ also has been 
used. As in the case of ultrasonic cleaning, particulate matter was not totally 
removed from assemblies. From a practical standpoint, these methods may remove 
most of the particulate contaminants, but the only way to guarantee their absence 
in an assembly is by detail cleaning the components prior to assembly, which 
must be executed in a controlled environment. Therefore, it is concluded that 
the only acceptable way to maintain propellant cleanliness of the turbopump 
passages is to ensure that they are not contaminated by the balancing procedure. 

In-Housing Balancing Speeds 

The Balance Simulation Studies showed that high speeds are necessary to obtain 
a high-quality balance when a limited number of balance correction planes are 
used. Evacuation of air from the turbopump environment will be necessary to 
avoid excessive heating by windage. How these speeds can be achieved in a 
vacuum without producing bearing damage is a critical issue that must be 
resolved if high-speed, in-housing balancing is to become a viable operation. 

A development program will be necessary to identify bearing operational limits 
and cooling/ lubricat ion techniques compatible with the vacuum enviionment and 
the propellant cleanliness requirements of operational turbopumps. Potential 
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c (sequences for bearings of high spin speeds in a vacuum include bearing 
overheating and surface degradation. 

At high spin speeds, heat is generated at a significant rate in the turbopump 
bearings, which are subjected to relatively high axial preloads. The potential 
for occurrence of negative internal clearance and thermal runaway increases at 
room temperature because inner race fits are generally tighter when compared to 
chilled conditions. Options available to avoid excessive bearing temperatures 
at high balance spin speeds are: 

1. Employ slave bearings (which defeats one of the purposes of in-housinR 
balance, i.e., eliminating relocation variation) 

2. Limit balance speed 

3. Limit balance run duration 

4. Provide fluid coolant (thus makit.g maintenance of vacuum more 
difficult) 

If balancing is performed in a vacuum environment, heat transfer away from the 
bearings mainly occurs by conduction. Without convective cooling, the bearing 
temperature rise will be rapid. Estimates of the permissible run time versus 
speed relationship for the high-pressure turbopump bearings are shown in 
Ti*. 2-71. Assumptions were: 

1. Maximum bearing temperature = 300 F 

2. Bearing initial temperature * 70 F 

3. No conduction of heat from bearing 

Surface degradation is currently controlled by providing the bearing races, 
balls, and cage with dry-film coatings. Determination of the durability of 
this treatment at higher speeds will require test and evaluation. Since 
dry-film coatings have limited wear life, development of some means of 
reinforcing or replenishing them may be necessary prior to implementing 
high speed, in-housing ualancing. 


RI/RD82-285 

2-98 






TASK III: BALANCE CONSISTENT ROTOR DESIGN 


Even the most ideal balance procedure will not produce satisfactory results if the 
initially balanced condition is not maintained. Examples of factors that can pro- 
duce post-balance response changes are: 

1. Rotor joints that permit relocation variations or shifting *1 rotor 
components 

2. Turbine blade coating, spalling, or uneven erosion of blade material 

3. Change in bearing stiffness due to wear 

4. Synchronous moment resulting from the effect of synchronous radial load 
on bearing tilt 

In Task III, these items were assessed for their effect on balance consistency. 
ROTOR JOINTS 

Elements that connect the individual rotor components into an assembly must provide 
location as well as transmit torque. 

To maintain an established balance condition, these rotor joints must be capable 
of minimizing: 

1. Mass eccentricity due to tolerance accumulation 

2. Mass mislocation at initial assembly due to clearances 

3. Relocation variations arising from operational strains or subsequent 
disassembly/reassembly cycles 

Achieving these goals can be approached by analyzing the capabilities of various 
rotor joints and selecting favorable configurations. 
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Tolerances 


Most turbopump rotors. Including the SSMF. HOPTP and FPFTP, incorporate several 
elements with appreciable mass, such as Impellers or turbine wheels. If all 
these elements can be assembled with perfect concentricity and normality of 
locating surfaces, the rotor could be balanced by performing a single-plane 
force balance of each component separately to correct its mass distribution. 

Good quality, low-speed balance machines can detect eccentricities of 25 p-in. 
or less. In the absence of shifts in the rotational center and bearing dead 
bands, a 1-pound rotating element could therefore be balanced to within 0.011 
gram-inch. However, from a practical standpoint, the eccentricities inherent in 
precision rolling element bearings, as well as those concentricities attainable 
with even the best machining techniques, are much larger, being in the range of 
100 p-in. Therefore, unless the same bearings are used for balancing and opera- 
tion, and unless angular indexing of components is maintained, centering ability 
tolerances may result in residual unbalance values of four or more times greater 
than the quoted balance machine capabilities. 

The radial eccentricities of pilots due to tolerances result in unbalance by dis- 
placing the components's mass center from the ideal geometric location. In addi- 
tion, nonnormality of axial mating faces of components in an axially stacked 
assembly, as shown in Fig. 3-1, produces a bowing of the rotor so that mass 
elements between the bearings -ire dJsplA-ce^. racially by a further increment, 
adding to the assembly's unbalance. Bending of the rotor in response to centt 
fugal force aggravates this effect, increasing unbalance response. 

Radial Splines 

Radial splines are one of the most efficient disassemblable torque transmission 
elements available. 
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Figuie 3-1. Effect of Non-Normal Axial Mating Surfaces 
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Side Fit Splines . The most popular spline configuration for turbomachinery 
application is the side-fit involute spline. When required, misalignment and 
axial motion can be accommodated by selection of appropriate tooth clearances. 

Even when these features are not required, the tooth clearances necessary for 
interchangeability detract from the spline's radial locating ability. The magni- 
tude of tooth clearances provided for manufacturing variations depends on the 
spline size. As an example, a 16-pitch, 1-inch-diameter, side-fit spline may have 
a diametral clearance ranging from 0 to 0.005 inch. Application of torque tends 
to center the mating splines; however, the center established is subject to tooth 
spacing tolerances and friction, so that centering accuracy even with torque 
applied of better than 0.0005 inch cannot be relied upon with confidence. 

0D or ID-Fit Splines . Splines with interference fit at tooth roots depend for 
location upon close control of major or minor diameters of the teeth. Location 
capability within 0.0005 to 0.0010 T1R can be attained by special machining pro- 
cedures. A disadvantage of these splines is that a flat root is required, pro- 
ducing a higher root bending stress. 

Tapered Root Splines . These splines are drawn together axially on a shallow taper, 
thus eliminating radial clearance. However, radial location is subject to machin- 
ing tolerances. They are rarely used for high-speed applications because present 
machining techniques cannot produce interchangeable parts with sufficiently small 
tolerance range to attain precise location while simultaneously avoiding over- 
stressing the inner diameter of the -otating element. 

Of the available types, involute side-fit splines can be machined most accurately. 
Tc ;ain radial location, they are generally used in conjunction with cylindrical 
radial pilots, thus separating the Joint functions into torque transmission and 
radial location. However, an. inaccurately centered spline member may produce 
radial loads and shaft bending if highly torque-loaded. 
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Radial Pilots 


Short cylindrical mating surfaces (pilots) between components of rotating machinery 
are the most common means of controlling alignment and concentricity. Piloted 
joints require an interference fit to ensure that radial position and, thus, 
balance condition, is maintained. Unbalance and shifted shaft critical speeds 
result if the joints loosen during operation. The designer must be informed of 
a..d consider all steady-state and transient conditions to ensure that pilot 
looseness never occurs. In general, operational transient conditions are not 
fully known during the design and initial development of a new turbopump. Dis- 
covery of the effects of transients will often creat ? the need for a redesign. 
Effects that must be evaluated to ensure maintenance of tight fits are: 

1. Elastic deflections due to centrifugal forces 

2. Thermal growths 

3. Radial growths due to Poisson's effect 

4. Growths from pressure loads 

5. Axial contraction due to Poisson's effect and centrifugal growth 

Centrifugal Growth . The radial growth of a solid disk, a disk with a central hole, 
and a thin disk or cylinder due to centrifugal force is defined by the following 
equations: 


Thin cylinder: 

6 


R(w) 


Solid disk: 


R(u) 


2 3 
pu)R 




2 3 

pu) R o 
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Disk with hole in center: 
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or. In dimensionlest form: 
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Simple approximation*: 


2-3 

x . P<*> R s . 

K(W) R ’ K 2 


R +R . 
o 1 


2 3 

Dividing each of the equations by to R q /E produces a normalized dimensionless 

parameter that is a function of the ratio of OD to ID. This relationship, plotted 

in Fig. 3-2, provides the centrifugal radial growth of the inner and outer radius 

of disks. This figure also shows that an approximation using the average radius 

and the thin cylinder equation is adequate for disks with an inside to outside 

radius ratio (R./R ) >0.25. 

i o 

Ther mal Growth . Thermally induced growth can be beneficial in pilot design if the 
mating materials are selected so that the relative growth tightens the joint. Care 
must be exercised to prevent yielding of the pilots through excessive thermal 
stress. Thermal radial growth, A, is expressed as a function of the temperature 
change, AT, and material thermal expansion coefficient: A = R “ AT. In evalu- 

ating the effects of thermal growth on the pilot, both steady-state and transient 
temperatures must be considered. 
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Figure 3-2. Centrifugal Growth of a Disk 
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A convenient reference (Table 3-1), which lists thermal expansion at various 
temperatures for several materials, is presented to aid in material selection. 

Using the information shown, materials and dimensions can be chosen to exploit 
differential thermal growth to maintain elastic preloads. For example, the pilots 
between the titanium impellers and the Inconel 718 sleeves in the HPFTP rotor are 
maintained under preload under all thermal and speed conditions of operation. In 
a similar manner, thermal properties are used to maintain adequate axial preload 
of impellers and sleeves by tension in the tie bolt. 

Preload Effects . Radial growth in axially loaded rotor components from the Poisson 
effect can lead to loss of pilot preload, permitting unbalance and changes in rotor 
critical speeds. These components experience a radial growth directly related to 
these loads: 



where 


A = radial growth 
A * stressed area 
P is the tensile load 
E = elastic modulus 
R = pilot radius 
U = Poisson's ratio 

This effect can cause radial growth of ball bearing inn^r races, producing radial 
clearance between the inner race and shaft, leading to an increase in residual 
unbalance . 

Errors in parallelism and normality of component axial faces subjected to axial 
preloads can introduce moment unbalances in addition to those arising from bowing 
of the rotor. Although the initial unbalance thus created could be corrected at 
low speeds, centrifugal loading in operation will create radial deflection and 
additional unbalance. 
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TABLE 3-1. MATERIAL THERMAL CONTRACTIONS 


a AT, MILS/INCH, 70 F TO TEMP 


MATERIAL 

-400 F 

-200 F 

10C0 F 

2000 F 

INCONEL 600 

-2.2 

m 

7.6 

mm 

INCONEL 625 

-2.3 


7.4 

KB 

INCONEL 718 

-2.4 

DU 

7.6 

KS9 

WASPALLOY 

-2.5 

-1.6 

7.2 


RENE 41 

-2.2 

-1.6 

7.7 

19.2 

HASTELLOY B 

-1 .8 

-1.4 

6.0 

15.7 

HASTELLOY C 

-1.4 

-1.0 

6.7 

16.0 

HASTELLOY X 

-Z. 1 ' 

-1.8 

7.6 

17.5 

MONEL 400 

-3.0 

-1.7 

8.5 

19.9 

MCNEL K500 

-2.5 

ESI 

7.8 

— 

ASTROLOY 

-2.6 

E9 

8.0 

23.5 

1018 

-1.8 

-1 .4 

7.5 

15.0 

4130 

-2.0 

-1.4 

7.2 

15.0 

4340 

-2.0 

-1.5 

7.2 

13.8 

302 STAINLESS STEEL 

-- 

— 

8.9 

20.0 

304L STAINLESS STEEL 

-3.2 

-2.4 

7.0 

20.2 

310 STAINLESS STEEL 

-3.2 

-2.4 

8.6 

20.0 

315 STAINLESS STEEL 

-3.0 

-2.3 

8.7 

19.2 

321 STAINLESS STEEL 

-3.2 

-2.4 

9.0 

21.0 

347 STAINLESS STEEL 

3.0 

-2.4 

9.0 

20.4 

410 STAINLESS STEEL 

-1 .8 

-1.3 

— 

-- 

416 STAINLESS STEEL 

-1.9 

-1.5 

6.0 

13.5 

17.4 PH 

— 

— 

6.0 

14.2 

A286 

-2.8 

-1.9 

8.8 

21.0 

INCONEL 903 

-1.6 

-1 .0 

4.5 

15.0 

440C 

-2.0 

-1 .4 

6.7 

15.4 
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Extended Pilot s 

In a configuration with an extended pilot, radial tightness is maintained at one 
axial station by utilizing the elastic properties of the piloted component, and 
radial growth at another station. The principle is illustrated in Fig. 3-3 in 
nondimens ional form; A, the radial deflection of the free end of the ring, can 
become negative, thus tightening the pilot. The action of the ring can be visu- 
alized by considering the ring as a beam on an elastic foundation with the left 
end subjected to a fixed radial displacement. The extended pilot principle is 
used in practice in the HPOTP main impeller (Fig. 3-4). 

To illustrate the benefit of an extended pilot, the results of a computer model 
of a drive turbine wheel were reviewed. The finite element computer model of the 
disk with its extended pilot is shown in Fig. 3-5. The deformation of the inside 
diameter due to centrifugal forces is indicated on this figure. The deflection 
of the pilot was compared to that predicted for the radial deflection occurring 
at one end of a cylinder when the opposite end is displaced radially. To make 
the comparison valid, the radial growth due to rotational speed must be added to 
those arising from imposed deflections at the opposite end. The finite element 
model was run with two additional extended pilot lengths to provide more data. 

The results of this comparison are as follows: 


Si 

U 

PREDICTED 

A/Ao 

(FIG. 3-3) 

SPEED 

CORRECTED 

A/A 0 

PREDICTED BY 
FINITE ELEMENT 
MODEL A/Ao 

0.50 

2.39 

-0.260 

-0.181 

-0.120 

0.75 

3.59 

-0.100 

-0.020 

0.014 

1.00 

4.79 

0.000 

0.080 

0.082 
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Figure 3-5. Extended Pilot Finite Element Model 
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The differences between th.'i A/A q predicted by the curve and by the finite element 
model are a result of toroidal rotation that exists at the "fixed" end of the 
cylindrical portion of the disk model. The finite element model confirms that 
by properly selecting the geometry of the cylindrical extension, a pilot can be 
obtained to tighten radially with speed. 


Limits on Interference Fits . Given the requirement that joints be designed so that 
preload remains constant or Increases under operating conditions, it is essential 
that the consequent stresses do not cause yielding or creep, with resultant pre- 
load loss. Interference fits must be selected so that resulting hoop stresses do 
not exceed material yield stress divided by some safety factor. An approximate 
expression for limiting the interference fit would be: 


Interference < 


(r)(f ty )* 

( F. ) ( F . S . ) 


where 

R “ radius at interference fit 
E * elastic modulus 
F ty * tensile vield stress 
F.S. = factor of safety 


The lower bound for interference fit is that value that will maintain joint tight- 
ness at all times. Machining tolerances of the fitted parts may make it difficult 
to select an interference range that can provide tightness throughout the load 
cycle for its minimum value and yet not exceed the allowable stress for its maxi- 
mum value. At present, practical machining techniques can achieve tolerances of 

+0.0000 , t r\ i +0.0002 , * 

-0 0003 ° n the ID and 0 0000 ° n tae °* mat ^ n 8 surfaces. The stress produced 

is size dependent; for example a 0 to 0.0005- inch interference fit on a 10- inch 
diameter would cause hoop stresses from 0 to 1500 psi. However, at a 1.0- inch 
diameter, this same fit range would cause stresses from 0 to 15,000 psi. 


*In conservative practice, it is assumed that the total interference must be 
accommodated entirely by one of the mating members. 
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This effect is especially serious for small turbopumps, where avoidance of pilot 
yielding and stress corrosion cracking of bearing inner races will require pre- 
cision tolerances proportional to size. 


- . ., + 0.0000 _ . . 

This, if a q 0003 to ^ erance can ° e 


achieved on a 3-inch inner diameter, an equivalent stress range from a 1-inch boiti 
fit will require a tolerance of + o*qo(H' 


An additional factor to be considered in interference fits is the effect of 
tolerances in thickness of protective coatings or platings. In the case of soft 
metal plating or fragile coatings such as MoS^, assembly techniques should include 
shrinking of interference fits by heating or chilling of parts, as opposed to 
axial pressing. Undesirable consequences of shearing away of the coating include 
reduction of interference fit, loss of radial position precision, and possible 
axial runout caused by the presence of sheared-off material lodging between axial 
mating faces. More precise fits can be obtained between plated surfaces if the 
after-plate dimensions are the result of machining or grinding. Otherwise, the 
plating thickness tolerances apply twice for each plated diameter and are additive 
to the me Mning tolerance specified for the underlying surface. 


Centering Bolts 

Rotating ass mblies consisting of a number of axially clamped parts (a ".'itack"), 
preloaded by through bolts should incorporate additional radial positioning 
features such as radial pilots or curvic couplings to aid in balance repeatability. 
"Body bound" piloted bolts are not the optimum choice for balance-consistent design, 
because an angular error can arise in the bolt threads. This error can be minimized 
by incorporating line-reamed matched holes in conjunction with press-fitted axial 
pins, as shown in Fig. 3-6. 
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Figure 3-6. Bolt Body Location 

Curvic Couplings 

T.i“ Fixed Curvic Coupling* is a curved-tooth face spline that incorporates pre- 
cision radial and axial location with high torque capacity. These couplings are 
often used, as in the SSME, to attach turbine wheels to mating rotor elements. 
Maintenance of proper axial clamping load is necessary to realize their radial 
location and bending stiffness potential. Therefore, the effects of centrifugal 
growth and thermal differential growth m,st be controlled in rotor designs. 

Proper preloading procedures must be incorporated Into assembly, reassembly, or 
balancing operations that require fastener loosening or removal in curvic coup- 
lings. The preload must be applied in a manner that precludes excessive differences 
in the circumferential distribution of clamping load. In a multifastener Joint, 
application sequence must be employed. A nonuniformly distributed preload can 
lead to circumferentially varying elastic deformations of the curvic teeth and 
resulting joint angularity. The design book published by Gleason (Ref. 5) recom- 
mends that the curvic' s bolts be stretched equally within 1% on boltup. 

*The Curv*.* Coupling was developed by the Gleason Works, Rochester, New York, in 
1942. 
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The joint, once pre loaded, should never have more than one fastener removed and 
installed at any time. Removal of several fasteners may Induce misalignment of 
the joint and alter the balance. These forms of control are exercised in the 
assembly of the Space Shuttle high-pressure turbopumps. The final preload of the 
turbine wheel joint in the HPOTP is attained by controlling the stretch of the 
12 turbine joint bolts. During assembly and balancing, an approximately uniform 
preload is applied by torquing each bolt to 210 lb-in. Subsequent in-housing 
balancing is accomplished by changing bolts of differing weights, removing and 
then retorquing only one bolt at a time. A potential improvement in balance 
repeatability might be achieved if the bolt preloads were always controlled by 
stretch rather than torque. This would reduce the variation in load caused by 
the bolt thread friction. 

Maintaining angular orientation of all rotor components is necessary to achieve 
balance consistency; it is necessary to match-mark all rotating components with 
respect to a reference. On subsequent assemblies, all parts must be reinstalled 
in the same clocked location as indicated by the reference marks, including all 
major and minor components, all bolts, nuts, washers, and a^y other attached 
parts. If the rotor is disassembled more than once, care should be takc-vi to 
realign marks. Fasteners or rotating flanges such as bolts, nuts, washers, and 
locking devices should be positioned to the outermost allowable radial location 
or be piloted radially, thus eliminating radial shifting of these parts due to 
centrifugal forces during rotor spin. 

ROTOR DESIGN ALTERNATIVES 

The balance consistency potential of a given rotor can be increased by reducing the 
number of joints or by improving the location capability of the joints. The follow- 
ing concepts show how these effects might be achieved for the SSME high pressure 
turbopumps. 

The HPFTP rotor incorporates sleeves initially selected to permit the use of 
toothed-labyrinth interstage seals requiring a material other than the titanium 
used for the impellers. Since these seals have subsequently been changed to a 
straight smooth configuration, incorporation of the seal surfaces into the impellers 
is considered feasible; based on this concept, two alternative arrangements of the 
HPFTP rotor were proposed for consideration. Both arrangements reduce the number 
of axial interfaces, thus improving the stack alignment capability. In the present 
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HPFTP rotor configuration, shown in Fig. 3-7, alignment of the pump rotors depends 
on ...e parallelism of three separate seal spacers and the mating impeller faces. In 
Alternate 1, shown in Fig. 3-8, the number of axial locating faces is reduced by 
extending the second- and third-stage impeller hubs, thus replacing separate spacers 
at the interstage seal and liftotf seal locations. Another method of eliminating 
separate spacers is shown as Alternate 2 (Fig. 3-9). Either of these alternative 
methods is capable of reducing composite error and runouts of the rotor. The 
consequences of implementing both alternatives are similar: 

Favorable : 

1. The potential for producing a bowed rotor is reduced when the number of 
axial interfaces is reduced. 

2. Rotor strbility may be improved by reduction of number of locations with 
potential hysterisis. 

Potentially Unfavorable : 

The elimination of Inconel 718 sleeves will require tighter assembly fits at 
radial pilots in which both members are titanium because thermal contraction 
differences between the titanium impeller and the Inconel 718 sleeve no longer 
exist. 

The HPOTP rotor configuration offers less opportunity for reduction in the number 
of joints used. However, potential improvement in piloting is co isidered achiev- 
able for the preburner impeller, and for the shims used to adjust the axial posi- 
tion of the main impeller and preburner impellers. As shown in Fig. 3-10, the 
preburner impeller is radially positioned by two interference fit pilots, (3) and 
both inboard of the drive spline, ( 7 ). It appears feasible to provide a 
pilot outboard of a shortened drive spline, as illustrated in Fig. 3-11. While 
preliminary evaluation shows adequate spline capacity, a complete analysis of this 
joint's torque loads and capacity is required prior to adopting this design modi- 
fication. An alternate method of achieving an outboard pilot is shown in 
Fig. 3-12, where the original spline length is retained and a secondary piloting 
surface is incorporated onto a threaded insert. While the preburner impeller is 
reworkable to this configuration, structural analysis of the impeller must be 
performed to determine if the cross section at (a) is adequate. 
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Piloting of shims can eliminate some sources of unbalance; if the shims atQ)and 
(5) in Pig. 3-10 are thin enough to move radially into the bore bottom reliefs 
at some unbalance will be produced. (The maximum estimated is 0.5 g-inch.) 

Single, piloted shims, as shown in Fig. 3-11 and 3-12, will have minimal radial 
shift. In order to use single shims only, they must be custom-machined or a 
large number of sizes must be stocked. From a practical standpoint, achieving 
the desired shim total thickness with fewer stock shim thicknesses can be obtained 
by combining shims. If this aspect is an important consideration, multiple shims 
can be piloted by using the alternate method of providing impeller and multiple- 
shim piloting with a ring insert as shown in Fig. 3-13. The pilot ring bore is 
finish machined after installation in the impeller, and becomes a part of the 
impeller. Although the piloting obtained is not as positive and accurate as 
that for the concept of Fig. 3-10, the insert pilot has the advantage of rework- 
ability from the present preburner impeller. 

Nonlndexable components, such as the large-diameter axial stack clamping nuts ( (^ 
and (§) in Fig. 3-10) and their locking devices are a source of residual unbalance. 
A long-term goal for achieving maximum balance consistency would include efforts 
to replace the individual nuts by a tension bolt, such as is used in the HPFTP 
rotor. 

TURBINE BLADE EROSION 


Changes in balance condition arise when the individual turbine blades are eroded 
at varying rates by the high-velocity drive gas. In the HPFTP, the first turbine 
stage blades are subjected to the highest gas temperature and, therefore, are the 
most subject to erosion. The second-stage blades do not experience significant 
erosion. Weight change can take place in two ways: the 'coating and substrate 

blade material can be gradually eroded as service time increases; a more discon- 
tinuous process is spalling of flakes of blade coating. 
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Erosion 


Turbine blade erosion is a graduated process that causes a relatively uniform 
weight loss among bludes and results in minor unbalance change. Weighing of 
representative blades before and after hoi-fire testing yielded the following 
estimate of the consequences of erosion. The average weight change was 0.0866 
gram, with individual blade weight changes ranging from 0.066 to 0.126 gram. 

Spread over the blade complement, this distribution results in a balance change 
of 0.13 gram-inch in 1000 seconds of operation. 

Spalling . Since elimination of a zirconium outer covering over the Nicrally blade 
coating, spalling has not been encountered. When the outer coating was used, 
large unbalance changes were sometimes observed. These changes were identified by 
an increase in turbopump vibration in a run-to-run progression. Vibration Increases 
of this type have not been observed since removal of t le outer coating. 

BEARING TILT 

To assess the effect of bearing tilt on rotor response after the balancing process 
has been completed, a finite element model of the HPOTP was established. Critical 
speeds were predicted, and rotor unbalance response calculated using the worst- 
case unbalance distribution. Correction weights were then simulated and a rotating 
moment applied. The applied moment was representative of the moment generated in 
the bearing response to synchronous radial loads from the residual unbalance. 

The single-level finite element model is shown in Fig. 3-14. Bearing stiffnesses 
used were 500,000 lb/in. for each bearing and 30,000,000 lb/in. for the hearing 
pedestals. Critical speed mode shapes are shown in Fig. 3-15. Predicted undamped 
critical speeds using this model are shown in Table 3-2 and are compared to the 
Rocketdyne calculated single-level undamped critical speeds This comparison is 
provided to show that the model generated is representative of the HPOTP modes 
and, thus, may be useful in assessing the effect of bearing tilt on rotor residual 
unbalance response. 
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Figure 3-15. SSME HPOTP Mode Shapes 






-n* 


ORIGINAL * 

OF POOR Q U Al ' 


TABLE 3-2. HPOTP CRITICAL SPEEDS 



FIRST CRITICAL 

SECOND CRITICAL 

THIRD CRITICAL 


SPEED, 

SPEED, 

SPEED, 


RPM 

RPM 

RPM 

ROCKETDYNE 

1 1 ,028 

34,379 

46,903 

MTI 

11,660 

34,328 

48,970 


A comparison of Fig. 3-16 and 2-62b shows fairly close agreement between the post 
balance 30,000 rpm response calculated from Rocketdyne Influence coefficient data 
and the MTI model. Differences between the two results can be attributed to the 
simplified model used for the bearing tilt study wherein all bearings are symmetric 
and a constant 500,000 lb/in. stiffness is used regardless of speed, possible dif- 
ferences in shaft flexibilities, and wheel inertia values used. 

The approach used to determine the effect of bearing tilt on the HPOTP rotor 
response was to: 

• Determine bearing loads from rotor dynamic response analysis for 
balanced conditions 

• Perform quasi-static bearing analysis of 57 mm bearing for various 
radial unbalance loads 

• Perform response analysis adding bearing tilt moment to the bear- 
ings in combination with preceding unbalance forces 

Table 3-3 shows bearing loads predicted for balanced and unbalanced states of the 
modeled HPOTP rotor. Using the radial bearing load for the No. 3 bearing as a 
basis and assuming an 1100 pound axial preload and a shaft rotational speed of 
30,000 rpm, the bearing resultant moment was calculated and plotted as shown in 
Fig. 3-17, where it is shown that for a 136-pound radial load, a moment of 76 
in. -lb is generated within the bearing, tending to cause the bearing to tilt. 
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Figure 3-17. HPOTP - 57mm Bearing Moment vs Applied Radial Load 



ORIGINAL PAGE 13 
OF POOR QUALITY 


TABLE 3-3. HPOTP BEARING LOADS AT 30,000 RPM 


BEARING 

NO. 

AXIAL 

LOCATION, 

INCHES 

RADIAL LOAD 

, POUNDS 

UNBALANCED 

BALANCED 

1 

3.69 

196.8 

54.8 

2 

4.44 

260.4 

28.4 

3 

14.32 

411.2 

136.3 

4 

15.44 

381 .8 

87.8 


This moment was then applied at the bearing synchronously in conjunction with the 
unbalance forces that generated the 136-pound bearing load to determine the rotor 
system response. This moment was applied as a synchronous rotating moment since 
it is generated by a synchronously rotating force. As seen in Fig. 3-16, the 
effect of the applied moment is minimal for the balance condition studied. Real- 
izing that unbalance forces do not constitute the total load experienced by the 
bearings in generating a moment, a second case assuming a 400-pound bearing load 
was completed. As seen in Fig. 3-18, even with the resulting 220 in-lb moment, 
rotor response is relative insensitive. 

BEARING WEAR 

Severe operating conditions and minimal lubrication will result in wear of bear- 
ing balls, races, and cage. Although evenly distributed wear will not increase 
residual unbalance, some change in rotor dynamic response is to be expected due 
to a loss of bearing radial stiffness. Ball wear will affect stiffness in two 
ways: the effective internal clearance increases, while axial preload is dimin- 
ished. An estimate of the combined effect of ball wear in the HPOTP turbine 
bearing is illustrated in Fig. 3-19. The bearing stiffness falls from a range of 
(0.428 to 0.510) x 10*’ lb/in. with no wear and 1000 pounds preload to (0.360 to 
0.430) x 10* 1 lb/in. with ball wear of 0.0005 inch. The reduced ball size permits 
an axial relaxation of the bearing of approximately 0.003 inch, reducing the pre- 
load to 840 pounds. If the individual balls stay round, but wear to differing 
diameters, nonsynchronous shaft forcing functions can be introduced at multiples 
of cage speed of approximately 43% of the rotor speed. Balls worn out of round 
will impose essentially random radial forces to the rotor. 
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Figure 3-19. Bearing Stiffness vs Ball Wear, 

HPOTP Turbine End Bearing 

In summary, bearing wear can produce rotordynamic effects, but is not a probable 
cause of unbalance change per se. 

Turbine Blade Loosene s s 

Turbine blades with fir tree roots tend to shift position slightly, thereby alter- 
ing balance. This effect is responsible for an unbalance indication of 2 to 3 

gram-inches that varies from run to run. Three means of dealing with the balance «j 

variation include: 

1. Base the correction on the average unbalance indication from a number 
of balance runs. The runs are sequential, with no physical changes 
being made to the rotor between the stop and start cycles. This 
method is optional on SSME turbopumps when varying balance is 
encountered . 
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2. Placing a tension band around blades to remove looseness. This method 
is not favored because the blades are moved radially inward in the 
attachment clearance and, therefore, are not in their running position 
which, under the influence of centrifugal force, is toward the outer 
limit of the clearances. 

3. "Brooming," or pressing a nonmetallic paddle against the outer surface 
of the blades as the wheel rotates to move all blades to the same end 
of their axial and circumferential clearances. This method has been 
used with some success in balancing the SSME turbopump rotors. 

Considerations of the effects of turbine blade looseness include: 

1. Centrifugal force on the blades will force the blades solidly to the 
outer limit of their clearance. 

2. Power loading by the drive gas will put a moment on the attachment 
that is impossible to simulate in a balance operation. 

These aspects of operation make high-speed balancing and operational balance data 
acquisition attractive in performing accurate balancing. 
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APPENDIX A 


DYNAMIC MODELS 


LINEAR MODEL 


Consider the rotor casing system represented by Fig. A-l; equations describing 
the system in physical coordinates are: 




Fj sin (ij/t -f a^) 

i < 

F 0 cos (ii»t + a 2 ) 

0 

X 

0 


(A-l) 


where the terms are defined as follows: 


m R = rotor mass matrix 

m ( , = casing mass matrix 

K r = rotor free - free stiffness matrix 

= casing stiffness matrix — casing may be tied to ground 
G = rotor gyroscopic matrix 

Kj = partitions of stiffness matrix defining coupling elements 

such as bearings elastic properties* 

Cj C 4 = partitions of damping natrix defining coupling elements 

such as bearings damping properties* 

Fj and F 9 = rotor unbalance matrices (a^ and a 2 = phase angles locating 
unbalances in shaft axes) 
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Introducing the following transformations. 


N| ‘ M M K! ■ MW 

1 u r) and |u c ) ■ Rotor and casing modal coordinates 

Rotor modal matrix in physical coordinates 
Casing modal matrix in physical coordinates 



where 

Rotor free-free natural frequencies 
Casing natural frequencies 




(A-2) 


(A-3) 


Note that the modal data could be determined experimentally if desired for an 
existing rotor casing. 


The generalized coordinate q r (t) and q ^ ( t ) is expressed in terms of 

|q r (t)|= [A R ] |U R (t)| so that the coordinate | q ^ ( t ) ^ is a linear combination of 

|u (t)|. Matrix [A r ] can be determined so that all the matrixes such as the 

mass matrixfm,, o "Ibecome diagonal; this decouples the system equation and it 
L 0 R 

becomes a group of independent linear equations which are easier to solve. 
Coupling is not an intrinsic property of the physical system, but depends on 
the type of generalized coordinates used. 


RI/RD82-285 

A-3 


ORIGINAL PAGE IS 
OF POOR QUALITY 



This simplifies to the following by Eq. A-3: 
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^ A F 




A l K l \ 

A R K 2 A C 


A C K 3 'St 
A l K 4 A C 


\ C 2 


A C C 3 \ 
A l C 4 A C 


- A r (F 1 cos a 1 - F 2 sin a 2 > 
T 

■ A r (Fj sin -f F 2 cos a 2 > 


Now assuming synchronous motion with 

\/j » fa) 


d iU)t 

* R e 

V 


“ /rr ) 


Substitute from Eq. A-6 into Eq. A-5 along with appropriate derivatives of 
Eq. A-7 and retain only the real part. 


2 

“r 

l-». 1 

[V 

l r 

0 

Ljd 

k 

1 L 

K l. 

5] 

\\\ 

i 

K 3 

K J 


[ 


C 1 C ? R 

C 3 C 4 J °C 


F 2 - F 1 1 


(note that e has been cancelled from both sides) 
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(A-9) 


Solution of Eq. A-9 yields the modal responses that can be transformed to physical 
responses by Eq. A-2. 


This linear approach treats the fundamental problem and forms the basis for com- 
puting the rotor unbalance from measured responses under conditions in which the 
deadband is neutralized or modified in such a way that small rotor motions will be 
linear in character. 

NONLINEAR MODEL WITH DEAD BAND 

Consider a rigid rotor on bearings with dead bands, as shown in Fig. A-2: 



Figure A-2. Rotor With Deadbands 
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Now, for at eady-st at »• synchronous conditions, there arc two possible equilibrium 
states, as shown In Fig. A-3. 


j CASE 1 


f CASE 2 


z 



m (R>r)u) 2 


Z 



rr. <R.|>u>2 
K |RA) 


R-DISTANCl FROM ORIGIN TO S 


Figure A- I. Steady-State Equilibrium Positions 
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Assuming R > A 


• Case i 
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m(R + r)o) 


K(R-A) 


-mruj - KA 

2 

mw - K 


• Case 2 


m(R-r)u) 


K(R-A) 


_ mru> - KA 

R - r* — 

row - K 

This simple case indicates that there are two possible values of R. Accordingly, 
two values of unbalance eccentricity r could be calculated. For a very simple 
rotor this is not a serious problem, but for a flexible rotor the ambiguous results 
greatly complicate the solution. 

The general flexible rotor response problem with deadband in the bearings can be 
treated by simulation method as follows. Starting with the same system used for the 
linear models, the equation of motion in physical coordinates, including clearance. 


K 1 0 1 M Mol [ q R 1 . r MIK 
L° KJ|q"| " "MvIKJ lMK 


caJ q c 


Fj sin + Oj) 


F 2 cos (ii» t + a 2 ) F ( 

+ _ + - 


(A-10) 
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whe r«* 


m 


K’ 


m 


C 


CA 


F l— F 2 


C 


rotor and cuslng mass matrices 
rotor free-free stiffness matrix 
casing stiffness matrix 
rotor damp i hr matrix 
casing damping matrix 
rotor unbalance matrices 

forces on rotor due to connecting elements 
forces on casing due to connecting elements 
rotor gyroscopic matrix 


The C and C_. are for linear damping effects Internal to the rotor and casing 

K LA 

respectively. F and F are force vectors which are due to connecting effects 
K l 

between rotor and casing which Include the following: 

1. Linear damping additions 

2 . Linear stiffness additions 

1. Nonlinear hearing forces (deadband) 

4. Kuhhing contact 

5. Fluid film hearings or squeeze film dampers 
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Recall the transformat ions of Eq . A-2 and A-3; transform the equation of motion 
to modal coordinates in a manner similar to that used for Eq . A- 1 ): 



Define the following using the previous modal transf ormations: 


r.i 

T r. i T\ 1 - 

r 2 ] 

N 

L K J M 

L "J 


I 1 N M ■ 


r i 

l T r,i Ta i - 

rd 

l>J 

1 L J l «J 

L J 

[s 

r h m • 

H 

M 

' MM • 

1 — T* 1 
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[ A rJ | F 1 008 “l ' F 2 008 a 2 } " | F l) 

f R ] T | F 1 8ln n l F 2 008 a 2| " I F 2 } 

H m " n 

[v]' l>,l • l>.l 

Using the above relations, Eq. 3-11 can be rewritten as follows: 



(A-12) 


This relation is the modal differential equation of motion. Generally, digital 
or analog simulation rethods are used to integrate it it yield u^ and u^, which 
are transformed back into physical coordinates. Figure A-4 shows a typical 
flow chart for solving Eq. A-12. 
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Figure A-4. Simulation Flow Chart 


While this simulation approach permits calculation of a system's response to 
unbalance, it does not guarantee that all the solutions will be found. If measure- 
ment of unbalance is desired in a system with dead bands, an iteration scheme 
could be used to match experimentally measured rotor motions with analytical pre- 
dictions for a number of unbalance distributions to yield the actual unbalance. 
However, ambiguity would still be present so that additional experimenntal trials 
would be needed to verify the unbalance prediction. 
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APPENDIX B 


BEARING LOAD CAPACITY 

The calculated bearing fatigue lives as a function of speed and loading are shown 
in Fig. B-l (HPFTP), B-2 (HPOTP preburner end bearings), and B-3 (HPOTP turbine 
end bearings). These plots were used as the basis for establishing rotor balance 
limits. 
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RADIAL LOAD. POUNDS 


Figure B-3 
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APPENDIX C 


ACCURACY 

Attainable balance accuracy depends on: 

1. Instrument scale factor 

2. Noise level 

3. Angle resolution 

4. Correction increment 

Weight 

Angle 


Choice of system and usage will determine the magnitudes of these variables; an 
estimate of achievable balance accuracy for the high-pressure turbopumps is 
listed in Table C-l. 
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TABLE C-l. BALANCE ACCURACY ESTIMATE 

RESOLUTION 
50 u-IN. 

15 
0.05 

15 

FOR 20 mV NOISE LEVEL: 

DETECTABLE UNBALANCE = (0.0001 + 0.00005) IN. X 25 LB X454 G/LB 

= 1.7 G-IN. 

ATTA, ScT ANCE * (oIMtot) ’ ‘- 97G - IN - 

FOR 4 mV NOISE LEVEL: 

DETECTABLE UNBALANCE ’« (0.00002 + 0.00005) IN. X 25 LB X454 G/LB 

= 0.8 G-IN. 

ATTAINABLE BALANCE = r n'Q7V/n~Q7 \ = 1-00 G-IN. 

ACCURACY 


ACCURACY 


0.00005 IN. 


0.02 V 
200 V/IN. 


0.0001 IN. 


0.004 V 
200 V/IN. 


0.00002 


COS 15 


(ASSUME 3- INCH RADIUS) 
0.05 X 3 = 0.15 


ITEM 

SCALE FACTOR 200 mV/MIL 
NOISE LEVEL 
20 mV 

4 mV 

ANGLE RESOLTUION, DEGREES 
CORRECTION INCREMENT, GRAM 

ANGLE AVAILABILITY, DEGREES 
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APPENDIX D 


CURRENT SSME BALANCE PRACTICE 


The SSME high-pressure turbopump rotor components are balanced at 1400 to 1500 rpm 
using rigid rotor techniques. Since the rotors consist of planes, they are 
balanced sequentially in two planes as the components are added. Subsequent to 
the low-speed balance, the rotors are disassembled and reassembled into the 
turbopump casings. 

During turbopump assembly, the angular orientation of all rotating elements is 
preserved to minimize relocation variation. High points of bearing eccentricity 
are aligned to the same reference for both the slave balance bearings and the 
actual turbopump bearings. 

The assemblies are then given ar. in-housing trim balance. The HPFTP is corrected 
in the two outboard planes: at the inlet end bearing slinger and at the first- 

stage turbine wheel. The HPOTP is corrected only at the turbine end. 

HPFTP 

The HPFTP (Fig. D-l) component balancing procedure incorporates the following 
steps: 

1. The trial assembly is performed (Fig. D-2) to verify the TIR* of the 
stacked components. The various parts are reoriented, if necessary, 
to achieve low runouts. Tic ,arts are indexed and marked so the 
angular orientation is pr» wed in subsequent assemblies. 

2. Temporary corrections are made to the outboard components consisting 
of the first-st e impeller and the unbladed turbine wheels (Fig. D-3). 

A tooling sleeve is substituted in the part stack for the second- and 
third-stage impellers and their associated sleeves. 


*Total indicated runout 
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MATERIAL REMOVED FROM SlINGER 
FOR IN-HOUSING TRIM BALANCE 
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3. Permanent corrections are applied in a two-plane balance on the second- 
and third-stage impellers. For this operation, the temporary correc- 
tions of step 2 remain with the balance assembly returned t e the 
configuration of Fig. D-2 . 

4. Permanent two-plane balance corrections are made to the first-stage 
impeller and the turbine plane. The turbine blades are individually 
weighed and then are arranged to correct the unbalance of the turbine 
plane determined in step 2. Final component two-plane balance is 
achieved in the assembly shown in Fig. D-4. 

HPOTP 

The rotating components and rotor assembly of the HPOTP (Fig. D-S) are balanced 
in the following sequence; 

1. The preburner impeller, shaft, and second-stage turbine disk with 
no blades as shown in Fig. D-6, are component balanced on a Schenck 
hard bearing balance machine. The initial unbalance is determined 
for planes A and D* plane A is balanced by removing material from the 
preburner and plane D is balanced with temporarily attached weights. 

2. The first-stage turbine disk is added to the assembly as shown in 
Fig. D-7 with no blades. A two-plane balance is conducteJ on planes 
A and E. Plane E is balanced with a temporary correction weight; no 
change is made to plane A. 

3. Main impeller (Fig. D-8) is force and moment balanced as a component 
on a tooling arbor by removal of material. 

4. The assembly of Fig. D-7 is dismantled and reassembled with the main 
impeller installed on the shaft. Alignment marks are added to all the 
parts to ensure they are reassembled in the same orientation as when 
balanced in steps (1) and (2). At the same time, first- and second- 
stage turbine blades are positioned bv weight with a computer program 
to reduce the unbalances of the turbine disks, which were measured 

in steps (1) and (2) . 
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Figure D-6. HPOTP Turbine/Preburner Balance Assembly 
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Figure D-7. HPOTP Turbine/Preburner Balance Assembly 
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Figure D-8. HPOTP Main Impeller 
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5. The entire assembly of Fig. D-9 is then balanced in two planes, B and C. 
Plane B is balanced by material removal. Plane C is balanced by 
exchanging special turbine bolts chosen by weight from a selection of 
bolts within a specific weight range. The rotor must then be dis- 
assembled prior to reassembly in the turbopump casing. 

6. The pump assembly is then in-housing balanced by adjusting bolt weights 
in plane C (see In-house Balance section for further discussion). No 
correction is made to the preburner impeller. 


IN-HOUSING BALANCE 


In-housing trim balance operations are presently performed on all HPFTP and HPOTP 
turbopump assemblies, using ON^ gas to spin the rotors to 5000 rpm. A Hofmann 
UGA 2000 vibration analyzer is used to detect unbalance; angular reference is 
obtained by a photocell signal keyed to a mark on the rotor. The vibration analy- 
zer, control equipment and instrumentation for in-housing balance is shown in 
Fig. D-10. 

HPFTP rotors are in-housing trim balanced on slave bearings because the outboard 
position of the bearings permits bearing exchange without disturbing the rotor 
assembly joints. Thus, the turbopump starts its service life with new bearings 
without lubricant residue. Angular indexing of balance and final assembly bear- 
ings is maintained to minimize the relocation variation potential of the bearing 
exchange. Slave bearings are identical to the bearings used in service, but are 
prelubricated with WD40 and cooled with GN,, during balance operation. The HPFTP 
setup for in-housing balance is shown in Fig. D-ll. Trim corrections are applied 
to the turbine plane by addition of washers under the turbine bolt nuts. Per- 
manent correction by material removal is performed on the pump inlet slinger. 

Trim balance corrections of 1/2 to 4 grams have been made on the pump inlet 
slingers, while corrections of 1/2 to 6 grams have been made in the turbine plane. 
Five to eight runs of approximately 1-minute duration are normally required to 
define the corrections needed in both planes. The vibration data are processed by 
a program stored in a calculator. The same data are fed into the Univac computer 
system, which provides permanent printed records of the unbalance conditions and 
corrections applied. 
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HPOTP rotors are trim balanced by exchanging turbine bolts of varying length 
(weight) in turbine plane, which is effectively located midway between turbine 
disks. Dry-film lubricant (MoS_) has recently been added as a treatment for the 
turbopump bearings. Assemblies with these bearings are given no further lubrica- 
tion during in-housing balancing. A total balancing rotation time limit of 
2400 seconds is observed for in-housing trim balance operations. 

The bearings in assemblies without dry-film lubrication applied were prelubricated 
at turbopump assembly with Krytox 143AZ, a LOX-compat ible fluorinated hydrocarbon 
oil. Its purpose was to prevent excessive heat buildup and cage scorching during 
the balance. A trail weight is used to calibrate the analyzer readout. A minimum 
of three runs of short duration is necessary to obtain and check the trim balance 
with the average operation requiring four to five runs of approximately 1 minute 
each. The purposes of these runs are to: 

1. Establish machine reading corresponding to the residual unbalance 
vector; its actual value is not known at this time. 

2. Determine the effect of an added known unbalance vector. Resolution 
of this resultant by vector analysis programmed on a calculator 
determines the initial unbalance magnitude and direction. An equal 
and opposite correction is then applied bv selection of a turbine bolt 
of proper weight from the nine available dash number bolts with a 
weight range of 35 to 40 grams each. 

3. Check that the resulting balance is within the unbalance allowance of 
2 g-in. 

Figure D-12 shows the HPOTP assembly mounted in the drive facility for in-bousing 
balancing operation. 
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Figure D-12. HPOTP In-Housing Balance Configuration 
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’’’here is no provision for balance correction at the preburner pump plane. The 
vibration amplitude was monitored for an extended period to confirm that no sig- 
nificant changes occur there when trim balance changes are applied to the turbine 
end. The trim balance operation, along with a tightening of component balance 
requirements has produced a reduction in observed vibration levels and in syn- 
chronous radial loads as estimated by bearing load path analysis. Trim balance 
corrections of 1 to 13 grams have been applied, with the average correction in 
the 3- to 5-gram range. 


Balance procedures preserve the LOX-clean condition of the interior of the turbo- 
pump. Only LOX-compat ible , LOX-clean materials are permitted in the operation. 
Including drive and cooling gas. 


EVALUATION OF CURRENT TRTM BALANCE 


To assess the value of tie low-speed trim balance procedure, a simulation study 
was conducted wherein the worst-case unbalance disti ibut ions were assumed to be 
corrected as required in the current practice. The wet linear response of the 
as-corrected rotors were then determined. Table D-l presents a summary of the 
conditions assumed. Figure D-l 3 indicates that the corrected HPKTP rotor response 
bearing loads are reduced sufficiently to meet the bearing load criteria. 


TABLE D-l. TRIM BALANCE SUMMARY 


TURBOPUMP 

ASSEMBLY 

CONDITION 

I 

II 

III 

IV 

■ 

WET RESPONSE 
(Figure No.) 

HPFTP 

AS ASSEMBLED 



+ 7 

+5 

+5 

D-l 3 


AFTER TRIM 
BALANCE 

H 

B 

+7 

+5 

-6 


HPOTP 

AS ASSEMBLED 

2 

5 

5 

5 


D-14, BEARING 
NO. 1 AND 2 








D-15, BEARING 
NO. 3 AND 4 


AFTER TRIM 
BALANCE 

2 

5 

-1.3 

-1.3 
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Figure D-13. HPFTP Trim Balance Effect 
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The HPOTP trim balance applies corrections to the ti ie end of the rotor only. 
The unbalance distribution identified with the worst-cas "wet" response (which 
produced the highest bearing loads) was case 5, as noted in Task I. The unbal- 
ance directions are + + + - for the four main masses, with the two turbine disk 
unbalance vectors directed oppositely (see Table D-l). As a consequence of this 
orientation, the response to the first critical speed at 10,000 rpm is small and, 
also, there would be only a small connection indicated by the trim balance proce- 
dure. For case 5, the most active response occurs at the second critical speed 
at approximately 35,000 rpm; the trim balance procedure would therefore not have 
a significant overall effect on response. 

To illustrate rotor response to a case with significant trim balance correction, 
the unbalance orientation of case 1 wa c considered. Like case 5, case 1 produces 
high bearing loads at 35,000 rpm; in addition, it also shows response peaks at 
10,000 rpm. The orientation of + + + + for the four planes would result in a 
trim balance correction of approximately 12 gram-inches directed at 180 degrees. 

The effect of in-housing trim balance for case 1 is illustrated in Fig. D-14 for 
the preburner end bearing pair, while Fig. D-15 shows the loading reduction for 
the turbine end bearing pair. The response curves of Fig. D-14 and D-15 indicate 
that the single-plane turbine end trim balance is effective in reducing first 
mode response at 10,CC-w rpm, but has little effect on the second critical speed 
at 35,000 rpm. This result is in igreement with the conclusions drawn from the 
balance simulations of Task II. To reduce the bearing loads as the rotor speed 
approaches its second critical, a correction must be made either at the main 
impeller, which is inaccessible, or at the preburner impeller. The potential 
benefit of this latter procedure is shown for dry response in the Task II Balance 
Simulation Studies. 
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APPENDIX E 


PROXIMITY PROBE CHARACTERISTICS 

Description of Several Proximity Probes and Additional Background Information . 
Fiber Optic and Capacitance Sensors 

The principle of operation of these two sensor types is sketched in Figures 
E-l and E-2. The fiber optic sensor uses two fibers (or fiber bundles) 
that are coupled by a reflective target. The coupling factor for light 
between fibers is a function of the gap, X, the fiber spacing, Y, and 
target reflectivity. The response curve offers two operating ranges; the 
more sensitive is marked in Figure E-l. 

The capacitance sensor in Figure E-2 has one insulated plate (center elec- 
trode) surrounded by a guarding electrode. The voltage across the gap, for 
a fixed alternating current into the center electrode, is an intrinsically 
linear function of the gap, X. The voltage-to-gap function is independent 
of temperature over a wide range of temperatures. 

Some of the characteristics of existing available capacitance and fiber 
optic sensor systems have been tabulated in Table E-l. 

A special version of the fiber optic sensor is shown schematically in 
Figure E-5 and its response in Figure E-6. Figure E-5 shows a fiber 
pair behind a lens. The distance between the lens and the fiber tips is 
equal to twice the focal length (F) of the lens. If a reflective target is 
placed on the other side of the lens, also at 2F, then the fibers are 
imaged sharply. If the target is moved out of focus by an amount + Xj or -X., , 
the image of the illuminating fiber becomes diffused to a diameter d^, which 
touches upon the diffuse reception circle (of diameter d f ) . If the target 
is moved further out of focus, the two circles start to overlap and the two 
fibers are coupled via the target. This is a schematic representation of 
the fiber coupling, and the signal resulting is shown in Figure E-6 (solid 
line). Due to internal reflections the actual response is as shown by the 
dashed line in Figure E-6. 
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Figure E-l. Fiber Optic Proximity Probe, Principle of 
Operation and Response Curve 
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) Linear Resoor.sa of Capacitance Probe Obtained with 
Fixed Current into “ro'oe Center Electrode 


. Cylindrical Capacitance Sensor With Guard Electrode 
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TABLE E-l 


COMPARISON OF CAPACITANCE AND FIBER OPTIC SENSOR FEATURES 
Existing Systems 



* 

Capacitance 

Sensor 

** 

Fiber Optic 
Sensor 

Measurement Range (in.) 

.040 

.014/ .055*** 

Mean Cap ( in . ) 

.020 

.010/ .090*** 

Diameter of Active Area (in.) 

.120 

0.047 

Minimum Axial Dimension (in.) 

0.020 

0.5 

Overall Probe Diameter (in.) 

.250 

0.062 

Sensitivity Stability (%FS) 

0.02 

0.2 

Noise (ZFS) 

0.05 

0.05 

Signal Transmitted by 

Coaxi-'.l Cable 

Fiber Bundle 

Cable Disconnect at Sensor 

Optional 

Not Advisable 

Cable Noise (ZFS) 

0.1 

0 

Frequency Range (kHz) 

3 

50 

Causes of Failure 
or Impaired Output 

Dirt in Gap or 
Loose Connector 

Dirt in Gap or 
Fiber Breakage 


*MTI Accumeasure System 

**MTI Fotonlc Sensor with .062 H probe 


***The larger stand-off and measuring range are available on the response curve 
to the right of the peak (see Figure E-l) which has a lower slope than the 
front of the curve (marked with operating point) . 
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Figure E-3. Flat Capacitance Sensor (Button Probe) With Guard Electrode 
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Figure E-4 . Fiber Optic Proximity Sensor 
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Figure E-5. Fiber Optic Sensor With Extender Lens 
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Figure E-6. Response of Fiber Pair With Extender Lens Near Focus 
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The major benefit of using the "extender" lens is the increase in distance 
(« 2F in Figure E-5) between the target and probe, without a loss of sensi- 
tivity.* The response slope obtained is comparable to that of the front 
slope of the sensor curve obtained without the benefit of a lens. The 
operating point on the curve in Figure E-6 can be chosen on either side of 
focus. The complexity of adding a (small) lens to a fiber bundle may be 
well worth the increased standoff. If a liquid is present in and around 
the lens and probe assembly, it will upset the functioning, the degree of 
upset depending upon the refractive index of the liquid. 

There is no published information available on operation of fiber optic 
sensors with or without lens totally immersed in cryogenic liquids. To 
avoid the immersion problem, it may be feasible to shield the optics from 
the liquid by a hermetically sealed flat window. If the front element of 
the lens is given a flat face, it can also serve as the hermetic barrier. 

It may also turn out that the refractive index difference between air and 
liquid 0 2 or H 2 is not great enough to disturb the lens function signifi- 
cantly. 


O 






The effect of LOX** on the response of a fiber optic probe is a small loss 
of peak signal (2-12Z) and a displacement of the peak (to a larger gap) by 
about 15Z. 






* Because the angle under which the fibers and target set the lens is about 
the same as that shown in Figure E-l for fibers looking directly at the 
target. In other words, in practice the fiber aperture in Figure E-l 
and the lens aperture in Figure E-5 are comparable. A lens aperture 
greater than the fiber aperture is of no benefit to sensitivity. A 
common fiber aperture is 54 (included angle) with 50%-100X attenuation 
occurring outside an included angle of 20 degrees. 

** As measured by Mr. Jim Rhodes of Marshall Space Flight Center at 
Huntsville, Alabama, and communicated to the author by telephone 
on 19 March 1982. 

RI/RD82-285 
E-8 






J 


ORIGINAL 18 

OF POOR QUALITY 


INDUCTIVE PROXIMITY PROBES 


Application Information 

I. Eddy Current Proximity Probe s 

The principle of the Inductive proximity sensor (eddy-current type) is shown 
in Figure E-7. This sensor consists of a coil powered with alternating current, 
resulting in an alternating EM field around the probe. If a conducting 
(nonmagnetic) surface is brought near the coll, the relationship between 
E and I changes. This change is the result of alternating currents (eddy 
currents) generated in the target by the external EM field of the probe. 

The eddy currents in turn generate an OI field which opposes that of the 
coil. 

The opposing effect of the eddy current's field causes the total field inside 
the target to become weaker with increasing distance from the surface. The 
depth at which the intensity is about 60X de*n from that on the surface is 
defined as the depth of penetration. 

The effect of the target is to reduce the inductance (L) and to Increase the 
AC series resistance (R) . 

1: part of the gap is filled with a conducting medium, the field of the coil 
is attenuated in the same way as by the target. If the field attenuation in 
the partially filled gap is small (< 10X), the probe can still function. This 
s of importance when using the probe in a corrosive environment where hermetic 
sealing is necessary. For that application, the probe is mounted in a hermetic 
envelope, typically of austenitic stainless steel. The thickness of the stain- 
less steel in the gap must be such that no excessive attenuation of the field 
occurs. Thus, it is important to know what factors determine the depth of 
penetration. The depth of penetration is dependent upon the electrical con- 
ductivity (p) and magnetic permeability (y) of the target material and on the 
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frequency (f) of Che AC field. These three factors are all of equal Importance 
as is shown by the formulas and graphs In Figure E-9. 

At a given operating frequency, the depth of penetration can vary two decades 
between the nonmagnetic resistive materials and the magnetic materials. Figure 3 
can be used ns a guide In selecting a usable combination of target material, 
probe envelope material and operating frequency. 

II. Magnetic-Reluctance Probe 

The principle of this prone Is Illustrated In Figure E-8. The field generated 
by the alternating current through the coil causes magnetic flux in the ferro- 
magnetic core. The flux density in the core is from two to three decades 
higher than without core. 

The flux passes from the core through the air gap back into the core. The 
air gap constitutes a major part of the magnetic reluctance of the probe. 

The typical increase in probe inductance when placing a ferromagnetic target 
in contact with the probe is So - 100Z. 

The effect of a nonmagnetic target on this probe (Figure E-8) is similar to 
that on the probe In Figure E-7. That is, a nonmagnetic target causes a 

decrease in inductance (L) and an increase in the AC series resistance of 
the probe with the core. When operated with a nonmagnetic target, the probe 
with core is an eddy-current probe. When operated with a magnetic target, 
it is a magnetic-reluctance probe. 

Applications of Inductive Probes 

In many applications, either type of probe will work satisfactorily. To judge 
where one system would be superior to another, some of the factors affecting 
probe operation are listed together with the affected areas: 
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Probe Parameter Affected Area 

Operating Frequency 1. frequency response 


Operating 

Temperature 


Stray Magnetic 
AC Field 


Operating 

Pressure 


target thlcknese 

depth of penetration ^envelope thickness 


type of medium in gap - [conductive 

(sodium, mercury)] 

[nonconductlve 
(water, air, 

steam) ] 


core 


1 . 

2 . 

3. 


curie point of ^ 

^target 

resistivity of target and envelope of transducer 
core 


permeability of^ 


1. permeability of^ 


target 

core 


target 


1. permeability modulation of; 


core 


' target 


The effects listed above are all of a magnitude that can affect probe performance 
considerably. 


When a prcbe system performs well under given circumstances, a change in a single 
operating parameter can reduce performance considerably. 

For example, the hermetic aase of a probe operating at 1 mHz is limited in face 
cover thickness to .010 inch if the cover is made out of 304 SS. The reason is 
that the depth of penetration for that probe in 304SS is about .020 inch. A cover 
thickness of .010 inch absorbs about 40Z of the external flux of the probe, which is 
an acceptable amount. Thus, when using a probe in a case (when normally it is not 
encased), the presence of the case material in the gap results in a decrease of 
probe signal and an Increase of thermal probe drift due to temperature effect on 
case material resistivity. 
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If the carrier frequency of this probe were reduced by a factor of five, the 
case material in the gap would hardly affect probe output. If a thick cover 
would be desirable because of high operating pressure, a low operating frequency 
is advantageous . 

When considering target thickness required for a given operating frequency, it 
may be advantageous to select a high operating frequency to reduce the necessary 
target thickness. 

• 

The commercially available inductance probes tend to operate at 1 mHz or above. 

The target thickness required for those is in the order of from .001 to .020 inch. 
Using a ferromagnetic target reduces the target thickness required, but introduces 
the pressure and temperature dependence of the permeability as a source of instability. 

These examples illustrate how the proper balance of operating frequency and pro- 
perties of target material and probe case materiel is of importance in probe per- 
formance. 

The high-temperature inductive proximity probes available from MTI have a magnetic 
core. Depending upon the probe type, the maximum c«.rrier frequency recommended for 
these probes is 3 kHz or 50 kHz. When used with a magnetic target they sense dis- 
placement through a considerable thickness (.050 to .100 inch) of liquid metal 
(sodium, mercury, or NaK) or austenitic stainless steel. 

The differences between a probe with core versus one without a core have a special 
significance with respect to high temperature (500*F and higher) probes. In those 
probes, the insulating materials are ceramics. The coefficient of thermal expansion 
of the ceramics and stals commonly used, differ by a factor of 2 or more. The 
result is that the coil windings on a higi temperature coil are not rigidly positioned 
by the coil form or the cement used on the windings. Therefore, individual windings 
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can move appreciable amounts (10 - 100 microinch) when the probe ia vibrated. In 
the probe without a core, the relative position of the windings with respect to the 
target (and housing) determine the output signal. In a probe with a core, the 
relative position of the core with respect to the target determines output and the 
position of the windings has only a minor effect. As it is often easier to locate 
the core firmly with respect to the probe housing than to locate the ceramic coll 
form (and/or windings), the probe with core is much less prone tc> be vibration 
sensitive. Especially when a proximity probe is used as a sensor in a high temp- 
erature load cell r pressure probe, the vibration sensitivity s of critical impor- 
tance. In such applications, the probe with core is likely to be far superior. 


RT/RD82-285 

E-13 






DDY CURRENT 


ORIGINAL PAGE IS 
CF POOR QUALITY 

3NDUCTING SURFACE (TARGET) 
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Figure E-7. Principle of Inductive Proximity Sensor 
(Eddy Current Type) 
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Figure E-8. Principle of Inductive Probe 
(Magnetic Reluctance Type) 
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PROGRAM LISTING - INTERROGATE MULTIPLE SENSORS 


• GETDAT 

• Sample 
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□ETDAT.FTN 26-MAY-82 FACE NO. 1 


SUBROUTINE GETDAT ( I BUF . NSENS » MUXA . SRATE » NPT . ARRY r IERR ) 
INTEGER*2 IBUF(NPT) 

DIMENSION MUXA < NSENS ) 

DIMENSION 4RRY(NPT. NSENS) 

C» + 

C***************** ******** **********$ ** *************** ********** 

C**»*«******«M**«M«*«*t*»*»t«*»«*t**.:.«M4.0*M. **«»«* «««*X *«*»**« .»««■* 

SUBROUTINE GETDAT — INTERROGATE MULTIPLE SENSORS 

version number ■ VOl.Ol lsneueae ■ FORTRAN 

•GETDAT* WILL SEOUENCE THE INTERROGATION OF A GROUP OF DATA 
CHANNELS AND RETURN THE ACQUIRED TIMES SERIES DATA SETS. 


cal line seauencel 

CALL GETDAT < I BUF . NSENS * MUXA , SRATE » NPT » ARRY r IERR ) 
inputs 1 

I BUF INTEGER ARRAY TO BE USED BY SAMPLING ROUTINE TO 
STORE THE ACQUIRED DATA POINTS 

NSENS NUMBER OF SENSORS TO INTERROGATE 

MUXA LIST OF A/D MULTIPLEXER ADDRESSES FOR EACH CHANNEL 
SRATE A/D SAMPLING RATE IN HERTZ TO BE USED 
NPT NUMBER OF DATA SAMPLES TO ACQUIRE FOR EACH CHANNEL 

outputs : 

ARRY ARRAY TO RECEIVE THE SAMPLED DATA 

•ARRY* MUST BE DIMENSIONED IN THE MAIN PROGRAM 
AS ’ARRY (NPT.NSENS)*, THE DATA IS CONVERTED TO VOLTS. 

IERR STATUS RETURN 

+1 = SUCCESS 

<0 » ERROR CODE AS REPORTED BY ’SAMPLE* 

lieitst ions X 

ANY LIMITATIONS OF ‘SAMPLE* APPLY TO THIS ROUTINE ALSO 

author! D. SMITH 

date! 24-MAY-82 

customer! MTI 

access! UNRESTRICTED 

external references! 
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GETDAT.FTN 


26-MAY-02 PAGE NO. 2 


C SAMPLE — HARDWARE DEPENDENT A/D CONTROL ROUTINE 

C 

C update history (cod*. date. name. reason) { 

C 

C*** ********************************* ************ *********** ***«X*«*K* 

c*» *********** ** ********************* ******** ****■»'»'♦ tr* **»»*»* ******** 

Cl- 

C 

C 

C LOOP OVER DATA CHANNELS 
C 

DO 100 I«1 » NSENS 
C 

C CALL SAMPLING ROUTINE 

C 

C SAMPLE IS A HARDWARE DEPENDENT ROUTINE THAT INTERFACES TO AND 
C CONTROLS THE OPERATION OF THE ANALOG TO DIGITAL CONVERSION 
C HARDWARE. GIVEN THE A/D MULTIPLEXER ADDRESS . SAMPLE RATE. 

C NUMBER OF DESIRED SAMPLES. AND A DATA BUFFER . 'SAMPLE* WILL 
C ACQUIRE THE SPECIFIED NUMBER OF DIGITIZED SAMPLES FROM THE 
C SPECIFIED A/D CHANNEL AT THE SPECIFIED SAMPLING RATE AND 
C RETURN THE DATA ALONG WITH A STATUS CODE TO INDICATE SUCCESS 
C OR FAILURE. 

C 

C CALULATE CLOCK PRESET - -l.EA/SRATE 
C 

KWPR— 1 .0E6/SRATE 

CALL SAMPLE (MUX A < I ) » IBUF .NPT .KWPR.O. I ERR) 

IF(IERR.LE.O) RETURN 
C 

C CONVERT THE DATA TO VOLTS 
C 

C THE A/D FORMAT IS AS FOLLOWS 2 
C 

C 12 BIT OFFSET BINARY 

C 4093 - +5 VOLTS -(1 LSB) 

C .2048 = 0 VOLTS 

C 0 »■ -5 VJLTS 

C 

DO 50 J-l.NPT 
C MASK OUT UPFER 4 BITS 

IBUF(J)-IBUF<J> .AND. *7777 
C OFFSET AND CONVERT TO VOLTS 

ARRY( J. I ) -FLOAT < IBUF <J)-2048) *3. 0/2048. 

50 CONTINUE 

100 CONTINUE 

IERR-1 
RETURN 
END 
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1 

. NL 1ST 

BIN 

1 

. NL 1ST 

LOC 

1 

• NL 1ST 

GYM 


.LIST 

TTM 


.TITLE DATA SAMPLING ROUTINE 
5 

I DTI 762 AND KUV-llA SUPPORT FOR DATA ACQUISITION 

I 

I DEFINE ADDRESS 

I 

ADC SR "177000 
ADUEC-130 
I 

KUCSR-170420 

KUVEC-440 

1 

I A/D CSR BIT DEFINITIONS 
I 

AD. GO-1 I A/D GO BIT 

AD . DMA-2 I DMA ENABLE BIT 

AD.GSO-4 I FROG . GAIN BITS 

AD.GS1-10 

AD.RTC-40 I CLOCK TRIGGER ENABLE 

AD.EXT-20 I EXT TRIGGER ENABLE 

AD.INT-100 I A/D INTERUPT ENABLE 

AD. DON-200 I A/D DONE FLAG 

AD. INC-40000 f CHANNEL INCREMENT MODE ENABLE BIT 
AD . ERR-100000 I A/D ERROR BIT 
I 

» CLOCK CSR BIT DEFINITIONS 
I 

KU .GO-1 I CLOCK GO 

KU.MDO-2 I MODE BITS 

KU.MD1-4 

KU.RT0-10 I CLOCK RATE BITS 

KU.RT1-20 

KU.RT2-40 

KU • INI- 100 80VERFL0W INTERUPT ENABLE 

KU . 0VF-200 r OVERFLOU FLAG 

KU • FOR-IOOOO 8 FLAG OVERRUN FLAG 

KU . 2G0-20000 8 ST2 GO ENABLE 

KU. IN2-40000 * ST2 INTERUPT ENABLE 

KU . ST2- 100000 8 ST2 FLAG 

I 

8 EXPANSION OF MODE BITS JSED 
I 

KU . M 1-2 8 REPEATED INTERVAL 

KU.M2-4 {EXT EVENT TIMING FROM ZERO BASE 

I 

I EXPANSION OF RATE BITS USED 
I 

KU.R1-10 II MHZ CLOCK RATE 

I 
I 

; ******** * *«>*«*** ************ ***»**»****««»» ««*****« 
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DTI 762 A/D SUPPORT AND KUV-11A SUPPORT TO ACQUIRE DATA 
CALL SAMPLE ( ICHN . I DUF » NPT > KUPR » IRCU » IERR ) 


2<R3> 

ICHN 

4(R5) 

IDUF 

6 ( RS ) 

NPT 

10(R5> 

KUPR 

12(R3> 

IRCU 

1 4 ( R5 ) 

IERR 


A/D MUX ADDRESS (0-A4) 

DATA DUFFER ADDRESS 
NUMDER OF DATA POINTS TO ACQUIRE 
CLOCK PRESET (MICRO SECONDS) 
REQUEST CONTROL UORD (NOT USED) 
ERROR UORD RETURNED 


A/D IS SUPPORTED ONLY IN DMA SINGLE CHANNEL MODE 
NO PARAMETER CHECKING IS DONE 


storage: : 


flags: 

.UORD 0 

F.ERR-1 

I A/D COORDINATION FLAGS 

error: 

1 ERROR 
1 

F . DONE “ 2 
F. MODE-4 
.UORD 0 

CODES 

TERROR UORD 

E.OK-1 

1 SUCCESS 


E.KU--8. 

1 CLOCK ERROR 


E.AD— 5 

IA/D ERROR 


E.TMO— 26. I TIMEOUT ERROR 

E.KUT— 9. 
f 

) LOCAL DEFINITIONS 
I 

CHN-2 

BUF-4 

NPT-6 

KUP-10 

RCU-I2 

ERR- 14 

{ 

I 

sample: : 


CLR 

FLAGS 

{CLEAR FLAGS 

CLR 

(?#KUCSR 

{ STOP CLOCK 

MOV 

♦E . OK .ERROR 

1 ASSUME SUCCESS 

TST 

<?#ADCSR+2 

{CLEAR A/L 

CLR 

9# ADC SR 


TSTII 

96ADCSR 


BPL 

.-4 


TST 

9#ADCSRf2 


MOV 

<*CHN(R5> »RO 

{GET MUX ADDRESS 

DIC 

•177700. RO 


SUAB 

RO 


MOV 

RO . <?♦ ADCSR 

{LOAD MUX ADDRESS 
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TSTB 0#AOCSR ICON*/ IS STARTED 

BPL .-4 i CLEAR IT 

TST (?♦ ADCSR42 


» 

» ASSUME RTC TRIGGER AND SYNC TO TACH (KU.2G0) 

I 

MOV l?KUP<RS) »t?#KUCSR + 2 » PRESET CLOCK 

BIS ♦Kfcl.Ml I KW.Rl »0#KUCSR .MODE-1 RATE-1 

BISB #AD . RTC , 04ADCSR i ENABLE RTC TRIGGERING OF A/D 

» WE ARE ASSUMING DMA MODE ONLY 


I 

MOV 

COM 

MOV 

MOV 

MOV 

MOV 

BISB 

I 

» START THE A/D 
BIS 


0NPT(R5)'RO I GET WORD COUNT 

RO 

RO r 0tADCSR + 4 » LOAD WORD COUNT 

BUF < R3 ) . 84ADCSR+6 » LOAD DMA ADDRESS 

♦DMA I NT » 04ADVEC+4 . LOAD DMA VECTOR 

♦ 340 » B4ADVEC+6 

♦ AD . DMA , 84ADCSR > ENABLE DMA 


PR0CE33 

♦KW • 2G0 » 04MJCSR i START THE CLOCK 


I 

» THE CLOCK WILL START AT THE NEXT RISING EDGE OF THE TACH INPUT 
l DATA ACQUISITION WILL BE DRIVEN BY THE CLOCK OVERFLOWS 
» DATA ACQUISITION WILL END WHEN ALL DATA POINTS ARE ACQUIRED 
I 

» NOW WE MUST WAIT FOR A/D PROCESS TO END 
t 


MOV 

R5 » - < SP ) 

JSR 

PC. ADWAIT 

MOV 

< SP ) + » R5 

MOV 

ERROR. 8ERR(R5) 

RTS 

PC 


9 

9 

I DMA INTERUPT SERVICE ROUTINE 
I 


dmaint: : 



TST 

84ADCSR 


BPL 

It 


BIS 

♦F. ERR. FLAGS 


MOV 

♦E. AD. ERROR 

it: 

CLR 

94KWCSR 


CLR 

9BADCSR 


BIS 

♦K. DONE, FLAGS 

9 

RTI 


9 

i A/D 

t 

WAIT ROUTINE 

r 

» ASSUME RT-11 

F/B OPERATING ! 


I ANY A/D ERRORS 
t SET ERROR FLAG BIT 

» SET DONE FLAG 


9 


MCALL . TWAIT 


RI/RD82-285 

F-6 


ORIGINAL PAG£ 13 
OF POOR QUALITY 

SAMPLE . MAC 26-MAY-02 PAGE NO. 


ADUAIT 

: : 





MOW 


•LCOUNT. COUNT ISET 

TIMEOUT COUNT 

AD2 : 

BIT 


• F. DONE. FLAGS 

1 DONE YETT 


BEQ 


10$ 

i EQ-NO 


BIT 


♦F. ERR. FLAGS 

i ANY ERRORS 


BEQ 


St 

»EQ-NONE 


CLR 


04KWCSR 



CLR 


etADCSR 


5*: 

1 

RTS 


PC 


1 

» WAIT 
t 

i LOOP 
t 

iot: 

CODE 

LOOP 

1 

OVER 

2 

MSEC WAITS 



MOV 


tEMTLST »R0 



• TWA IT 




DEC 


COUNT 



BGE 


AD2 



MOV 


tE.TMO. ERROR 1 TIME 

OUT ERROR 


BIS 


tF.ERRIF. DONE. FLAGS 



BR 


AD2 



I 

> TIMER ARG LIST 
I 

EMTLSTI.BYTE 0.24 
.WORD TIME 

time: .WORD 0.2 12 MSEC timer 

LC0UNT-500. ITOTAL TIME-1 SEC 

count: .word 0 

I 

.END 
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S UBROUT I NE LOP AS ( ARRY r NPT . FREQ . S I r Ni. > 

CI + 

C****************** ************************************************* ****«« 

c 

C SUBROUTINE LOPAS VERSION 1.0 

C 

C THIS SUBROUTINE WILL DESIGN A LOW-PASS FILTER OF NS GECTIONS t 
C AND RETURN THE FILTERED DATA. 

C 

C INPUTS? 

C 

C ARRY 

C FREQ 

C NPT 

C NS 

C SI 

C 

C OUTPUTS? 

C 

C ARRY 

C 

C EFFECTS? 

C 

C AUTHOR? ROBERTA M • BRUNNER DATE? 09/25/81 

C 

C OTHER SUBROUTINES REFERENCED? 

C 

C LPDES 

C 

C REVISION HISTORY? 

C 

C*** ***********************»*********************'***«******** 

c»- 
c 

DIMENSION F ( 10 » 3 ) r ARRY ( 1 ) 

DIMENSION A(10)rB{10)»C(lv) 

C 

C DESIGN THE FILTER 

C 

CALL LPDES < FREQ » SI »NS.A»B»CfGR) 

C 
C 

C WRITE OUT FILTER COEFFICIENTS 
C 

DC 35 L-lrNS 

35 URrTE(6.300)A<L) rB(L) rCCL) 

300 FORMAT < / » ' A- ' r 1PE12 . 5 . 3X t ' B* ' » E12 . 5 f 3X » 7 C- ' . E12 . 5 ) 

C 

C SET INITIAL VALUES TO 0 

C 

DO 1 N-lrNS+1 
DO ! M-lf2 
1 F(N»M)>0. 

C 

C FILTER F(T) AND GET NPT VALUES 


DATA ARRAY 

CUTOFF FREQUENCY < HZ ) 

NUMBER OF DATA POINTS 

NUMBER OF FILTER SECTIONS DESIRED <MAX«9> 
SAMPLING INTERVAL (SECS) 


FILTERED DATA ARRAY 
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DO A M-l.NPT 
F(1.3)-ARRY<M> 

GO THROUGH NS FILTER SECTIONS 
DO 2 N-l.NS 

TEHP-A<N>*<F<Nr3)+2.#F<N»2>+F<Nrl > ) 
F<N+l,3)-TEMP-B<N>*F<NUr2>-C<N>*F<N+l»l > 

CONTINUE 

UPDATE ALL PAST VALUES OF SIGNALS 

DO 3 N-l.NS+1 
DO 3 MM-1'2 
F(N»MH)»F(N.HH+1) 

CONTINUE 

SET ARRY ( H ) “PRESENT VALUE OF FILTER OUTPUT AND CONTINUE 


ARRY < H ) *F ( NS+ 1 » 3 ) 
CONTINUE 


RETURN 

END 


RI/RD82-285 



nnonnnnononnnonnnnononnnnoo 


ORIGINAL PAGE L3 
OF POOR QUALITY 


FILTER. FTN 


24-MAY-82 PACE NO. 8 


SUBROUT I NE LPDES ( FC t T i NS . A » B . C » GRAF ) 

C*********** ********** **** ****************************** $*********** 

LOUPASS BUTTERUORTH DIGITAL FILTER DESIGN SUBROUTINE 
INPUTS t 

FC - FILTER CORNER FREQUENCY IN HERTZ 

T - DATA SAMPLING INTERVAL IN SECONDS 

NS - NUMBER OF FILTER SECTIONS 

OUTPUTS l 

A . B . C - FILTER COEFFICIENTS 

GRAF - (NOT USED) 

THE DIGITAL FILTER HAS ‘NS* SECTIONS IN CASCADE. FOR THE KTH 
SECTION t THE TRANSFER FUNCTION IS S 

A(K)*<Z**2+2*Z+1 ) 

H<*>- r 

Z**2+B<K)*Z+C(K> 


THUS IF G(K) AND F < K ) ARE THE INPUT AND OUTPUT OF THE KTH SECTION 
AT TIME “> M*T t THENJ 

G(M>«A(K)*F(M)+2*r < M-l > +F < M-2 ' -B < K > *G ( M-l ) -C < K > *G < M-2 > 


C**********************************************«******************** 

DIMENSION A(l).B(l)rC(l> 

PI *3. 1415 ‘.’26536 
WCP-SIN<FC*PI*T)/COS<FC*PI*T) 

C 

DO 120 K*1 »NS 

CS»COS( FLOAT (2C(K'fNS)-l )*PI /FLOAT <4*NS>) 

X*1 . /( 1 . +WCP*WCP-2 • *UCP*CS ) 

c 

A(K)-UCP*UCP*X 
B(K)-2.*(UCF*UCP-1. )*X 
120 C<K)-(1.+UCP*WCP+2.*WCP*CS>*X 

RETURN 
END 
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THIS SUBROUTINE WILL 
THE FILTERED DATA. 

inputs: 

ARRY 

FREQ 

NPT 

NS 

SI 

outputs: 

ARRY 

effects: 


DESIGN A HIGH-PASS FILTER AND RETURN 


DATA ARRAY 

CU T OFF FREQUENCY(HZ) 

NUMBER OF DATA POINTS 

NUMBER OF FILTER SECTIONS DESIRED(MAX-9) 
SAMPLING INTERVAL (SECS) 


SUBROUTINE HIPAS ( ARRY . NPT , FREQ . S I » NS ) 

c: + 

C ****** ******************** ************************* ********* ************* 

C 

C SUBROUTINE HIPAS VERSION 1.0 

C 
C 
C 
C 
C 
C 
C 
C 
C 
C 
C 
C 
C 
C 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 

C* *************************************** ******************** 

C>- 

D I MENS I ON F < 10 » 3 ) . ARRY ( 1 ) 

DIMENSION A(10) . B( 10) rC(10) 


FILTERED DATA ARRAY 


AUTHOR: ROBERTA M. BRUNNER 


OTHER SUBROUTINES REFERENCED 1 


DATE: 09/23/81 


HPDES 


revision history: 


C 

C DESIGN THE FILTER 

C 

CALL HPDES(FRi:QrSI.NSrAf B.C>GR) 

C 

C WRITE OUT THE FILTER COEFFICIENTS 
C 

DO 35'L-lrNS 

33 WRITE ( 6 > 300 )A(L)»B(L)»C(L) 

300 F ORMAT < / » ' A* / »lPE12.5»3X» / B" , »E12.5»3Xf / C*'»E12.3) 

C 

C SET INITIAL VALUES TO 0 

C 

DO 1 N-l'NS+1 
DO 1 M-lt2 
1 F(N.M)-0. 

C 

C FILTER F ( T ) AND GET NPT VALUES 

C 


DO 4 M-l.NPT 
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C 

F< 1»3)-ARRY<H) 

GO THROUGH NS FILTER SECTIONS 
DO 2 N-lrNS 

TEMF»A<N)*<F<N.3>-2.*F<N.2>iF<N» 1 ) ) 

F<Nf 1 .3)»IEMP-B(N)*F(N+1 » 2 > -C < N ) *F < Nt 1 » 1 ) 

CONTINUE 

UPDATE ALL PAST VALUES OF SIGNALS 

DO 3 N>1'NS+1 
DO 3 HM«1»2 
F < N »Mh ) “F ( N » HH+1 ) 

CONT IN JE 

SET ARRY< H ) "PRESENT VALUE OF FILTER OUTPUT AND CONTINUE 

ARRY(M)-F(NS+1.3) 

CONTINUE 

RETURN 
END 
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SUBROUTINE HPDES < FC . T . NS » A » B r C » GRAF ) 

C*«* ************************************************ **************** 

HIPASS BUTTERWORTH DIGITAL FILTER DESIGN ROUTINE 


INPUTS! 

FC - FILTER CORNER FREQUENCY IN HERTZ 

T - DATA SAMPLING INTERVAL IN SECONDS 

NS - NUMBER OF FILTER SECTIONS 


OUTPUTS! 

A»B» C - FILTER COEFFICIENTS 

GRAF - (NOT USED) 

1 HE DIGITAL FILTER HAS *NS* SECTIONS IN CASCADE. FOR THE KTH 
SECTION » THE TRANSFER FUNCTION IS! 

A<K)*(Z**2-2*Z+1> 

H(Z>- 

Z**2+B(K)*Z+C(K) 

THUS IF G(K) AND F(K) ARE THE INPUT AND OUTPUT OF THE KTH SECTION 
AT TIME - M*T t THEN! 

G(M>-A(K)*F(M)+2*F< M-l ) +F < M-2 ) -B ( I ) *G < M-l ) -C ( K ) *G ( M-2 ) 


C ******************************************************************* 

DIMENSION A<l)'B(l)rC<l> 

PI-3. 1415926S36 
UCP-SIN(FC*PI*T)/COS(FC*PI*T) 

C 

DO 120 K-l »NS 

CS-COS ( FLOAT ( 2* ( K+NS ) -1 )*P I/FLOAT <4*NS> ) 

A(K)-1,/(1. +UCP*UCP-2 . *UCP*CS ) 

B(K>-2.*< UCP*UCP-1 . )*A(K) 

120 C < K ) - ( 1 . +WCP*WCP+2 • *UCP*CS ) *A ( K ) 

RETURN 

END 
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SUBROUTINE BANPAS : ARRY . NPT . FI » F2 » S I . NS > 

C» + 

C*)M )M4********t***«*************************U*t****Xt *******>>**»** t ******* 

C 

C SUBROUTINE BANPAS VERSION 1.0 

C 

C THIS SUBROUTINE WILL RESIGN A BAND-PASS FILTER OF NS SECTIONS. 

C AND RETURN THE FILTERED DATA. 

C INPUTS t 
C 

C ARRY 

C FI 

C F2 

C NPT 

C NS 

C SI 

C 

c outputs: 
c 

C ARRY 

C 

c effects: 
c 

c author: rokerta m. brunner date: 09/25/81 

c 

c other subroutines referenced: 
c 

C BPDES 

c 

C REVISION history: 

c 

c********************************4*************************** 

c»- 

DIMENSION F( 10. 5 ) .ARRY ( 1 ) 

DIMENSION A(10).B(10).C(10).D(10).E(10) 

C 

C DESIGN THE FILTER 

C 

CALL BPDES(F1.F2.SI.NS.A.B.C.D.E.GR) 

C 

C WRITE OUT THE FILTER COEFFICIENTS 
C 

DO 35.L-1.NS 

33 URITE( 6. 300) A(L) . B(L).CCL) »D(L) .E(L) 

300 FORMAT </ » ' A-' . 1PE12.5.2X. 'B-' .E12.5.2X. 'C»' .E12.5.2X. 

1 'D-' .E12.5.2X. 'E-' .E12.5) 

C 

C SET INITIAi. VALUES TO 0 

C 

DO 1 N-l.NS+1 
DO 1 M*1 » A 
F <N»M ) “0 • 

1 CONTINUE 

C 


ARRAY OF DATA TO BE FILTERED 
LOWER CUTOFF FREQUENCY: H.” 

UPPER CUTOFF FREQUENCY:HZ) 

NUMBER OF DATA POINTS IN THE ARRAY 
NUMBER OF FILTER SECTIONS DESIRED :MAX»9> 
SAMPLT«"G INTERVAL :SECS> 


ARRAY OF FILTERED DATA 


PkGE W 
quality 
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FILTER F ( T ) AND GET NPT VALUES 
DO 4 M-lrNPT 
F< 1 »5)-ARRY(h) 

GO THROUGH NS FILTER SECTIONS 
DO 2 N-lfNS 

TEMP-A(N>*<F(N»5)-2.*F(N»3)+F(N.l) >-B <N >*F < N+l . 4 > 

F<N+l»5)«TEhP-C<N>*F<N+l»3>-D(N>*FiN-U »2)-E(N)*F(N+l . 1 ) 
CONTINUE 

UPDATE ALL PAST VALUES OF SIGNALS 

DO 3 N-lfK'S+1 
DO 3 MM= 1 f 4 
F,'N»MM)-F(N»MM+1 ) 

CONTINUE 

SET ARRY < M ) "PRESENT VALUE OF FILTER OUTPUT AND CONTINUE 


ARRY(H)«F<NS+1 »5) 
CONTINUE 


RETURN 

END 
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C 

SUBROUTINE BPDES < FI »F2 » T » NS » A . B »C » D . E » GRAF > 

C** «*»*** t ********************************************************** 


BANDPASS BUTTERUORTH DIGITAL FILTER DESIGN SUBROUTINE 


inputs: 

FI - LOWER PASS FREQUENCY 

F2 - UPPER PASS FREQUENCY 

T - DATA SAMPLING INTERVAL IN SECONDS 

NS - NUMBER OF FILTER SECTIONS 

OUTPUTS 1 


A.B.CfD.E - FILTER COEFFICIENTS 
GRAF - (NOT USED) 


THE DIGITAL FILTER HAS *NS* SECTIONS IN CASCADE. FOR THE KTH 
SECTION r THE TRANSFER FUNCTION T SJ 

A(K)*<Zt*4-2*Z**2+l> 

H (Z>- 

Z**4+B(N)*Z**3+C<K)*Z**2+D(K>*Z+E<K) 

THUS IF GOO AND F (K ) ARE THE INPUT AND OUTPUT THE KTH SECTION 
AT TIME - M*Tt THEN « 

G(M)-A(K)*(F<M)-2*F(M-2>+F(M-4) )-B<K)*G<M-l ) -C(K)*G(M-2) 
-D<K)*G(M-3)-E<K)*G(M-4> 


C 


C 


C 


C 


«***M******t********************************t*****«*************** 

DIMENSION A(l).B(l)/C(l).D(l)rE(l> 

PI-3.1413926536 


U1-SIN(F1*PI*T)/C0S(F1*PI*T> 

U2-SIN(F2*PI*T)/COS(F2*PI*T> 

UC-W2-W1 

Q-UC*UC+2.*W1*U2 

S-U1*U1*U2*U2 

DO 150 K-l.NS 

CS-COS ( FLOAT < 2* ( K+NS ) - 1 ) *P I /FLOAT < 4#NS ) ) 
P— 2.*WC*CS 
R-P*U1*U2 
X-l • +P+Q+R+S 

A(K)-UC*UC/X 

B(K)-<-4.-2.*P+2.*R+4.*S)/X 

C(K)-(6.-2.*Q+6.*S)/X 

D<K)-(-4.+2.*P-2.*R+4.*S)/X 
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E(K)»(l.-PfQ-R+S)/X 

RETURN 

END 
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SUBROUTINE FFTR(FR. FI tK» INVRSE. IWNDW. SCALE) 
DIMENSION FR(l).FKl) 




SUBROUTINE FFTR — SINGLE PRECISION REAL FFT 


version number • UOl.Ol language » FORTRAN 

•FFTR* PERFORMS AN FFT COMPUTATION (FORWARD OR INVERSE) ON A SET 
OF DATA COMPOSED OF A REAL ARRAY < *FR* ) AND AN IMAGINARY ARRAY <*FI*> 
USING TIME DECOMPOSTION AND INPUT BIT REVERSAL. COOLEY-TUKEY 
ALGORITHM. THE ROUTINE ALSO WILL OPTIONALLY APPLY A HANNING WINDOW 
TO THE DATA SET BEFORE PERFORMING THE FFT COMPUTATION. 


M**«***t******«************»********~*******»«***tt**************)t* 


call ins seauence! 

CALL FFTR < FR . F I . K . INVRSE » IWNDW . SCALE ) 

Inputs ! 

FR ARRAY OF REAL PART OF DATA 

FI ARRAY OF IMAGINARY PAR OF DATA 

K ORDER OF TRANSFORM 

NPTS - 2*«K 

K MUST BE AN INTEGER POWER OF TWO 

INVRSE CONTROL CODE 

0 - FORWARD FFT 1 - INVERSE FFT 

IWNDW 0 - NO WINDOWING 1 - USE HANNING WINDOW 

outputs! 

THE FFT IS DONE IN PLACE. THE OUTPUTS OVERWRITE 
THE INPUT ARRAYS *FR* AND 'FI*. 

SCALE COEFFICIENT SCALING FACTOR 
limitations! 


SINGLE PRECISION ONLY 


author! DAVID A. SMITH 

data! 23-APR-82 

customer! MTI 

access ! UNRESTRICTED 


external references! 
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C NONE 

C 

C update hivtorw < cod* » data • name » reason ) » 

C 

C* ***************************************** *************************** 
C**********C** ******************************* ************************* 

Cl- 

c 

c 

PI-ATAN2<0. #-1.0) 

N-2**K 

C 

C CALCUALTE SCALE FACTOR 
C 

8C ALE-2. 0/FLOAKN) 

C 

C APPLY HANNING WINDOW 
C 

IF< IUNDU.NE. 1 ) GO TO 20 
DO 10 I-l.N 

Z-2 . 0*PI*FL0AT ( 1-1 ) /FLOAT (N ) 

FR<I>-FR<I)*.5*(1.0-C0S<Z> > 

FI(I)-FI(I)*.5*(1.0-C0S(Z) ) 

10 CONTINUE 

C 

C MUST ADJUST SCALE FACTOR FOR HANNING WINDOU 
C 

SCALE-SCALE*2.0 

C 

C DO FFT COMPUTATION 
C 

20 MR-0 

NN-N-1 

DO 2 M-l.NN 

L-N 

1 L-L/2 

IF < < MR+L ) . GT . NN ) GO TO 1 
MR-MOD < MR. L>+L 
IF (MR.LE.M) GO TO 2 
TR-FR(M+1 ) 

FR < M+ 1 >-FR<MR + l ) 

FR<MR+1 >-TR 
TI-FKM+l ) 

FI(M+1 )-FI (MR+l > 

FKMR+D-TI 

2 CONTINUE 
L-l 

3 IF(L.GE.N) RETURN 
ISTEP-2*L 

EL-L 

DO 4 M-l.L 
A»PI*FL0AT(1-M)/EL 
UR-COS < A ) 

WI-SIN(A) 

C NEGATE UI FOR INVERSE FFT CALCULATION 
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IF ( INVRSE . NE . 0 > UI — UI 
DO 4 I-h.N.ISTEP 
J*I+L 

TR-WR*FR< J)-UI*FI ( J) 

TI-UR*FI< J)+UI*FR( J) 

FR< J ) *FR ( I ) -TR 

F I < J ) »F I < I ) - 1 1 

FR < I ) *FR< I ) +TR 

FI < 1 5 »FI < I )+TI 

L-ISTEP 

GO TO 3 

END 
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SUBROUTINE RUNSUB ( XR1 . X 1 1 . XR2 . XI 2 . NPT . XR » X I ) 


SUBROUTINE RUNSUB -- SUBTRACT TWO SPECTRA 

XR1 - ARRAY OF REAL PART FOR SPECTUM 61 

XII - ARRAY OF INAG KART FOR SFECTRUN *1 

XR2 - ARRAY OF REAL PART FOR SPECTUM #2 

XI2 - ARRAY OF IMAO PART FOR SPECTRUM 62 

NPT - NUMBER OF FFT LINES 

XR - REAL PART OF RESULT 

XI - IMAO PART OF RESULT 

note: the data is the arrays must be as follows: 

1. SAME NUMBER OF FFT LINES 

2. COVER SAME FREQUENCY RANGE 

NOTE: THE SPECTRA ARE SUBTRACTED AS FOLLOWS: 

XR< I )«XR2< I ) -XR1 ( I > XI ( I >«XI2< I )-XIl < I > 


************** 

DIMENSION XR1< 1 ) »XR2< 1 ) . X 1 1 ( 1) , X 1 2 < 1 ) . XR ( 1 ) . X I ( 1 ) 

C 

DO 100 I-l.NPT 
X»XR2(I)-XRia> 

Y-XI2(I)-XI1(I) 

C WRITE(6rl000>IrXR2<I>fXI2<I>rXRl(l)rXIl<I)»X»Y 

C1000 FORMAT (I5»6F12.4> 

XR< I )»X 
XI(I)-Y 

100 CONTINUE 

RETURN 
END 
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PROGRAM LISTING - BALANCING ROUTINE 

• BALANCE 

• CMATEQ 
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PROGRAM BALANCE 

C***********************************»*****«********************************* 

C * WEIGHTED LEAST SQUARES APPROACH TO INFLUENCE COCFF. BALANCE 
C 

C * E.S.Z0R2I 

CM*«*im«ttMm**M«M*»Mt*»m«»l«»«***«*M**M**M**t**MMI4*Htt*Mt 

LOG I CAL * 1 BLK 
C 

COMPLEX ETAOt ANETAr ALPHA » Z » T » ALPHAC t AHOLD r AP rANETAC » ETA. 

1EATC .EATF . AAA . BLPHA 

ETAO - COMPLEX BASELINE DATA 

T - COMPLEX TRIAL WEIGHTS 

ANETA - COMPLEX RESPONSE DUE TO TRIAL WEIGHT 

ALPHA - IC MATRIX 

BLPHA - INVERSE OF 'ALPHA* 

COMMON/EZ1/ETAO ( 25 ) . T( 15 ) . ANETA < 13) * ALPHA! 15 . 13) . BLPHA ( 15 . 13 > 

W - WEIGHTING MATRIX 

ALPHAC - TEMP STORAGE (SEE CODE) 

AHOLD - TEMP STORAGE (SEE CODE) 

AP - TEMP STORAGE (SEE CODE) 

COMMON/E Z2/W (15.15). ALPHAC ( 1 5 . 1 5 ) * AHOLD (15. 15). AP <15.15) 

ETA - CALCULATED CORRECTION WEIGHTS 

EATC - TEMP STORAGE (SEE CODE) 

EATF - PREDICTED RESIDUAL RESPONSE AFTER BALANCING 
ETR1 - REAL PART OF OUT-OF-ROUNDNESS 

ETI 1 - IMAG PART OF OUT-OF-ROUNDNESS 

OMMON/EZ3/ET A ( 25 ) »EATC(25) 

1.EATF(25>.ETR1(25>.ETI1(25> 

IA ■ DIM OF ALPHA 
IA-15 

ZERO W ARRAY 

DO 1 J-l.IA 
DO 1 1 = 1 » I A 
W(I»J>«0. 

PRINT PROGRAM ID 

t 

OPEN < UNIT-2 .NAME- 'Til' ) 

OPEN (UNIT-3. NAME- 'TI \ ' ) 

OPEN( UNI T-6. NAME-' Til ' ) 

WRITE ( 3 . 899) 

9 FORMAT (//.10X.' *PROGRAM BALANCE*' »/. 

1' WEIGHTED LEAST SQUARES BALANCING PROGRAM' »//> 

GET SPEEDS . PROBES . PLANES 


M - ♦ OF SPEEDS 
N - • OF PROBES 
K - ♦ OF PLANES 
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URITE ( 3. 900 ) 

900 FORMAT (5X.' ENTER NUMBER OF BALANCE SPEEDS: ' *•> 

RE AD ( 2 r 903 ) M 
URITE < 3 » 901 ) 

901 FORMAT (5X*' ENTER NUMBER OF BALANCE PLANES I '»*) 

RE ADC 2.903 ) K 
WR I TE ( 3 . 902 ) 

902 FORMAT ( 5X » ' ENTER NUMBER OF PROBES: '.*> 

RE AD ( 2 » 903 ) N 

903 FORMAT (110) 

C 

C GET OUT OF ROUNDNESS DATA 
C 

URITE(3.910) 

910 FORMAT (//» IX 'BOUT OF ROUNDNESS DATA*' »/) 

DO 3000 J-l.N 
WRITE(3.905) J 

90S FORMAT <5X» 'PROBE' * 12* ' ( RADIUS . PHASE ) J '*♦> 

READ (2. 20) RAD. PHI 
20 FORMAT (2F20.0) 

PHI “PHI *3. 1415926/180. 

ETR1( J)-RAD*COS(PHI ) 

ETI 1 ( J >»RAD*SIN( PHI ) 

3000 CONTINUE 
C 

C GET BASELINE DATA 
C 

WRITE (3. 907) 

907 FORMAT (//.IX. '*BASELINE DATA*') 

DO 100 I-l.M 
WRITE! 3*904)1 

904 FOR' AT(/.3X» '**SPEED NUMBER '* 13 .'**' ) 

DO ICO J-l.N 

URITE ( 3 » 12 ) J 

12 F0RMAT(5X* 'PROBE' *13* ' (RADIUS .PHAGE) '.*) 

READ <2. 22) RAD. PHI 
22 FORMAT ( 3F20 . 0 ) 

PHI=PHI*3. 1415926/180. 

NN-((I-1)*N)+.J 
ETR* ( RAD*COS ( PHI ) )-ETRl(J) 

ETI-(RAD*SIN(PHI) )-ETIl(J> 

ETAO(NN)-CMPLX(ETR.ETI) 

UW-1.0 
W(NN.NN)»UU 
100 CONTINUE 
C 

C READ DATA WITH TRIAL WEIGHTS (T) AT POSITION (K> FOR (M) SPEEDS 
C 

URITEC3.921 ) 

921 FORMAT ( // * IX * ' *TRI AL WEIGHT DATA*') 

C 

C LOOP GVER PLANES 
C 

DO 200 IJ-l.K 
WR I TE ( 3 . 922 ) iJ 
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922 FORMAT (/. 3X . ' **PLANE NUMBER ' . 1 3 » ' »« ' ) 

WRITE ( 3*923) 

923 F0RMAT(5X» 'ENTER TRIAL WEIGHT ( MAGNITUDE * PHASE ) J ' »«) 
READ(2.20)RADT,ANGT 

ANGT-ANGT*3. 141926/180. 

TR-RADT*COS(ANGT) 

TI-RADT*SIN(ANGT) 

T < I J ) -CMPLX C T R » T I ) 

C 

C GET VIBRATION DATA 
C 

C LOOP OVER SPEEDS 
C 

DO 200 I-l.M 
WRITE (3*925) I 

925 F0RMAT(/.5X» 'ENTER VIBRATION DATA (SPEED NUMBER 13 .')' > 

DO 200 J-l.N 

UR I TE ( 3*924 ) J 

924 F0RMAT(8X. 'PROBE' .12. ' (RADIUS. PHASE ) S '.•> 

II-((I-1>*N)4J 

READ (2.20)RAD. PHI 
PHI -PHI *3.1415926/1 80. 

ANR-(RAD*C0S(PH1 > )-ETRl ( J) 

ANI-(RAD*SIN( PHI ) >-ETIl( J> 

ANETA( II) -CMPLX ( ANR . AN I > 

ALPHA ( 1 1 . 1 J ) - ( ANETA (II) -ETAO ( 1 1 ) ) /T ( I J ) 

200 CONTINUE 
C 

C PRINT OUT IC VALUES 
C 

WRITE (3. 926) 

926 FORMAT ( //» ' * INFLUENCE COEFFICIENT MATRIX*') 

DO 310 I-l.M 

URITE(3.927> I 

927 F0RMAT(/.5X. 'SPEED NUMBER' . 13./) 

DO 300 J-l.N 

II-( I-1)*N+J 
BLK-32 

WRITE (3. 928) ( BLK . ALPHA ( I I . I J ) . I J-l . K ) 

928 FORMAT (/.IX. Alt' ('.1PE10.3.'»'.E10.3.'>'»2(1X.A1.'('.E10.3»'»'. 

1 E10.3. ' ) ' ) ) 

300 CONTINUE 

310 CONTINUE 

C 

C CALCULATE BALANCE SOLUTION 
C 

c 

C - INVC TRAN( CON JG ( A ) )«U*A3 * TRAN ( CON JG ( A ) *U * ETAO 

C 

NM-N*M 

C 

C * CALCULATION OF ALPHA ( CON JG ) TRANSPOSE- ALPHAC 
C 

510 DO 520 II-l.K 
DO 520 IJ-l.NM 
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ALPHAC ( 1 1 > I J ) *CON JG< ALPHA ( IJtII) > 

320 CONTINUE 

DO 530 II-lrK 
C 

C * CALCULATION OF ALPHAC*W» AHOLD 
C 

DO 530 I-’-l.NM 
AHOLD < 1 1 * I J ) *0 . 

DO 530 KK-l.NM 

AHOLD ( 1 1 f I J ) *AHOLD < 1 1 » I J ) + ALPHAC ( 1 1 » KK ) *W < KK v I J ) 

530 CONTINUE 
C 

C * CALCULATION OF AP* AHOLD* ALPHA 
C 

DO 540 II-l.K 
DO 540 JJ-1»K 
AP < I * f J J ) "0 . 

DO 540 KK»1 » Nh 

AP < II t JJ) “AP ( II > JJ ) f AHOLD ( 1 1 r KK ) *ALPHA<KK t JJ) 

540 CONTINUE 
C 

C GET COMPLEX INVERSE OF ALPHA 
C 

DO 400 IR-1'K 
DO 401 IC-lfK 
401 BLPHA ( IR » I C ) *0 . 

400 BLPHA ( IR»IR)*1. 

CALL CMATEQ<AP, BLPHA. DETER . K » K » I A » 0 . 1 . E-l 0 . IERR . 0 > 
IF(IERR.LE.O) WRITE <3»914 ) 

914 F ORMAT < / » ' MATRIX IS SINGULAR' > 

IF ( IERR.LE.O) STOP 
C 

C BLPHA CONTAINS INVERSE 

C MULTIPLY BLPHA*AHOLD» RESULT IN AP 

C 

DO 570 II-1'K 
DO 370 JJ-lfNM 
AP ( 1 1 » JJ ) *0 . 

DO 370 KK-lrK 

570 AP< 1 1 » JJ )*AP< 1 1 » JJ )+ BLPHA ( 1 1 > KK ) *AHOLD(KKf JJ ) 

C 

C CALC CORRECTION HEIGHTS » -AP*EATO 
C 

DO 700 II-lfK 
ETA (II) *0 • 

DO 700 JJ-1»NM 

700 ETA<II)"ETA<II)+< - AP < 1 1 » JJ ) *ET A0( JJ ) ) 

C 

C OUTPUT CORRECTION WEIGHTS 
C 

WRITE<3»930) 

950 FORMAT <//' *CALCULATED VALUES OF CORRECTION WEIGHTS*'//) 
DO 850 1*1 »K 

ANG*57. 29378<ATAN2<AIMAG<ETA< I ) > » REAL < ETA < I ) ) ) 
ANG*AMOD< ANG+720. »3A0. ) 
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VALUE-CABS(ETA(I>> 

WRITE (3.932) I. VALUE. ANG 

932 FORMAT ( ' PLANE NUMBER ='»I3»' WEIGHT - '»F10.3, 
1' ANGLE - ' F7 . 1 ) 

830 CONTINUE 
C 

C CALCULATION OF CORRECTED (EATC). AND FIND (ETAF > RESPONSE 
C 

DO 980 11=1, Nh 
EATC ( 1 1 ) *0 . 

DO 881 JJ=1»K 

EATC(II)=EATC(II )+ALPHA< II»JJ)*ETA(JJ) 

881 CONTINUE 

EATF (II )=ETAO( 1 1 )+EATC< II ) 

880 CONTINUE 

C 

C PRINT EXPECTED RESIDUAL VIBRATION 
C 

WRITE<3. 1010) 

1010 FORMAT ( //, IX ' *PREDICTED RESPONSE AFTER BALANCING*') 

DO 1900 1=1, M 
WRITE< 3,913) I 

913 FORMAT ( / , 3X , ' ** SPEED NUMBER' , 13, '**' ,/,/, 

1 ' PROBE',' AMPLITUDE',' PHASE') 

DO 1900 J-l.N 
NN«<(I-1)*N)+J 
C 

r. ADD IN RUNOUT 
C 

EATF < NN)=EATF ( NN) +CMPLX(ETR1 ( J) ,ETI1(J) ) 

1820 PS I =37 . 29578*ATAN2( AIHAG (EATF (NN ) ) , REAL < EATFtNN ) > > 
PSI=AM0D(PSI+720' ,360. ) 

1860 RAD=CABS(EATF (NN) ) 

WRITE(3.913) J.RAD.PSI 
913 FORMAT (I6,F10.3,F10.1) 

1900 CONTINUE 
RETURN 
END 
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SUBROUTINE CMATEQ ( A, B. D. III. JJJ. ID. NPRT » EPS . IERR r LUN ) 


COMPLEX 


A(ID.l). B< ID. 1 ) 

inputs: 

A 

_ 

INPUT MATRIX C A(III.III) 3 

B 

- 

SOLUTION MATRIX C B(III.JJJ) 3 

D 

- 

DETERMINANT (NOT USED) 

III 

- 

SIZE OF MATRIX A 

JJJ 

- 

CLOUMNS IN MATRIX B 

ID 

- 

DIMENSIONED SIZE OF MATRICES 

NPRT 

- 

DEBUG PRINT FLAG (O-NO PRINT 1-PRINT) 

EPS 


MINIMUM VALUE FOR PIVOT ELEMENT 

(ANY PIVOT .IE. ’EPS* — > SINGULAR MATRIX) 

IERR 

- 

ERROR RETURN (1-OK -1-SINGULAR MATRIX) 

LUN 

- 

NOT USED (FOR CALL COMPATIBILITY WITH OUT OF 
CORE SOLUTION ROUTINES) 

COMPLEX 


S» D 


WHEN CMATEQ IS USED AS A COMPLEX MATRIX INVERSION ROUTINE 
MAT B MUST BE A N X N IDENTITY MATRIX 

NOTE: DEBUG PRINT GOES TO UNIT • 6 


IERR-0 

IF ( NPRT .NE. 1 > GO TO 88 
WRITE (6. 621 > 

621 FORMAT(//. 5X. 'START CMATEQ 6X BEFORE INVERSION' ) 
DO 91 I - 1. Ill 

91 WRITE ( 6 . 622 ) ( A(I.J). J-l. Ill ) 

622 FORMA T(2X. 10G12.3) 


88 CONTINUE 

NV-IABSt JJJ) 

NK-III 

FOR SPECIAL CASE OF A 1X1 SYSTEM 

IF < III ,GT. 1 ) GO TO 1100 
IF(NV.LE.O) RETURN 
DO 5 I-l.NV 
E>(1.I)-B(1.I)/A(1.1> 

CONTINUE 
IERR-1 
RETURN 

GENERAL CASE 

1100 CONTINUE 
C D » < 1 .0.0.0) 

C IF < JJJ . LT • 0 ) D-(0. 0.0.0) 
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C 

C LOOP OVER COLUMNS OP A<I.J> 

C 

KKH-KK-1 
DO 90I-1.KKM 
C 

C FIND LARGEST ELEMENT OF COLUMN • I 
C 

S-<0. O.O.O) 

DO 10J-I.KK 

R- CAPS (A(J.I) ) 

IF < R ,LT. CABS (S) > GO TO 10 

S-R 

L- J 

10 CONTINUE 

C 

C IF ALREADY ON DIAGONAL . DONE 
C 

IF ( L .EQ. I ) GO TO SO 
C 

C ELSE SWAP ROUS 
C 

DO 20 J-I.KK 
S-A(I.J) 

A< I. J)«A<L. J) 

20 A<L»J)-S 

IF (NU. LE.O) GO TO 40 
C 

C SAME FOR B MATRIX 
C 

DO 30J-1.NV 

S-B(I.J) 

B( I. J)-B(L. J) 

30 B(L.J)-S 

40 CONTINUE 

C D— D 

C 

C IF A(I.IXEPS . THEN A IS SINGULAR 
C 

50 


C 

C ELSE 
C 

S3 


60 


70 
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IF < CABS (A(I.D) . GT . EPS ) GO TO 53 
IERR — 1 
RETURN 


IPO-141 
DO 80J-IP0.KK 

IF ( CABS (A(J.I)) .EQ. 0.0 >G0 TO 80 
S-A( J. I )/A( I . I ) 

A ( J . I )■ ( O.OfO.O ) 

DO 60 K-IPO.KK 

A ( J.K ) -A (J»K)-A(I»K)*S 
IF (NV.LE.O) GO TO 80 
DO 70 K-l.NV 

B< J.K)-B< J.K)-B(I.K)*S 
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80 CONTINUE 

90 CONTINUE 

CAN'T DO DETERMINANT SINCE OVERFLOW WILL RESULT 

DO 100I-1»KK 
100 D-D*A(IfI) 

IF (NV.LE.O) 00 TO 130 
KMO-KK-l 
DO 120K-1.NV 

B ( KK . K ) • B ( KK . K ) /A < KK . KK ) 

DO 1201*1 » KMO 
N-KK-I 

DO 1 10J*N. KMO 

10 B<N»K)*B<N»K>-A<N. J+l)*B< J+1»K> 

20 B(N»K)»B<N»K)/A<N»N) 


DEBUG PRINT AFTER INVERSION 

IF ( NPRT .NE. 1 ) GO TO 99 
WRITE ( 6 • 623 ) 

623 F0RMAT</.5X» 'AFTER INVERSION' ) 
DO 92 I - 1. Ill 
92 WRI TE ( 6 . 622 ) ( B(I.J). J-l.JJJ ) 
99 CONTINUE 
130 CONTINUE 
IERR-1 
RETURN 
END 
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